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Abstract 
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We review theoretical approaches to explore the phase diagram of nuclear and quark matter at 
high baryon density. We first look over the basic properties of quantum chromodynamics (QCD) 
and address how to describe various states of QCD matter. In our discussions on nuclear matter 
we cover the relativistic mean-field model, the chiral perturbation theory, and the approximation 
based on the large-iV c limit where N c is the number of colors. We then explain the liquid-gas 
phase transition and the inhomogeneous meson condensation in nuclear matter with emphasis put 
on the relevance to quark matter. We commence the next part focused on quark matter with 
the bootstrap model and the Hagedorn temperature. Then we turn to properties associated with 
chiral symmetry and exposit theoretical descriptions of the chiral phase transition. There emerge 
some quark-matter counterparts of phenomena seen in nuclear matter such as the liquid-gas phase 
transition and the inhomogeneous structure of the chiral condensate. The third regime that is 
^T) being recognized recently is what is called quarkyonic matter, which has both aspects of nuclear 

and quark matter. We closely elucidate the basic idea of quarkyonic matter in the large- iV c limit 
J> and its physics implications. Finally, we discuss some experimental indications for the QCD phase 

diagram and close the review with outlooks. 
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1 Introduction 



The existence of our world as it stands today relies on peculiar properties of nuclei and ultimately the 
dynamics of quarks and gluons in quantum chromodynamics (QCD) at the microscopic level. Research 
on nuclear and quark matter at high baryon (quark) density is expected to anchor our empirical un- 
derstanding of the origin of matter in the Universe to a more fundamental language. Without strong 
medium effects quarks and gluons are confined inside of protons and neutrons (or hadrons in general) 
and this property is generally called "color confinement," which is a consequence of non-perturbative and 
non-linear dynamics of QCD. Another important feature of QCD is the generation of dynamical mass 
due to a condensate of quarks and anti-quarks, i.e. the chiral condensate. Theoretical understanding of 
color confinement and mass generation is one of the unanswered challenges in modern physics. 

In extreme environments such as the high temperature T, the high baryon density p, strong external 
fields, etc, the color confinement and/or the dynamical mass may be lost and a new state of matter 
out of quarks and gluons, namely, the quark-gluon plasma (QGP) could be formed. The relativistic 
heavy-ion collision experiments have aimed to create QGP in the laboratory and it is almost doubtless 
that Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory discovered QGP at 
high enough T and Large Hadron Collider (LHC) at CERN confirmed it at higher energy. With the axis 
of the baryon chemical potential /is in addition to T, one can draw the QCD phase diagram in the Hb~T 
plane. The beam-energy-scan of the heavy-ion collision is expected to explore the QCD phase diagram 
experimentally. The Facility for Antiproton and Ion Research (FAIR) at GSI and the Nuclotron-based 
Ion Collider Facility (NICA) at JINR are under construction to investigate the baryon-rich state of 
QCD matter as well as RHIC at lower collision energies. The theoretical understanding of the phase 
diagram, on the other hand, appears stalled. The purpose of this review is not to cover as many topics 
as possible on the whole phase diagram but to assemble theoretical discussions in the high density 
region accessible by the experiment. 

The study of QGP intrinsic properties is an interesting subject on its own, and furthermore, we 
could have inferred a deeper insight to confinement and mass generation mechanisms by approaching 
them not only from the vacuum but from the QGP side. From the theoretical point of view, it is quite 
non-trivial how to address "confinement" in the language of the quantum field theory, while it is possible 
to formulate "quark confinement" unambiguously in finite-T QCD with quarks made infinitely heavy 
(quench limit). This implies that the "order parameter" of quark confinement cannot have a strict 
meaning in the presence of light quarks and that of "gluon confinement" cannot be given in a simple 
way. In a gluonic medium gluons are screened by themselves which makes the meaning of confinement 
blurred. Mathematically speaking, quark confinement is well-defined only for a static color-charge in the 
fundamental representation, but gluons are associated with a color-charge in the adjoint representation. 
Indeed, it has been a long standing problem how to find a rigorous characterization of confinement and 
deconfinement for arbitrary systems with dynamical quarks and gluons. The strict order parameter for 
color confinement is, if any, still unknown, and it could be even conceivable to interpret the lack of 
the order parameter as indicating that confinement and deconfinement are connected. As we will see 
later, such a connection through smooth crossover between confinement and deconfinement may make 
it possible to approximate the QCD thermodynamics near crossover in terms of the hadronic degrees 
of freedom alone or of the quasi-particles of quarks and gluons 

As a matter of fact, it is still a challenging problem to extract analytical information directly from 
theory in order to investigate QCD matter even at asymptotically high temperature and/or baryon 
density. The running coupling constant, a s = g 2 /(Air), of the strong interaction becomes smaller at 
high energies owing to the asymptotic freedom, and this seems to suggest that confinement might be 
lost at high enough T or /ib- Then quark matter could have been realized in high-density environments 
such as inner cores of the neutron star [TJ [2]. For the purpose to investigate QCD matter at high T and 
/ib perturbative methods have been developed (see Refs. [3j H] for textbooks). It has been understood 
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by now, however, that the perturbative expansion of QCD thermodynamics breaks down badly at the 
ultra-soft magnetic scale ~ 0(g 2 T). One can immediately pin down the source of this incapability. Only 
the Matsubara zero-mode is dominant in the high-T limit and the dimensional reduction occurs [5], so 
that hot QCD is translated into magnetostatic QCD (MQCD), i.e. three-dimensional QCD with the 
magnetic coupling constant, = g 2 T. Because such three-dimensional QCD is a confining theory, 
non-perturbative information should be required for the QCD thermodynamics even in the high-T 
limit. Linde's problem of infrared catastrophe [HIE] is the most typical manifestation of the breakdown 
of the perturbation theory, and the magnetic screening mass, mu ~ 9 2 T, should be generated non- 
perturbatively from MQCD (see Ref. [8] for a concrete evaluation of the magnetic mass in a confining 
model). The introduction of /ib would not make this situation better. The magnetic mass is as vanishing 
perturbatively, which is seen concisely in the form of the hard-dense-loop (HDL) effective action [9]. 

We note that the perturbation theory can be successful in a color-superconducting (CSC) phase [TQl 
ITT] in which the unscreened transverse gluon enhances the gap energy significantly [12] and also the gap 
energy lead to interesting confining properties [T3]. Although CSC has been an important element of 
the QCD phase diagram since two seminal papers appeared [14, 15J, we will restrict the density region 
considered in this review before CSC is turned on. This is partially because it is unlikely to detect 
CSC in the heavy-ion collision experiment and partially because the phase structure involving CSC 
(including inhomogeneous states) is sensitive to model uncertainties at intermediate density and it is 
difficult to identify the robust part from such model dependent results. Interested readers can consult 
Ref. [TH] for general features of CSC, Refs. [T71 [TS] for CSC in the context of the phase diagram, and 
Ref. [19| for inhomogeneous CSC states. 

Although QCD has been a well-established theory, there is no reliable way to obtain any information 
on high-density QCD directly. The most powerful non-perturbative method, i.e. numerical Monte-Carlo 
simulation on the lattice, is of no practical use for the system at low-T and high-/XB because of the 
notorious sign problem. The Dirac determinant takes a complex value at /xb and neither its real 
nor imaginary part is positive semi-definite. The Monte-Carlo simulation based on the importance 
sampling therefore breaks down (see Refs. [201 EES] for various ideas to evade the sign problem and 
obstacles). There are many QCD-like models designed to mimic some part of QCD dynamics but 
we should keep in mind that each model has its validity limit. Nevertheless, it should be feasible to 
reach a consistent picture by collecting various pieces of knowledge from different approaches. This is 
actually, if not the best, the only possible strategy to construct the most presumable scenario for the 
QCD phase structure, unless the sign problem will be resolved. In this sense, we believe, it should 
be of paramount importance to gather known results widely from nuclear matter, quark matter, and 
also a new paradigm in between, namely, quarkyonic matter [21]. Because the notion of quarkyonic 
matter is relatively new as compared to traditional nuclear and quark matter, we shall pay attention 
to its definition and formulate it in a form of the "McLerran-Pisarski conjecture" aiming to clarify its 
confusing and sometimes misunderstood interpretation. 

To make this review article as self-contained as possible, we begin with basics of QCD in Sec. [2] and 
explain its global symmetries. Some of QCD symmetries are spontaneously broken or restored depending 
on external parameters such as T and /iB, which determines the location of the phase boundary on the 
Hb~T plane. In the theoretical research on the QCD phase diagram, center symmetry in the gauge 
sector and chiral symmetry in the flavor sector play the most important role. In particular, if the 
phase transition is of second order, the classification according to the universality is at excellent work 
to make model-independent statements for the physical properties in the vicinity of the critical point. 
The universality argument predicts that, when the temperature is raised, the deconfinement phase 
transition is expected to be of first order if the quark mass is infinitely heavy, while the chiral phase 
transition should be of second order (first order) with two (three, respectively) massless flavors. In 
reality the quark masses are non-zero and neither center symmetry nor chiral symmetry is exact. We 
overview some highlights of the recent lattice-QCD results to find that QCD undergoes crossover for 
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deconfinement and chiral restoration nearly simultaneously. 

We then proceed to the discussions on nuclear matter in Sec. |3j For the instructive purpose we 
revisit the calculations using a relativistic mean-field model, where most of the technical procedures are 
analogous to the mean-field treatment of quark matter. In fact, contemporary challenges in the quark 
matter research, i.e. the possibilities of the QCD critical point and the inhomogeneous condensates, were 
studied long ago in the context of nuclear matter, and it should be useful to flash back those theoretical 
considerations. We note, at the same time, that nuclear matter is a more complicated environment 
than quark matter; for example, it is not obvious at all how chiral symmetry should be assigned in 
terms of hadronic degrees of freedom. In other words, it is an unanswered question what kind of mesons 
and baryons would become lighter if chiral symmetry is (partially) restored. We elaborate theoretical 
background and some possible scenarios. In the final part of Sec. |3]we introduce an approximation based 
on infinitely large number of colors (N c — > oo), which is a crucial ingredient for the correct understanding 
of quarkyonic matter. In this context we shall follow holographic nuclear physics according to Ref. [22] 
and explain the most promising holographic QCD model called the Sakai-Sugimoto model [23J. In 
this review we will take a close look at the derivation of the phase diagram from this model since the 
Sakai-Sugimoto model is a genuine QCD dual having the correct physical degrees of freedom and the 
chiral symmetry. Besides, this holographic QCD model is so unique that it can describe baryonic and 
quark matter within the single framework. Any other (and more conventional) chiral models are not 
quite successful in incorporating baryonic and quark matter in a unified way. Section 3.5.3 may look a 
bit technical with many equations since we have tried not to skip details too much for convenience for 
readers who would intend to go beyond a sketchy knowledge. 

In Sec. [4] we consider the phase diagram beyond the territory of nuclear physics. The very first 
prototype of the QCD phase diagram was drawn by Cabibbo and Parisi in 1975 based on the physical 
interpretation of the Hagedorn limiting temperature [21]. Since the exponentially growing spectrum of 
hadrons plays an important role in the interpretation of experimental data and the nature of decon- 
finement crossover, we take a quick look at the statistical bootstrap model and the derivation of the 
Hagedorn spectrum. Then, we discuss the fate of chiral symmetry adopting well-used chiral models 
- the Nambu-Jona-Lasinio (NJL) model and the quark-meson (QM) model. These chiral models are 
designed to give rise to the spontaneous breaking of chiral symmetry and can properly reproduce critical 
phenomena if the chiral phase transition is of second order. Furthermore, the effect of the deconfinement 
crossover can be partially taken into account in the extended version of these chiral models. We next 
elucidate the underlying physics of the liquid-gas transition and the critical point in a parallel way to 
the previous section on nuclear matter. Also, we develop a qualitative argument on the mechanism to 
induce inhomogeneous condensates such as the chiral spiral configurations. 

The last part of Sec. 4J3 is devoted to the clarification of the characterization and the properties 
of quarkyonic matter [2TJ. To understand quarkyonic matter correctly, we need to know how nuclear 
matter should behave in the large- N c limit. Quarkyonic matter is, in a sense, not a new state of matter 
but just a possible view of large- N c nuclear matter. There are several theoretical proposals for the 
identification of quarkyonic matter, but what really features quarkyonic matter is the strength of inter- 
baryon interactions. It is not quite straightforward to give a clear definition of quarkyonic matter in 
real QCD with N c = 3. The intuitive interpretation would lead us to a picture of quarkyonic matter 
as (baryonic) matter whose pressure is sustained mainly by quarks. This explains how a nomenclature, 
quarkyonic (= quark + baryonic), makes sense. 

In Sec. [5] we make a quick overview of what has been confirmed in the relativistic heavy-ion collision 
experiment. The so-called beam-energy scan over the QCD phase diagram is still ongoing and there 
are already many suggestive results available. We will not (and cannot) cover all of interesting data. 
In addition to the heavy-ion collision data, neutron star physics has brought useful information on the 
equation of state (EoS) of dense matter [25]; the discovery of the heaviest neutron star with almost two 
solar-mass enabled us to dig out possible EoS's. It is, however, still an open question whether quark 
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matter may exist in the neutron star. The repulsive vector interaction could make the EoS hard enough 
to be consistent with such a heavy neutron star. Interestingly this issue of the vector interaction is 
closely related to the question of the QCD critical point. 

We close this review with outlooks in Section [6j In theory, as long as we cannot be equipped with 
any versatile tool, we must continue making a patchwork of various approaches that complement each 
other. Recent theoretical works indeed suggest that the baryon-rich state of matter may have rich 
contents than believed, but there are always counter-arguments that would favor less structures. We 
illustrate two possible scenarios with and without the first-order phase boundary and emphasize that 
the QCD phase diagram can be still non-trivial enough even without the first-order phase transition 
and the QCD critical point. 

2 Symmetries of the Strong Interaction 

The aim of this review is to discuss the phase diagram of nuclear and quark matter with solid boundaries 
associated with first-order, second-order phase transitions, or rather vague borders of smooth crossover. 
In most cases of the phase transition the manifestation of the global symmetries changes in accord with 
the state of matter. In this section we summarize the global symmetries of the strong interaction based 
on the fundamental theory. 

QCD is a non-Abelian gauge theory with N c colors and Nf flavors. It is known up to now that N c = 3 
and Nf = 6 in nature. Since only up, down, and strange quarks are relevant to the thermodynamic 
properties at the QCD energy scale, Aq CD ~ 200 MeV, we will limit ourselves to the case with two light 
and one heavy flavors, which is commonly denoted as "(2+1) flavors" in the convention. Furthermore, 
the system with (2+1) flavors is sometimes approximated by only 2 light flavors, which may look like a 
crude approximation, but can be absolutely legitimate to investigate the (pseudo) critical phenomena 
governed by the softest modes only. It is also a useful limit to take Nf = or to make all quarks 
infinitely heavy (though the Nf — > limit should be carefully taken [26]). Such a limit is often called 
the quench approximation, and utilized in the lattice-QCD simulation. In the quench limit the system 
of the pure Yang-Mills theory has no excitation of dynamical quarks but only gluon loops are allowed. 

The Lagrangian density of QCD consists of the pure Yang-Mills part, the quark part, and the 
CP-violating part, respectively, i.e. 

c = -\f; v f^ + pjyfd^ + g>fAp* - mJiM + o^f; v f^ . (i) 

Here, a refers to the adjoint color index from 1 to N% — 1 and i, j refer to the fundamental color indices 
from 1 to N c (namely, in the N c = 3 case, 1 =red, 2 =green, 3 —blue, particularly). Also, a is the flavor 
index (1 =up, 2 =down, 3 =strange, etc). The field strength tensor is F^ va = d ^A ua — d u A^ a + g f abc A b ^A c p 

and its dual is defined as F^ va = ^e^^F^ with the convention e 0123 = 1. This Lagrangian density 
involves 5 parameters for N c = 3 and Nf — 3. The latest estimate for the current quark masses is; 

m u = 2.3t° i MeV , m d = 4.8±g;| MeV , m s = 95 ± 5 MeV , (2) 

in a mass- independent subtraction scheme at a scale ~ 2 GeV [27]. The strong coupling constant, a s , 
is one of the physical constants and it runs with the energy scale. The world average value at present 
is [27], 

a s (m z ) = 0.1184(7) . (3) 

The CP-violating part in Eq. ([!]) originates from the structure of the ^-vacuum and there is no reason 
why 9 should be vanishingly small. The experiment for the neutron electric dipole moment yields an 
upper limit, \6\ < 10~ 10 [28], and 9 is so far consistent with zero. It is still a big mystery in theoretical 
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physics whether 9 ~ is just an accidental fact or a consequence from some unknown dynamics such 
as axions. Because 9 is such small, the CP-violating term does not have a phenomenological impact in 
normal circumstances. In far non-equilibrium situations, however, 9 may take a non-zero value locally 
(corresponding to pseudo-scalar condensation) and then the strong CP-violation could be possibly 
detectable in the relativistic heavy-ion collision [291 [30]. Furthermore, even though 9 itself is vanishing, 
the energy curvature with respect to 9 (i.e. the topological susceptibility) is non-zero and closely related 
to the m v > mass [SH Thus, the in- medium m v > mass could provide us with information on the 
topological structure of the QCD vacuum. 



2.1 Gauge symmetry 

The QCD Lagrangian density is invariant under the gauge transformation by construction. The gauge 
transformation changes the quark and the gluon fields, respectively, as 

i/>^Vi/>, — >■ V^A^ — j-d^jV^ , (4) 

where V G SU(iV c ) (i.e. V • — ■ V — 1 and det V = 1). The local gauge symmetry is never 
broken spontaneously owing to Elitzur's theorem |33j, while the global one can be which is realized in 
superconductivity leading to the Meissner mass for the gauge bosons as a result of the Englert-Brout- 
Higgs-Guralnik-Hagen-Kibble mechanism [34"l 1351 136] . The Meissner effect on gluons plays an important 
role in color superconductivity, but it is not the main subject of this article. 



2.1.1 Center symmetry 

A particularly important part of the gauge symmetry at finite T is center symmetry that characterizes 
deconfinement of quark degrees of freedom in a gluonic medium. In the pure Yang-Mills theory without 
quark field if), the genuine gauge symmetry is not SU(iV c ) but SU(A^ c )/Zat c . One can understand this 
from Eq. Q by choosing V = Zk = e 2mk ^ Nc x l7v cX iv c where e 2mk ^ Nc (k = 0, . . . , N c — 1) is an element of 
Zjy c (i.e. ( e 2m k/N c ^N c _ This Zk is certainly an element of SU(iV c ) and commutes with all elements 
of SU(iV c ). In other words, z k belongs to the center subgroup of SU(iV c ) that is nothing but Z Nc . 

Let us focus on the pure Yang-Mills part first and put aside the quark field ip for the moment. At 
finite T the imaginary-time direction is compact with a period (3 = 1/T. To maintain the periodicity of 
the gauge field, V should be a pseudo-periodic function with a twist at most by that is; V(x, t + (3) = 
ZkV(x, r). For a concrete example, one can choose V as 

V(t) = diag(e 2nikT/NclS , e 2nikr / N ^^ ( e -^i(N c -i)kr/N c ^ ^ ^ 

where the first (N c — 1) elements are identical and the last one is chosen to satisfy det V(r) = 1 for any 
r. This V(t) obviously belongs to SU(iV c ) and satisfies V(t + (3) = ZkYij). Likewise other forms with 
e -2-Ki{Nc-i)kr/N c p p erm utated with different diagonal components can also be the center transformation. 

The point is that the QCD Lagrangian density is invariant under general gauge transformation Q, 
but the boundary condition of the manifold on which the theory is defined may not. (One may prefer 
to adhere the general gauge transformation Q, and it is indeed possible to change the prescription of 
the imaginary-time formalism; see Ref. [37] for this example.) Center symmetry at finite T is thus the 
gauge subsymmetry with an aperiodicity by z&. 

Center symmetry controls the behavior of the temporal Wilson line or the Polyakov loop. We will 
use the following notation throughout this article; 

L{X) EE -pe^SidrMi^r) ? ^ = -L\xL(x) , $ EE (£) , (6) 

N r. 



7 



Center 

Symmetric •=> + * + f = Confinement 



Figure 1: Schematic picture of quark confinement in a gluonic medium. The color orientation is 
randomly distributed and the average is vanishing which is interpreted as the prohibition of single- 
quark excitation. 



where L(x) is called the Polyakov loop, £(x) is the traced Polyakov loop, and $ represents the Polyakov 
loop expectation value. 

For the physical interpretation, $ is usually considered to be related to a single-quark free energy 
f q (x) as $ oc exp[— (3f q (x)], but there are theoretical subtleties on this (see Refs. [38, 39J for details). 
The Polyakov loop correlation function is related to the heavy quark potential fq q {r) (in the singlet 
channel) in the following way; 

{l\ Xl )t{x 2 )) = Ce^[-/3f qq (\ Xl - x 2 \)} . (7) 

The inter-quark potential is an important measure to characterize whether the system is in the confined 
phase or in the deconfined phase. The quark confined phase should have a linearly rising potential, 
fqq{ r ) — <jr with a finite string tension a, so that the Polyakov loop correlation function decays expo- 
nentially at large separation. In the deconfined phase, on the other hand, the inter-quark potential is 
thermally screened and fq q {r — > oo) — > (const.) leading to a non- vanishing correlation function of the 
Polyakov loop. In summary, in the pure Yang-Mills theory, when we take the limit of \x\ — x 2 \ oo, 
the behavior of the correlation function should be 

(Confined Phase) (£ j {oo)£{0)) = $ = (Symmetric Phase) , 

(Deconfined Phase) (£\oo)£{0)) ^ $^0 (Broken Phase) , 

where we inferred <3> postulating the clustering decomposition, (ff(oo)£(0)) — > \(£)\ 2 - 

It is easy to make sure that $ is an order parameter for the spontaneous breaking of center symmetry, 
as labeled in Eq. (|8]). Under the general gauge transformation V, the Polyakov loop changes as L(x) — > 
V(f3)L(x)V^ (0). If the physical state bears symmetry under the gauge transformation that satisfies 
V{j3) = ZkV(0), the Polyakov loop expectation value transforms as 

$ e 2mk/N ^ , (9) 

that means that $ = is concluded. In other words, the center-symmetric state corresponds to the 
confined phase, whilst the deconfined phase is accompanied by the spontaneous breaking of center 
symmetry. 

In the definition according to the finite-temperature field theory, the Polyakov loop might look a 
mystical quantity. Intuitively, the Polyakov loop represents the screening factor of the fundamental color 
charge in a gluonic medium. Contrary to the naive picture, the confined phase is a highly disturbed 
state as exemplified in the strong coupling expansion on the lattice [40J, and there, the color does not 
have any preferred direction. Then all different Zn c sectors equally appear and the thermal weight for 
quark excitation in a gluonic medium picks up all Z Nc phase factors, that leads to ^ k e 2mk / N <= = o. The 
single-quark excitation is thus screened in the confined phase as a result of symmetric average in color 
space. In the N c = 3 case, for instance, this confining situation is illustrated in Fig. [TJ 

In the above, we have considered only the clean environment without dynamical quarks, but the 
situation becomes drastically different once light quarks are included. It is evident from Eq. (|4]) that 
center symmetry is lost because the boundary condition of quark field ij) directly reflects the aperiodicity 
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Figure 2: Schematic diagram of the RG flow in the dimensionless parameter space that characterizes 
the shape of the effective potential. The critical line defines the phase transition, on which the flow 
runs into the Wilson-Fisher fixed point. 



by Zk- Hence, the Polyakov loop $ is only an approximate order parameter for deconfinement and it 
always takes a non-zero value regardless of the state of matter. Nevertheless, $ effectively works well 
to indicate crossover behavior as a function of T as we will see in Sec. 12.31 

We should emphasize that it is a highly non-trivial statement that $ or any other operator cannot 
be an order parameter for deconfinement. Despite tremendous amount of theoretical efforts, there 
has been no way found to construct an exact order parameter for deconfinement in the presence of 
dynamical quarks, and all proposed candidates turned out to be unsuccessful (see Refs. [JU E2] for 
several examples). One logical interpretation for the lack of the order parameter seems to be the 
absence of strict distinction between confinement and deconfinement. In other words, for any T ^ in 
principle, the probability to find colored excitations in a thermal bath may not be strictly zero (though 
exponentially small like ~ e _/3 ^ q ). This may sound like a radical idea, but otherwise, one should question 
the theoretical framework of finite-T QCD itself. 



2.1.2 Deconfinement phase transition and critical phenomena 



The phase transition is described by the dynamics in terms of the order parameter and soft modes 
in general. Concerning the deconfinement phase transition in the absence of dynamical quarks, an 
effective theory for the traced Polyakov loop £(x) (or one can take higher-dimensional representations 
into account [43J) is useful. The effective theory can be generally expanded as 



p(glue) 



d d x 



+ c 2 \£\ 2 + c 4 |£| 4 + c Nc Re(£ N <) + 



(10) 



with T-dependent coefficients Cj(T). The fourth term is implied by Z^ c center symmetry. The effect of 



dynamical quarks can be implemented as discussed in Sec. 4.2. 1| The above effective action is supposed 
to capture the deconfinement phase transition of the (d+l)-dimensional pure Yang-Mills theory at finite 
temperature with d being the number of spatial dimensions. 

The order of the phase transition and the critical phenomena are characterized by the behavior of 



the effective action or the flow of the coefficients Cj in Eq. (10). Let us consider a coarse-grained action 



r| lue at some momentum scale k, which is defined by an effective action with the quantum and thermal 



fluctuations above k integrated out (see also Sec. 4.2.2 for a theoretical framework). The full effective 
action is therefore retrieved in the k — > limit. Then, naturally, the coefficients are dependent on k as 
well as T, that is, q(T; k). To take account of the rescaling process in the renormalization group (RG) 
equation, we make these coefficients dimensionless using k and denote them as q. 
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If the phase transition is of second order, with decreasing k, these dimensionless q's flow toward an 
infrared (IR) fixed point along the critical line (or hypersurface in general), as schematically illustrated 
in Fig. [2j As T changes, the initial point of the RG flow moves, and when it crosses the critical 
hypersurface, the destination of the flow drastically changes from the broken phase (with C2 — > — oo) to 
the symmetric phase (with C2 — > +00). If the initial point sits exactly on the critical hypersurface, the 
flow heads for the IR (Wilson-Fisher) fixed point that describes the second-order critical point, and the 
critical phenomena are uniquely determined by the flow pattern near the fixed point. That is, according 
to the general theory of critical phenomena, the critical behavior associated with the second-order phase 
transition (approaching T c from the ordered phase) is characterized by 

$ oc (£(x)) oc t p , X oc I"* , C oc t~ a , £ oc t' v , {£(0)e(x)) oc }_ 2+v (11) 

as a function of the reduced temperature, t = \T — T c \/T c . The critical behavior of the order parameter 
is specified by /3, the susceptibility x by 7, the specific heat C by a, and the correlation length £ by v. 
The two-point spatial correlation function (£(0)£(x)) has an anomalous dimension r\ at T = T c . These 
critical exponents are not all independent but they satisfy the scaling relations; 

(Rushbrooke scaling) a + 2(3 + 7 = 2 , 

(Fisher scaling) 7 = (2 — r\)v , (12) 

(Josephson scaling) dv = 2 — a . 

It has been established that the properties near the IR fixed point can be classified with the global 
symmetry of the theory and the spatial dimensions d and this idea of the classification is commonly 
referred to as the universality. In scalar theories in particular, d — 4 is known as the critical dimension. 
The IR fixed point moves to the origin when d — > 4, and eventually merges with a ultraviolet (UV) 
fixed point right at d = 4. When d > 4, the fixed point goes across the origin and is hidden in the 
unstable and unphysical region with 64 < 0. Then, the flow goes into the Gaussian fixed point at the 
origin that describes a free theory, so that the critical exponents take the classical (mean-field) values 
sec Tab. [I}. 

According to the universality argument, the deconfinement phase transition in (d+1) dimensions 
can be analyzed by a scalar theory sharing the same global symmetry (i.e. Zn c center symmetry) in d 
dimensions [44, 45J. The critical exponents can be deduced in a simpler theory than the Yang-Mills 
theory if the phase transition is of second order described by the same IR fixed point (Svetitsky-Yaffe 
conjecture). The key point is that this argument gives us the critical exponents, but it does not guarantee 
anything about the order of the phase transition. Table [T] summarizes the Svetitsky-Yaffe prediction 
for the critical exponents and the results of the numerical tests. 

The situation relevant to QCD of our interest is the SU(iV c ) case at d = 3. The universal behavior 
of the second-order phase transition has been confirmed in the numerical simulation of the SU(2) Yang- 
Mills theory. In the N c = 3 case, the effective scalar theory does not have a stable IR fixed point, 
which implies that the phase transition is possibly of first order. Indeed the SU(3) Yang-Mills theory 
has turned out to exhibit a first-order phase transition in the numerical simulation. The situation is 
subtle for N c > 4. One may think that the corresponding scalar theory with Z(N C ) symmetry may have 
a first-order phase transition, but this Z(iV c ) symmetry is dynamically enhanced to U(l) symmetry 
at the critical point, and then a second-order phase transition is theoretically possible as listed in the 
Svetitsky-Yaffe conjecture. It has been, however, established that the order of the phase transition is 
not second but first for the SU(6) and SU(8) Yang-Mills theories. Hence, the gauge theories at N c > 4 
seem not to fall into the Svetitsky-Yaffe universality class. References for the lattice simulations are 
listed in Tab. [Q 
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Spatial dimension 


Gauge group 


Svetitsky-Yaffe conjecture 


Numerical test 


a = 2 


U(l) 

r~7 / ~\t . ,i \ 

Z(iV > 4) 
SU(iV > 4) 


/ Turn's -1 1 a 

(Kl) 77 = 1/4 






Z(2), Z(4) 
SU(2) 


= 0.125 , 
K z/ = 1 

' 11 = 0.25 
a = ' 


(Universal Second Order) |16] 




SU(4) 


(Varying Exponents) 


(Second Order) |4"T] 




Z(3) 
SU(3) 


= (U1 , = 0.83 

T = L44 = 0.27 
a = 0.33 ' 


(Universal Second Order) |3H] 


a = 3 


U(l) 

Z(iV > 4) 
SU(iV > 4) 


p = 0.35 ^ = Q 6? 

7 = 1.32 n nQ 
' 77 = 0.03 
a =-0.01 ' 


(rirst Order) |49j 




Z(2) 
SU(2) 


/3 = 0.33 nfiQ 
K 1/ = 0.63 

' ri = 0.03 
a = 0.11 


(Universal Second Order) [50] 




Z(3) 
SU(3) 


(First Order) 


(First Order) [H] 


a > 4 


U(l) 

Z(iV 7^ 3) 
SU(iV ^ 3) 


71 r? = 
a = ' 






Z(3) 
SU(3) 


(First Order) 





Table 1: Svetitsky-Yaffe conjecture for various gauge theories and the numerical test in the lattice 
simulation. KT in d = 2 represents the Kosterlitz-Thouless transition. Table adapted from Ref. [4"5] . 



2.2 Flavor symmetry 

The quark part in the QCD Lagrangian density (JIJ has global symmetries in flavor space. Because the 
gauge part is not sensitive to the flavor structure, we can safely concentrate on the quark part in what 
follows. 

A quark field if) is decomposed into two chiral sectors, 

, _ 1 -75 , , _ 1 +75 , nQ , 

= — 2 — W , = 2 W ■ (13) 

The quark Lagrangian density is then expressed in terms of V'l.r, that reads, 

^{i'fDn - m)ip = ^rfB^ + ^r^A^r - fcmifa - ^r«^l , (14) 

where we suppressed the color and the flavor indices to simplify the notation. The covariant deriva- 
tive is defined as usual as {D^ij = d^Sij — igA a JT^ and the mass matrix takes the form of m = 
diag(m u , m d , m s , . . . ) in flavor space. 
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Here let us assume the chiral limit in which there are N{ massless flavors (m u = — m s — ■ ■ ■ — 0) . 
In the chiral limit the left- and right-handed sectors are totally disconnected from each other and any 
mixing between ip-^ and ipR is forbidden on the Lagrangian level. Therefore, the Lagrangian density is 
invariant under two independent chiral rotations by 



-» ipi = Lip L , ->■ ip' R = Rip R , 

L = exp(i0ft a ) G U(JV f ) L , R = exp(i6^t a ) e U(iV f ) R 



(15) 



where t a 's denote the su(iVf) algebras and the flavor index runs over a = 1, • • • , iV f 2 — 1. This means 
that the QCD Lagrangian density has the following symmetry, 

U(JV f ) L x U(iV f ) R ~ SU(JV f ) L x SU(iV f ) R x U(l) v x U(1) A , (16) 

apart from the discrete symmetry. (The center elements of SU(iVf) also belong to U(l), and so Z(iVf) is 
redundant in the right-hand side above; see Ref. [IS].) 

Realization of the global symmetry always leads to the existence of the conserved Nother current. 
Under the infinitesimal (coordinate dependent) transformations ijj'^ ~ (1 + iSO^x)^)^ and ~ 
(1 + i59^(x)t a )ipji, the change in the Lagrangian density should be at most the surface term, 5C = 
8l(x)d l ,j^ a (x) + 8^ l (x)d l ,j^ a (x), where 

jT = ^l7^>l, jr = ^R7^>R- (17) 

Then we can make appropriate combinations corresponding to the vector and the axial-vector transfor- 
mations, i.e. 

Qv = Q a R + Q a L, Ql = Q a K-Q a L, (18) 

where we defined R = j d?x V , L iR ^ a V'L,R- 
2.2.1 Chiral symmetry 



Among the global symmetry (16) the first part of SU(A^)l x SU(iVf) R is called chiral symmetry specif- 
ically. We postpone the discussions on the remaining part of U(1)l x U(1) r . Here we give a concise 
primer of spontaneous chiral symmetry breaking and associated low-energy theorems. 

So far, we have seen what symmetries the QCD Lagrangian density possesses. However, the chiral 
invariance is not manifest in the low-lying hadron spectra where any degenerate patterns between parity 
partners are absent. The resolution is that chiral symmetry is dynamically broken due to the strong 
interaction. The generator of the axial transformation, Q\, does not annihilate the ground state then. 
In a precise expression, some operator 0(x) exists so that 

(n\[iQl,O(x)]\Q)^0, (19) 

where \ is the physical state. The symmetry generated by Q\ is spontaneously broken in this state 



\Vl) if Eq. (19) holds, and the above non-zero expectation value is nothing but the order parameter. 



The simplest choice would be O = ipt a ip, leading to the chiral condensate as an order parameter; 

= (^r + $rA) ^ . (20) 

Since ipip is unchanged under SU(iVf)v transformation generated by Qy, the physical state is invariant 
and Qy\fl) = 0, which obeys the Vafa-Witten theorem [52]. In summary, chiral symmetry is sponta- 
neously broken by {'i/ji/j) ^ as 

SU(iV f ) L x SU(iV f ) R -> SU(JV f ) v . (21) 
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The physical manifestation of such spontaneous symmetry breaking is the presence of massless 
scalar particles, namely, the Nambu-Goldstone (NG) bosons. In QCD the pions are identified as the 
approximate NG bosons. We note that the original theory contains only fermions with spin 1/2 and 
gauge bosons with spin 1 and, therefore, the NG bosons must be composite states generated through 
non-perturbative dynamics. The number of NG bosons is specified by the dimension of the G/H 
manifold, where G = SU(iVf)L x SU(iVf) R is the original chiral symmetry whereas H = SU(iVf)v the 
unbroken vectorial symmetry; there appear as many NG bosons as dim(G/H) = Nf — 1. 

The axial current, J£ a , has a direct coupling to the NG boson Tr b , the strength of which can be 
parametrized by the decay constant f n as 

(n\J£ a (0)\n b (q))= l f w q»5 ab . (22) 

In the experimental point of view f n is a useful measure, and the empirical value is f n = (92.4±0.2) MeV. 
In a medium f n serves as an order parameter for (partial) chiral symmetry restoration. The divergence 
of the axial current is non-vanishing due to an explicit symmetry breaking by the presence of m q , i.e. 

<Vr = , (23) 

which is known as the partially conserved axial- vector current (PCAC). If the explicit breaking is 
turned off, the pion is precisely massless as dictated by the NG theorem and thus the axial current is 
conserved. Utilizing the PCAC and the soft pion theorems, one can derive low-energy theorems. The 
most well-known is the Gell-Mann-Oakes-Renner (GOR) relation, 

flml = -m q (^> , (24) 



and another important low-energy theorem is the Goldberger-Treiman relation, f^g^NN = niN gA, 
where g n NN is the pion-nucleon coupling, is the nucleon mass, and g& is the nucleon axial charge. 



2.2.2 U(1)a symmetry and the quantum anomaly 

The QCD Lagrangian density ([l| has not only chiral symmetry SU(A^)l X SU(A^)r but also U(l)v x 
U(1)a- The U(l)y symmetry corresponds to the baryon number conservation and should not be broken 
except in color-superconductivity at asymptotically high density. Then, according to the NG theorem, 
there should be one more massless boson associated with the U(1)a symmetry breaking. There is, 
however, no such light particle in the iso-singlet pseudo-scalar channel. Instead of the NG boson, the 
lightest meson in this channel is if whose mass is m v > = 958 MeV, that is too heavy to be the NG 
boson. The problem of missing U(1)a NG boson is called the U(1)a problem. 

Now it is a textbook knowledge that the U(1)a symmetry in the classical Lagrangian is explicitly 
broken via quantum effects. Such quantum anomaly appears ubiquitously in the gauge field theory. 
If we choose to adhere the vector gauge symmetry, the quantum anomaly arises in the axial current 
generally as 

^T^-^^tr^F^J , (25) 
which is non-vanishing only for a = (iso-singlet) because t a, s are traceless for a = 1, . . . , iV f 2 — 1. 



Therefore the chiral symmetry breaking pattern is truly the one given by Eq. (21 ) and rf cannot be the 
NG boson. We make a remark here that it is not such easy to disprove the existence of the massless pole 
because j£° can be modified into a conserved but gauge- variant form. Interested readers can consult 
Ref. [53]. 

From the microscopic point of view U(1)a symmetry is broken through the instanton-induced inter- 
actions. Let us take a look at the derivation with an assumption that an instanton with the size p is 
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placed at z [54J. Then, the partition function can be expressed in a form of the functional integration 
over the collective coordinates as 



d*z dp p' 1 (det D[A})- 1/2 {det G[A}) e 



(26) 



after the one-loop integration around the instanton background A(x). Here D[A] and G[A] are the 
gluon and the (massless) quark propagator inverse, respectively, in the presence of the background 
A(x). 

It is important to note that a fermionic zero-mode is accompanied by the instanton gauge configu- 
ration, which has definite chirality in accord with the sign of the winding number Qw! f° r A(x) with 
Qw — 1 for instance, there exists a wave-function that satisfies G f [y4.]'0Ro = (apart from the mass 
term). Therefore, the partition function Z should be zero in this case. Now we deform the theory 
slightly so that Z stays finite, that is, we insert a source term, ?/>l J(x)i/jr, in the fermionic part G. 
Then, at low energies, we can approximate the Dirac determinant as 



det[G + J] ~ det / g^x-^lo^Oz)^: 



^ro 



(27) 



Here det in the right-hand side is taken in flavor space. Now we shall consider an effective theory of 
QCD with all gluonic fields integrated out. Then this effective theory should be written in terms of 
quark fields ip(x) and ip(x). The question is how to incorporate the instanton effect with the quark 
degrees of freedom only, or how to find an effective interaction vertex induced by instanton. We can 



find such an instanton-induced interaction to reproduce Eqs. (26) and (27); the following vertex can 



mimic Eq. (27) and thus embody the instanton effect; 



C 



U(l)/ 



n{p, T = 0) [det ^(1 - 7 5 )^ + det ^(1 + T 



(2* 



including the latter term coming from an anti- instanton. As expected, the interaction (28) breaks 
the U(1)a symmetry explicitly, while chiral SU(iVf)L x SU(A^)r, symmetry is kept unbroken. This 
determinant form had been postulated earlier from the symmetry consideration [55] and the above 



interaction (28) is now called the Kobayashi-Maskawa-'t Hooft (KMT) interaction. 

In the above derivation it is clear that the interaction strength is proportional to the instanton 
density that is given in the one-loop perturbation theory by 



n(p,T = 0) 



57T 



P 



-5 p -8tt 2 /9 2 



(29) 



originating from det_D[*4.]. Thus, the strength of the U(1)a breaking interaction depends on how many 



instantons the system would accommodate, though the anomaly itself is never diminished and Eq. (25) 



is not altered by the finite-T effect. This fact opens an interesting possibility. That is, at finite T, the 
one-loop instanton density is modified as [56J, 



n(p, T) 



57T 



g z 



p 



-5 p -8ttV3(p) 2 



cxp 



-7T 2 p 2 T 2 



2N C N t 



(30) 



The latter exponential factor suppresses the instanton excitation significantly. Similarly, at finite den- 
sity, the exponential suppression factor has been found [57J, and after integrating over the instanton 
size, the suppression factor should be T -14 or p~ 14 . In any case these expressions are perturbative ones, 
and any quantitative estimate is not really trustworthy near T c . Nevertheless, it is quite suggestive that 
n(T ~ T c ) in the above expression is already a half of n(T = 0). Then, one may well speculate that 
the U(1)a symmetry is effectively restored before the critical temperature/density is reached [58], and 
if so, as we will see in the next subsection, the universality class and the critical properties would be 
significantly affected. 
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2.2.3 Chiral phase transition and critical phenomena 



The order of the chiral phase transition was systematically investigated first by the RG analysis [59J. 
As we explained, if the phase transition is of second order, the fixed-point structure yields the critical 
exponents not relying on the microscopic dynamics. One can also say that, reversely, the absence of 
the stable fixed point suggests that the theory should exhibit a first-order phase transition. 

To clarify the nature of the phase transition we need not to solve QCD at finite temperature but we 
can make use of an effective description with chiral symmetry. The simplest model is the linear sigma 



model composed from the meson field matrices, M, 



M -> iSMR 



~ ip-Rjipu and its conjugate M' that transform as 
M f ->■ R ] M ] L. (31) 



The general Lagrangian density consistent with chiral symmetry is decomposed into £ = £ <J +£u(i) A +£ / ' 
where 



1 /„ +s 1 o Aln /,r,/tM 2 7l2 92 



to = -tr (^M^Mt) + - / u 2 tr(MM t ) 



[tr (MM f )]' 



tr (MM f MM f ) 



(32) 



and the instanton- induced KMT interaction £u(i) A = c(det M+det M^j which is reminiscent of Eq. (28), 
and the explicit breaking part £' = /o tr (t°M) + / 3 tr (t 3 M) + / 8 tr (t 8 M) + ■ • • ranging over all the Cartan 
subalgebras. Here, in C a , the first is the kinetic term, the second is the mass term, and the third and 
the fourth are potential terms that results in the spontaneous symmetry breaking, respectively. From 
the symmetry property M can be decomposed as 



Nf- 



M 



(33) 



a=0 



with the scalar fields, £ = t a a a = |(M+ M f ), and the pseudo-scalar fields, IT = t a rr a = i(M - M f ). 
(Note that M changes to under the parity transformation.) Because M transforms as V'MV under 
a vectorial rotation (V = L = R), one can classify the matrix components in terms of the hadron 
language. In the N c = 3 case, for example, 



v/2 



V2 



a 



K*- K*° 




n 



V 



K~ 



V3' 



71 



V2 



^6 

K° 



K + 
K° 

2 ~ + :I 



(34) 



with two independent condensates a and £ (isospin symmetry is assumed) corresponding to two chiral 
condensates, (mm) = (dd) and (ss). In this case rj Q is anomalous, while there are eight NG bosons, tt°, 
7T , -ft' , X , and ?7 S , in the chiral limit. In reality SU(3)v symmetry is substantially broken by the 
strange quark mass m s , so that the masses in the strange sector, K°, K°, K ± , r)g, are lifted up by m s . 
Besides, rjg and r/o can mix together by m B , which amounts to physical states of rj and rj'. 

One could put this linear sigma model directly in a finite-T environment and investigate the phase 
transition of chiral restoration. For the present purpose to study the order of the phase transition and 
the critical properties, however, the universality argument is extremely useful. It is sufficient to begin 
with a dimensionally reduced description after the integration over temperature fluctuations, as already 



discussed in Sec. 2.1.2 , and then we can simply focus on the 3-dimensional linear sigma model with 
T-dependent couplings gi(T; k), g2(T; k), c(T; k), etc, at the RG scale k. Since we are interested in the 
critical phenomena only, we shall work in the chiral limit, /o = f% = fs = ■ ■ ■ = 0. 
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Flavors 



Pisarski-Wilczek conjecture 



N { = 2 



N { >3 



(3 = 0.38 
7 = 1.44 
a = -0.19 

(First Order) 



v = 0.73 
r) = 0.03 



N r = l 



/3 = 0.35 
7 = 1.32 
a = -0.015 



v = 0.67 
7] = 0.038 



N t > 2 



(First Order) 



Symmetry breaking pattern 



SU(iV f ) L x SU(JV f ) R SU(iV f ), 



SU(iV f ) L x SU(JV f ) R x U(1) A -> SU(iV f ) 



v 



Table 2: Pisarski-Wilczek conjecture for the chiral phase transition in d = 3 space. 



Without the axial anomaly Let us first consider the case without the U(l ^-breaking term (i.e. 
c = 0). This is possible if the effective restoration of U(1)a symmetry occurs slightly below T c due to 



instanton suppression as argued in Sec. 2.2.2 In this case the symmetry breaking pattern is 



U(iV f ) L x U(JV f ) R -> U(JV f ) v 



(35) 



and there are iV f 2 massless NG bosons. The critical phenomena can be analyzed by the e expansion, 
in which the spatial dimension is taken as d = 4 — e and e is assumed to be a small number. Because 
d = 4 is the critical dimension, the coupling constant around the fixed point is ~ 0(e), which justifies 
the perturbation theory. 



One can then carry out the one-loop calculation using the Lagrangian (32) to find the following (3 
functions; 



dgi 
dk 



-m + 



iV f 2 + 4^ 4N { 



9i + -^9192 + 9 2 



I3 2 = k 



dk 



2N { 2 

-eg 2 + 2g x g 2 + —g 2 



(36) 



for general iVf. We introduced a notation for the dimensionless couplings, gi = g\k e and g 2 = g 2 k e . 
The first term in the right-hand side of Eq. (36) involving e appears trivially from k~ e in g~\ and g 2 . One 



can solve this set of differential equations to draw the RG flow pattern. The fixed points can be found 
readily from the conditions, (3\ = (3 2 = 0. Two solutions always exist, one at (§l,g 2 ) = (0,0), and the 
other at 



m,92) 



3e 







and two more fixed-points can appear for iVf < a/3 at 



3e 



2(27 



8iV f 2 



9-iV/TA r f\/2(3-^ 



-5JV f + iV f 3 ± 3 



2(3 - Nf) 



(37) 



(3* 



From explicit analysis it turns out that Eq. (37) is an IR-stable fixed point for < iVf < \/2, and then 



this fixed point would represent the critical phenomena of the 0(2iVf) universality class (which is rather 
close to the situation with the axial anomaly). Therefore, when iVf = 1, the phase transition can be of 
second order, and if so, it is characterized by the 0(2) critical exponents. For N { > a/3 this fixed-point 
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becomes unstable in the cfe direction and there is no IR-stable fixed-point at all. Thus, when Nf > 2, it 
is impossible to have a second-order phase transition and it should be of fluctuation-induced first order. 

One might have had an impression that the leading-order e expansion is insufficient to conclude 
whether the phase transition is of first order. It should be mentioned that this conclusion of the 
fluctuation-induced first-order transition has been confirmed non-perturbatively in a more sophisticated 
framework of the functional RG method [60l 161] . 



With the axial anomaly If the instanton-induced U(l)A-breaking interaction remains non- vanishing, 
as is a more conventional scenario, the critical properties are drastically changed especially in the case 
with Nf = 2. Because, in the iV f = 2 case, the original group is SU(2) L x SU(2) R ~ 0(4), the critical 
phenomena should belong to the 0(4) universality class. This can be analyzed in the e expansion in the 
same way as previously. The instanton term det M + det yields substantial mass terms for the scalar 
particles only, and E decouples from the low-lying dynamics. Then, with only II, two terms involving 
gi and g 2 in Eq. (32) are no longer independent, and the (3 function is reduced to (3i in Eq. (36) with 
g~2 — 0, which has a stable fixed-point unlike the previous case without the anomaly. 

For Nf > 3 no IR-stable fixed point is found, and it is likely that the phase transition is of first 
order. In fact, even at the tree level, the KMT interaction at Nf = 3 involves terms like a 3 in the chiral 
limit where ( = a. The presence of such cubic term strongly implies that the phase transition should 
be of first order. 

To summarize the expected critical properties of the chiral phase transition in d = 3 space, we shall 
make a table based on this Pisarski-Wilczek conjecture [59] as in Tab. [2j 



2.3 Highlights of the lattice-QCD results 

In this subsection let us look over the numerical results from the Monte-Carlo simulation of QCD on 
the lattice. Although there are some data available for spatially 2-dimensional case, we will limit our 
discussions here to the realistic situation with d = 3 in space. For details on the lattice-QCD simulation, 
readers may consult Refs. [62, 63J. 



Pure Yang-Mills theories The left panel of Fig. |3| shows a figure adapted from Ref. [64J for the 
renormalized Polyakov loop in various representations in color space. It is the renormalization effect 
that renders $ exceed the unity. We see that the order parameter clearly shows a discontinuous jump 
at T = T c . (We should note that $ = below T c where center symmetry is unbroken which is omitted 
in Fig. |3j) In the pure Yang-Mills theory the physical scale is set with the string tension a, and the 
critical temperature is found as T c = (0.630 ± 0.005) yfa ~ 280 MeV (65]. In the right panel of Fig. [3] a 
figure adapted from Ref. [66] presents thermodynamic quantities such as the internal energy density e 
and the pressure p for the pure Yang-Mills theories with various gauge groups. From this figure we can 
understand that the thermodynamics has only weak dependence on the number of colors once properly 
normalized by the Stefan-Boltzmann value (for more general large- N c studies, see Refs. [67, 68J). It 
is also evident in the behavior of s(T) that the physical degrees of freedom rapidly change at T = T c 
corresponding to the liberation of colored excitations. On the other hand, the pressure p is always 
continuous regardless of the order of the phase transition. The state-of-the-art lattice data at high 
precision is available in Ref. [61]. The blue curves in Fig. |3j represent the potential fit results for later 



convenience, which we will discuss in Sec. 4.2.1 



Full QCD Figure [4] is the counterpart of Fig. [3] for the full QCD simulation with dynamical quarks. 
In this case not only the Polyakov loop $ but also the chiral condensate (V>V) i s an interesting quantity 
and they both show smooth crossover as seen in the left panel of Fig. [4} To conclude that they are 
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Figure 3: (Left) Renormalized Polyakov loop in various representations (with the dimension D) in the 



SU(3) Yang-Mills theory. The blue curve represents the potential fit from f2 glue [$] discussed in Sec. 4.2.1 
Figure adapted from Ref. [64J. (Right) Energy density and pressure in the SU(iV ( 
normalized by the Stefan-Boltzmann value 
again. Figure adapted from Ref. [66J. 



Yang-Mills theories 
The blue curves represent the potential fit from f2 g i uo [$] 



crossover in a rigorous way, the finite- volume scaling is necessary [70]. For the practical purpose, it is 
more appropriate to use a subtracted order parameter at T defined by [71] . 

Ai, s = ^ , M(T) = m s ((^)i - ^(« 8 )iV T 4 , (39) 

rather than the naive chiral condensate, where m s is the strange quark mass in the lattice unit and iV r 
is the site number along the temporal direction. This combination is chosen to be free from an additive 
divergence ~ m q A 2 inherent in the chiral condensate. This subtracted order parameter, Ai s , is plotted 
in the left panel of Fig. |4j 

Interestingly, as seen from the left panel of Fig. |4j two crossovers of deconfinement and chiral 
restoration happen almost simultaneously; the renormalized Polyakov loop and Ai iS start increasing and 
decreasing, respectively, in the same temperature range. One may define the pseudo-critical temperature 



(see Eq. (40)) but it is not meaningful to ask whether two pseudo-critical temperatures coincide or not 
quantitatively. As we will discuss later, the chiral crossover is pretty close to the second-order phase 
transition and critical scaling properties are expected, while the deconfinement crossover spreads over a 
wide temperature range. In this sense, it is often emphasized that only the chiral restoration is critical 
and there is no phase transition associated with deconfinement. This statement is indeed true, but we 
would insist that it is still a highly non-trivial problem how these two crossovers can be such correlated 
(see Ref. [71] for a possible scenario). In view of the internal energy density, e(T), in the right panel of 
Fig. |4| deconfinement crossover occurs in a narrower window of the temperature than suggested from $ 
in the left panel. We see that the physical degrees of freedom rapidly increase within the temperature 
range < 100 MeV and the phenomenon of color deconfinement makes reasonable sense though it does 
not come along with critical scaling. In the same way as the pure gluonic case, the pressure p(T) is a 
continuous function in the whole T range. 

Sometimes, to quantify the crossover location, the pseudo-critical temperature is defined. The 
pseudo-critical temperature should coincide with the critical temperature when crossover turns to a 
second-order phase transition in the chiral limit. One should keep in mind that the pseudo-critical 
temperature is a quantity only for convenience and it does not add any physical insight. There are 
several different prescriptions which lead to different pseudo-critical temperatures. 
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Figure 4: (Left) Renormalized Polyakov loop and the subtracted and normalized chiral condensate 



in the (2+l)-flavor QCD simulation. The definition of Ai )S is explained in Eq. (39). Figure adapted 
from Ref. [72] and originally two separate plots are superimposed with the horizontal axis properly 
adjusted. (Right) Energy density and pressure in the (2+l)-flavor QCD simulation. Figure adapted 
from Ref. [73] and originally two separate plots are superimposed with the vertical axis adjusted by the 
Stefan-Boltzmann limit. 

A frequently used prescription for the pseudo-critical temperature is to see the peak position of the 
chiral susceptibility, 

T8 2 \nZ 

X ^=V^nf> (40) 

which gives T c = 151 (3) (3) MeV for chiral restoration, where the first error comes from the T ^ 
and the second from T = analyses [75j. We note that this estimate gives the systematic error but 
does not have uncertainty from the choice of the prescription. In view of the order parameter behavior 
in Fig. |4| we should consider that the pseudo-critical temperatures from various definitions could not 
avoid substantial uncertainty more than AT C ~ ±20 MeV. 

It is obvious at first glance from Fig. [4] that there is no prominent phase transition at all associate 
with deconfinement in full QCD data, while $ serves as a good order parameter in the pure Yang-Mills 
theories. Thus, the pseudo-critical temperature for deconfinement is less meaningful than that for the 
chiral phase transition. One way to locate the deconfinement pseudo-critical temperature is to look for 
the inflexion point of $(T) or the peak position of the temperature derivative of the Polyakov loop, 
d&(T)/dT. This prescription gives an estimate T c = 176(3) (4) MeV [75]. Needless to say, the real 
uncertainty width should be much larger than the errors, that is, AT C ~ ±50 MeV at least from the 
curve in the left panel of Fig. [4] Another possible way to define the pseudo-critical point is to use the 
peak position of the Polyakov loop susceptibility in the same manner as for chiral restoration. The left 
panel of Fig. [5] shows an example at various light quark masses. We can notice a general trend that two 
pseudo-critical temperatures stick to each other, which is clearer for larger quark mass (mi = 0.4m s ) 
and becomes vague for small (physical) quark mass (mi = 0.1m s ). The latest lattice data can be found 
in Ref. [76] . 

Instead of the Polyakov loop, it is also possible to utilize the strange quark number susceptibility to 
define the pseudo-critical temperature; 

Td 2 \nZ 

Xs = v^T' (41) 



where /i s is the chemical potential for strange quarks. This quantity, Xs, behaves in a similar way to the 
Polyakov loop (see the right panel of Fig. [5]), which can be easily understood in a model as commented 
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Figure 5: (Left) Peaks in the Polyakov loop susceptibility xl and the chiral susceptibility xi, from which 
the pseudo-critical temperatures are defined for various quark masses. Figure adapted from Ref. [TTJ 
and originally three separate plots are superimposed. (Right) Strange quark number susceptibility in 
the (2+l)-flavor lattice simulation together with the renormalized Polyakov loop. The pseudo-critical 
temperature for deconfinement is defined from the inflexion point of the susceptibility with respect to 
the temperature. Figure adapted from Ref. [75] and originally two separate plots are superimposed. 



below Eq. (137). The inflexion point of Xs gives the pseudo-critical point consistent with the definition 
by means of the Polyakov loop, i.e. T c = 175(2) (4) MeV [75]. 

In the application of the lattice-QCD data to the heavy-ion collision experiment it is crucially 
important to check how close to criticality the QCD phase transition is. If u- and <i-quarks are massless, 
the chiral phase transition should be of second order that belongs to the 0(4) universality class as argued 



in Sec. 2.2.3 If the physical quark mass is situated in the critical region, the subtracted order parameter 



M in Eq. (39) should obey the critical scaling, which can be parametrized in terms of the universality 
scaling function and scaling variables. The lattice data as reported in Ref. [751 179] is consistent with the 
0(4) or 0(2) universality class. (The 0(2) universality class is also expected for the coarse staggered 
fermion.) Therefore, the pseudo-critical temperature for chiral restoration conveys a definite meaning 
if it is fixed such that the magnetic equation of state is fitted. 



Finite density Because the chiral pseudo-critical temperature is relatively well-defined, it makes 
sense to draw a phase boundary on the QCD phase diagram using this pseudo-critical temperature. 
The Monte-Carlo simulation is not applicable so far at finite density, but it is feasible to include the 
density effect by means of the Taylor expansion in terms of yU q /T. In this way one can estimate the 
/^-dependence of the pseudo-critical temperature or the curvature of the phase boundary with respect 
to /i q . That is, the pseudo-critical temperature is expanded as 



(42) 



and recent lattice-QCD results yield the estimate, n q = 0.059 based on xm m Ref- [I2]> K q = 0.059 
based on Xfip and K q = 0.080 based on Xs in Ref. [80] (see also Ref. [81] for former results). These values 
are smaller than the curvature expected from the freeze-out line of the heavy- ion collision experiment, 
but more than twice larger values are reported based on the EoS [82] • Also, there are suggestive data 



on the QCD critical point, which will be introduced later in Sec. 5.1.2 



Effective U(1)a symmetry restoration Apart from thermodynamics, there are useful data available 
from the lattice simulation, and among many, let us mention on suggestive data for effective U(1)a 
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symmetry restoration. It is not quite feasible to measure the strength of the KMT interaction or 
the instanton density directly, and instead, the (pure) topological susceptibility Xtop can detect the 
instanton fluctuations and should be connected to the KMT interaction. The lattice simulation at 
finite T has resulted in a rapid drop of Xtop around T ~ T c [83], which is a circumstantial evidence 
for the U(1)a restoration above T c . More direct analyses on the U(1)a symmetry encoded in the Dirac 
spectra are ongoing at finite T but zero p q |84j. It might be even possible to probe the in- medium 
7/ mass experimentally [85], and thus quantitative outputs from the lattice-QCD simulation should be 
more and more important in the future studies. 



3 Nuclear Matter 

In this section we make a review over fundamental properties of nuclear matter and raise some theoretical 
issues at high baryon density before quark matter emerges. Among various observables in nuclear 
matter, the saturation density and the binding energy are the most basic physical quantities. The 
binding energy is inferred from the semi-empirical mass formula, i.e. the Bethe-Weizsacker mass formula, 
M(Z, N) = Zrrip + Nm n — B(Z, N), with the binding energy parametrized as 

B(Z, N) = a v A - a s A 2 ' z - a c -^ - a sym (A + 5 A , (43) 

where A is the mass number. The first, the second, the third, and the fourth terms represent the volume 
energy, the surface energy, the Coulomb energy, and the symmetry energy, respectively. The last term 
is the pairing energy. For nuclear matter where A and the volume V are infinitely large with p = A/V 
fixed, the surface term and the pairing term are irrelevant. For symmetric nuclear matter [Z = N) the 
symmetry energy also drops off. (In this article our concerns are limited to symmetric nuclear matter 
only.) Then, the empirical value for the binding energy per nucleon in symmetric nuclear matter is 
B I A ~ 16.3 MeV, and this is realized at the saturation density, po — 0.17 nucleons/ fm 3 . The Fermi 
momentum corresponding to po is kp — 260 MeV in a free Fermi gas, which is less than one third of 
the nucleon mass — 939 MeV. Therefore, in this sense of kp/m^ < 1, normal nuclear matter can 
be regarded as a dilute system of nucleons which allows for an expansion in terms of kp. This fact is 



essential to distinguish it from quarkyonic matter as we will address in details in Sec 4.5 From this 



observation it should be legitimate to approximate the binding energy in the following expanded form; 

fl w=!=-^ +< 4- /5 i- (44) 

where the first term is the kinetic energy of a free Fermi gas and the cubic and the quartic terms 
represent the interaction effects with dimensionless parameters a and 0. In Ref. [HE] it is argued that 
a = 5.27 and = 12.22 can reproduce the binding energy and the saturation point, and also the 
nuclear matter compressibility, K = k F (d 2 B/dk F ) = 236 MeV, within a reasonable range. This simple 
parametrization works fine to describe the ground state properties of nuclear matter, but it does not 
suffice for the finite-T research on the phase transitions. It is, thus, indispensable to develop systematic 
instruments to evaluate a and (5 from microscopic dynamics. 



3.1 Relativistic mean-field model 

In the Walecka model or the cr-u model [HZ], nuclear matter is approximated by a Fermi gas of in- 
medium nucleons that feel mean-fields in the scalar channel (denoted by a) and the vector channel 
(denoted by u>). This is an old and simple model but is still useful to test theoretical ideas and derive 
an EoS even in contemporary contexts (see Ref. [HH] for example for the application to the relativistic 
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heavy-ion collision). Moreover, the calculation techniques for the mean-field treatment in the Walecka 
model are in complete parallel to those in quark matter. This analogue is quite useful not only on 
the technical level but also for the physical interpretation of the first-order phase transition and the 
associated critical point at the terminal of the first-order phase boundary. 



3.1.1 Saturation properties 

The bare nucleon-nucleon (NN) potential has an origin in the pion exchange and also heavier mesons 
for a shorter distance. In the mean-field approximation it is assumed that this interaction effects can 
be renormalized in a one-body potential represented by the mean-field variables, so that we can drop 
the pionic degrees of freedom until they are thermally excited at high enough temperature. Then, the 
Walecka model is defined by the Lagrangian density, 

£ = ^{yy^ + /i B 7o - 9*1^ - m N + g a o)$ + ^ (d^ad^a - m 2 (X 2 ) - -u^ u + -mlu^ (45) 

with the following model parameters; the bare nucleon mass mN, the cx-meson mass m a , the cu-meson 
mass the coupling constants in the scalar channel g a , and in the vector channel g^, which are 
determined to fit the nuclear matter properties. We note that ip above represents the Dirac spinor 
for nucleon as if it were a point particle. Since the nucleon is a composite state having a complicated 
internal structure, the nucleon mass at high density cannot avoid ambiguity, which will be the central 



subject of Sec. 3.3.2 



In the mean-field approximation, both a and ojq can take a finite value induced by the non-zero 
chiral condensate (ipip) and the baryon density {ifj^ip). If we are interested in the isospin non-symmetric 
situation such as neutron matter, we have to include the p meson as well. In the above we choose a 
convention in which a is negative (i.e. the same sign as the chiral condensate). In this approximation 
the mean-field Lagrangian density reads, 

£mf = V 5 [iJnd" + (Vb - 9J*>°)7o - (% - g a <r)] V> - \ m l° 2 + ^ m ^( w °) 2 • ( 46 ) 

From this we can immediately identify the renormalized or in-medium nucleon mass and the effective 
baryon chemical potential as 

m* N = m N - g a a , p* B = fJ, B - 9^° • (47) 

Then, the thermodynamic potential divided by the volume V (denoted simply by Q throughout this 
review) for a free quasi-nucleon gas is expressed as 



Q = -2 • 2Tj j0Jfo I 1 + e-^-^} + In [l + e ~^ + ^] j 



mliml - m N ) 2 _ m 2 (^ - p B ) 2 

2& l9l ' 1 



where 2 comes from the isospin degeneracy and another 2 from the spin degeneracy. We note that 
the zero-point oscillation energy is dropped, which corresponds to the approximation in the nuclear 
hydrodynamics [89J . It should be mentioned that one should be careful of the treatment of the zero- 



point oscillation energy if the fermion mass is small, as explained in Sec. 4.2.2 The mean- field variables, 
m N and p B , are to be fixed from the stationary conditions; dfl/dm^ = dQ/dfi B = 0, leading to the gap 
equations, 

m* N = m N - 4^ f ^ P miS [n F (u p - p* B ) + n F (u p + p* B )] , (49) 
9^ f cftp 

fx* B = n B - 4-^ J [n F (up - p* B ) - n F (u p + p* B )] . (50) 
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Figure 6: (Left) Mean-field variables as functions of fi B - There is a first-order phase transition from 
the low-density state in which the nucleon mass is the bare one to the high-density state in which the 
in-medium nucleon mass drops to almost a half. The green curves represent the unstable solutions of 
the gap equation. (Right) Baryon density as a function of the chemical potential. The first-order phase 
transition occurs so that the area in the region I is equal to that in the region II. 



With the solutions of these equations substituted to Q in Eq. (48), we can get the thermodynamic 
potential as a function of T and fiB, from which all thermodynamic quantities are calculable. The 
baryon density is expressed by p = —dQ/dp>B and the internal energy density by s = f2 + Ts + fi-Qp- 

Two coupling constants are chosen so as to satisfy the following saturation properties of nuclear 
matter, that is, 

d(s/p) 



P 



m N - 16.3 MeV 



P=P0 



dp 



. (51) 



P=P0 



We make a remark that the latter is understood as a stability condition for nuclear droplet. In fact, 
this density derivative at T = leads to d(e/p)/dp = Pb/p ~ t/p 2 = p/p 2 = 0. We note that this p is a 
pressure difference between the nuclear droplet and the empty vacuum. When we discuss the liquid-gas 
phase transition of quark matter, we will come back to this point. A set of parameters consistent with 
the nuclear matter properties is 



m N = 939 MeV , m a = 550 MeV , m w = 783 MeV , g s = 10.3 , g u = 12.7 . (52) 



The behavior of the in-medium mass is plotted in the left panel of Fig. [6] together with the 
effective chemical potential p B . It is clear from Fig. [6] that there is a first-order liquid-gas phase transition 
of nuclear matter (or the liquid- vacuum phase transition, strictly speaking). The nuclear matter should 
exhibit a first-order phase transition at fiB = — B/A ~ 923 MeV. In thermodynamics, equivalently, 
we can find this critical value of p, B from the condition that the thermodynamic potential Q becomes 
balanced in two states. In view of the curve of p(pb) m the right panel of Fig. [6] the phase transition 
point is determined so that the area in the region I is equal to that in the region II. This can be 
understood from Af2(/iB c ) = Q(pb c ) — Q(pb c ) = ~ f MBc p(p)dfi = along the path from the lower dot 
to the upper dot in the right panel of Fig. [6j Such a way to find /xb c is reminiscent of the Maxwell 
construction for p and V. 

The left panel of Fig. [7] is the saturation curve of the energy per nucleon as a function of the baryon 



density p. Obviously the saturation point and the binding energy agree with the values in Eq. (51) as 



a result of the fit. In this simplest setup, however, the compressibility cannot be reproduced well, and 



a potential term V(a,u) should be added in Eq. (45) [90], which we will not discuss, for we are mostly 



interested in the qualitative aspects of the phase structure. 
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Figure 7: (Left) The saturation curve of the energy per nucleon as a function of the baryon density 
p. The global minimum is located at the saturation density. The dotted curve represents the unstable 
solution of the gap equation. (Right) First-order phase boundary and the critical point in the Walecka 
model and in the chiral perturbation theory in the three-loop order |91j . 



3.1.2 Liquid-gas phase transition of symmetric nuclear matter 

The existence of the first-order phase transition is an inevitable consequence from the saturation prop- 
erties of symmetric nuclear matter. This world as it is, in other words, embodies the mixed state, and 
in general, a self-bound system of fermionic particles should have a first-order phase transition. 

Differently than the second-order phase transition governed by the global symmetry and its spon- 
taneous breaking, the first-order phase transition in nuclear matter does not change the symmetry 
pattern at all. Because two states before and after the transition have a different density as seen clearly 
in the right panel of Fig. [6j one can regard this type of the transition as associated with the density. 
Then, typically, the dilute phase corresponds to a gaseous state and the dense phase corresponds to a 
liquid state, and nuclear matter hence results in the liquid-gas phase transition. Right at the first-order 
phase transition a mixed state should be formed, during which the gas and the liquid parts can coexist 
in space. The most optimal spatial shape should emerge from the competition between the Coulomb 
energy and the surface energy (see Ref. [92] and references therein), and the resulting peculiar shapes 
are referred to as various "pasta" phases. 

To fit the saturation properties the Walecka model was solved at T = 0, and it is straightforward to 



extend such model calculations to the finite-T situation. We have already written the gap equations (47) 
down with the matter terms. Then, solving them numerically, we can make a prediction for the first- 
order boundary elongated toward T ^ 0. In the Walecka model meson fluctuations are totally neglected, 
and thus, we should not apply this model to the analysis at high T where abundant pions are thermally 
excited. As noticed in the numerical results in Fig.[7j the terminal point of the first-order phase boundary 
(i.e. the critical point) lies around T ~ 20 MeV, which is small enough to justify our neglecting pion 
fluctuations. 

The phase boundary tends to bend toward smaller chemical potential with increasing temperature, 
which can be understood from the Clausius-Clapeyron relation [93]. On the first-order phase transition 
the thermodynamic potential should be balanced between two phases. If we move along the phase 
boundary from T and pe by infinitesimal dT and d/j,B, the thermodynamic potentials, ^dilute in a dilute 
state and fi dense in a dense state should remain equal, which immediately leads to 

^dilute Jrn . ^dilute, C^dense Jrn , ^dense , dT Ap 

-dT- dT+ ^^ d ^ = ^T- dT+ ^^ d ^ d^^-As' (53) 

where Ap = pdense — Pdiiutc is the density difference and As = s den se — ^dilute * s ^ ne en t ro Py difference. 
The slope dT/dji^ is negative in Fig. [7] and this implies that Ap and As should have the same sign. 
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In most cases, indeed, the entropy density becomes greater in the dense state than in the dilute state, 
which means Ap/As > 0. 

From the theoretical point of view the underlying mechanism in quark matter to produce a first- 
order phase transition is exactly the same as that in the Walecka model we elucidated here. In the case 
of quark matter, however, there is no strong support for the existence of the first-order phase transition. 



3.2 Chiral perturbation theory 

One of the most advanced methods to deal with nuclear matter in a relativistic framework is the chiral 



perturbation theory. We have already introduced a Lagrangian (32) that describes low-energy chiral 
dynamics. It was necessary to utilize the linear sigma model because a can be a soft mode at chiral 
restoration. At low temperatures, however, chiral symmetry is significantly broken, and the low-energy 
dynamics is completely ruled by the light NG bosons. Hence, it should be the most appropriate to 
formulate an effective theory in terms of the NG bosons only, which is obtainable by integrating a out 
from the linear sigma model. In such a way, one can come by a non-linear representation of the low- 
energy effective theory. Actually, thanks to the low-energy theorems, one can adopt any representation 
only to get the same physics results (Haag theorem) [M], which provides us with a systematic and 
model-independent tool to disclose the properties of hadronic matter. 

3.2.1 Effective Lagrangian and the power counting 

In the vacuum we do not have to consider the scalar particles because they are much heavier than the 



NG bosons. Then, as seen in Eq. (33), the theory is defined as a function of r • 7r where r is the Pauli 
matrices in isospin space in the iVf = 2 case. It is the most useful choice to employ a representation in 
terms of unitary matrix, 

U = e *r-*/U } jj —tjJuR, (54) 

having the same transformation property as M with the a meson dropped off. The chiral-symmetric 
Lagrangian density is to be expanded in terms of U and derivatives; the lowest-order is 

4 2 > = ^ftv(D^UD^) + &tr( X U^ + U X ] ) , (55) 

with two derivatives and and one symmetry breaking term. The parameter in the latter term is x = 
2Bq(s + ip) where s is the quark mass matrix. Using Bq one can express the chiral condensate as 
(qq) = —f%Bo and the pion mass as m 2 = (m u + ui^Bq which is consistent with the Gell-Mann-Oakes- 
Renner relation (24). The covariant derivative is defined as D^U = d^U + il^JJ — iUr^ that transforms 



as D^IJ — > lJ{x)(DJJ)R(x) with properly defined and r M for local chiral rotations L(x) and R(x) 
as in the gauge theory. 



The chiral effective Lagrangian (55) can be regarded as the lowest-order 0(p 2 ) contribution in the 
expansion with respect to pion momentum scale p as compared to the typical chiral scale A x = 4:7rf n ~ 
1 GeV. Then, the pion mass m n is 0(p) and then x ^ m 1; is @{p 2 )- We can perform the theoretical 



calculation for physical observables up to 0(p ) using this Lagrangian (55) at the tree-level. To go to 



the next order of 0(p 4 ), one should include one-loop diagrams from cW together with the tree-level 
contributions from at 0(p 4 ). In the SU(iVf = 3) case there are 10 independent terms (including x, 
^At> r fi) an d 10 low-energy constants are required at 0(p 4 ). In the SU(2) case there are less parameters. 
These constants are fixed to fit with the experimental data of the meson masses, it-it scattering, rare 
pion decays, decay, fx/ fn, etc, and then one can carry theoretical predictions out. This systematic 
program of the low-energy expansion of the effective theory of QCD in terms of the NG bosons is called 
the chiral perturbation theory (ChPT). 
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Figure 8: (Left) Chiral condensate at finite T and p in the limit of the free pion gas and the linear 
density approximation given in Eq. (58). (Right) Chiral condensate calculated in the NLO three- loop 
order. Figure adapted from Ref. [TOO] and the linear density approximation results (dashed lines) are 
overlaid for reference. 



Let us consider how to apply the ChPT to figure out the behavior of the chiral condensate at finite 
temperature and density. Because the operator ipip and the quark mass m q are conjugate to each other, 
the chiral condensate (ipip) is derived from the derivative of the thermodynamic potential Q with respect 



to m q , that is, (ipvp) = —dVt/dm^. Thanks to the GOR relation (24), it is a straightforward procedure 
to convert the m q -derivative to that in terms of the pion mass squared m%.. Then, in a medium at finite 
T and density p, the deviation of the in-medium chiral condensate from the vacuum value (V'V')o can 
be expressed as 

(iH>)o ^ \dml ml^m N ) , 

where the second contribution has the nuclear sigma term that is defined and estimated as [95j EE] , 

ovn = m q 5^ ~ 45 MeV . (57) 
om q 

The lowest-order correction originates from contributions induced by T from the free pion gas and 
p from the free nucleon gas. If the pion mass can be negligible (i.e. ~ 0) and p is sufficiently small 
enough to approximate the scalar density p s = —dVt/dm^ with the density p, the analytical integration 
is possible, leading to [97, 98] , 



(TOO 8/2 m 2/2 

at the lowest order in the linear density approximation. The finite-m^ correction is not really negligible, 
in fact, and the expansion breaks down around T ~ 150 MeV due to the neglected excitations of other 
mesonic states. Because the a meson decouples from the dynamics, one has no chance to describe the 
chiral phase transition using the non-linear representation. Nevertheless, it would be suggestive to see 
explicitly where the chiral condensate vanishes with increasing T and p according to the lowest-order 



expression (58). The results are depicted in the left panel of Fig. 8 It seems that the chiral condensate 



disappears to draw a sort of the phase boundary on the T-p plane, though the approximation in 



Eq. (58) breaks down there. Even though we cannot reach the chiral phase transition in the ChPT, 



such in-medium reduction of the chiral condensate or /„- can be an experimental measure for the partial 



restoration of chiral symmetry [99] (see also Eq. (66)). 
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3.2.2 Application to nuclear matter 

It is necessary to take account of the heavy baryonic degrees of freedom into the dynamics at high 
density. The fermion part is expressed by the Lagrangian, C = ^[i^ — g(a — it ■ tt)]^, which is 
converted into the non-linear representation after dropping the scalar field. The usual procedure to 
write the baryon theory down is to perform the field redefinition as 

£ = e W2A j iV L = eV L , Wr = £Vr, U = tt. (59) 

Then the fermionic Lagrangian is translated into the form similar to the gauge theory, C = N(ip — 
— m)N, with = <9 M + iV^ and the vector and the axial-vector fields, 

V, = - l -{t}d»t + , A, = -\{gd»t - . (60) 

Including the renormalization constants that change the axial-vector coupling from unity to — 1.26, 
the nucleon effective Lagrangian is expanded in terms of the pion fields up to the quadratic order as 

£ N = N{i$ - m N )iV + ^NYl^N ■ - j^N-fr ■ tt x d^N + ^§tt 2 NN (61) 

Jir 2 4/^ Zj n 

with the pseudo-vector 7rA^A^-vertex and the Tomozawa- Weinberg 7T7riViV-contact vertex. With these 
interaction vertices the thermodynamics and the chiral condensate are derived; the next-to-leading- 
order (NLO) results are calculated and shown in the right panel of Fig. [8] adapted from Ref. jlUUj . The 
linear density approximation comes to work fine for T > 50 MeV, while the T-dependence at NLO 



shows a significant deviation from Eq. (58). 



From the point of view of the systematic expansion with respect to the energy scale, it is not straight- 
forward to formulate perturbative expansions involving baryons with the consistent power counting. A 
new scale, i.e. the nucleon mass m^, is not small as compared to the chiral symmetry breaking scale 
A x and appears in higher loops that contribute to lower order. Thus, the apparent one-to-one corre- 
spondence between the loop and the derivative expansions is, at first glance, lost. Yet, within the fully 
relativistic framework one can retrieve the consistent power counting for multi-loop diagrams using ap- 
propriate renormalization conditions |101[ 1102] . Another resolution to avoid this drawback is to adopt 
a heavy baryon reduction yielding the non-relativistic limit of the original theory, that is, the heavy 
baryon chiral perturbation theory (HBChPT) where the power counting is restored |103j . We write the 
nucleon momentum as 

p» = m NV » + p , (62) 

where is the four- velocity with v 2 = 1 and is the residual momentum of order Aqcd, so that one 
can perform the derivative expansion systematically in energy range below A x . The Lagrangian density 
is written in terms of the nucleon field B carrying which is related to the original field via 

B{x) = expOim^v ■ x) N(x) . (63) 

With this the nucleon part of the Lagrangian density becomes 

£ N = N (i'fd^ -m N )N = B iv^d^B . (64) 

Therefore no hard scale appears in the theory. Higher-order terms are arranged in powers of k^/A x 
and P/fflm. Loop calculations in the HBChPT framework are simpler as compared to the relativistic 
treatment and the HBChPT has been applied to a broad range of the nuclear physics problems. For 
reviews, see Refs. [TM IT051 IT06] . 
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In the nuclear medium the methodology of chiral effective field theory can also be formulated. The 
in-medium ChPT Lagrangian is derived from the generating functional after integrating the nucleon 
fields out but keeping the particle-hole part of the fermion determinant. The in-medium pion decay 
constants, 



(N\J A (x)\7r(q))=iq°fle^ , (N\J A (x)\n(q)) 
are found at one-loop order in symmetric nuclear matter as [107J 



flip) = u 



1 - —(0.26 ±0.04) 
Po 



fM = U 



p 



1 - —(1.23 ±0.07) 
Po 



(65) 



(66) 



A model-independent analysis valid at low density leads to a scaling relation between /* and the in- 
medium quark condensate |108j . with which one sees that the condensate decreases in nuclear matter, 
indicating partial restoration of chiral symmetry (see Ref. |109] for a review). 



3.3 Mesons, baryons, and exotica 

From the ChPT one can calculate how the chiral condensate decreases. If the chiral condensate signif- 
icantly drops, then, it should have sizable influences on the hadron spectra such as the nucleon mass. 
Indeed we have already elucidated that the in-medium effective mass in the Walecka model is different 
largely from the bare value. The question is, then, how chiral symmetry becomes manifest for the 
hadronic degrees of freedom if (partial) chiral restoration occurs in dense matter. 



3.3.1 Mended symmetry 

In general a large variety of representations are possible for composite states made from quarks that 
belong to the fundamental representation. When chiral symmetry G were not spontaneously broken 
down into H, the physical states would be classified by irreducible representations of G. With broken 
chiral symmetry, the physical states belong no longer to those irreducible representations, but instead 
are superposition of all the possible representations. The physical state |a) is expressed in a series of a 
complete set of states | J) characterized by good quantum numbers J; 

i«) = E c ^ J )- ( 6? ) 
j 

The coefficients cj contain information about the broken chiral symmetry. The state \a) can further 



be constrained by applying the large- N c approximation (see Sec. 3.5 also) as well as current algebraic 
relations |110[ 1111] . The states for the low- lying excitations with zero helicity are classified by a limited 
number of the irreducible representations and form a closed algebra which is referred to as mended 
symmetries. 

As shown in Ref. |110j . assuming that the Adler-Weisberger sum rules for the scattering process, 
7r a + a(X) — > n b + for helicities A and A', can be saturated with narrow one-particle states, the set 
of the sum rules is put into the Lie-algebraic form, 

[X a ,X b ] = if abc T c , (68) 

where f abc is the structure constant, T c is the generator of the SU(iVf)v group that satisfies 

[T a ,T b ] = ifabcTc , (69) 

and X a is the axial- vector coupling matrix defined from the matrix elements at zero invariant momentum 
transfer of the axial-vector current between states with collinear momenta via 

(q\'(3\(J° Aa + J 3 Aa )\p\a) = (2n)*2p+5 3 (p- q)5 yx [X a (X)]^ a . (70) 
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Thus, from this definition, some non-trivial calculations lead us to the following commutation relation, 



[T a ,X b ] =if abc X c . (71) 



Putting together, Eqs. (68), (69), and (71) close the algebra and one sees that the X a together with the 
unbroken SU(iVf)v generators T c satisfies the commutation relations of the Lie algebra of the broken 
symmetry group G. 

In the broken phase the X a is not symmetry generators. In fact, it does not commute with mass- 
squared matrix m 2 . We note that m 2 here refers to the general hadron masses and is not necessarily 



vanishing in the chiral symmetric phase. The non-commutativity which follows from Eq. (68) along 
with the Jacobi identity defines m|, that is, the chiral-symmetry breaking part of m 2 as 

[X b , [m 2 ,X a ]] = -m\b ah . (72) 

In the iVf = 2 case, for example, 777.4 transforms as the fourth component of a chiral four-vector (like the 
a meson). This tells us that the mass matrix is a sum of a chiral-singlet m\ that satisfies [X^m 2 ,] = 
and a non-singlet m 2 , i.e. 

m 2 = ml + ml . (73) 

When m\ vanishes, i.e. chiral symmetry is restored, X a becomes true symmetry generators, so that X a 
commutes with m 2 and hadrons form degenerate multiplets filling out a complete set of representations 
ofG. 

In the chiral broken phase the algebraic representations of X a do not always coincide with the mass 
eigenstates as accounted for by [X a , m 2 } ^ 0. In general, there, physical particle states are expressed as a 
sum of all possible elements of various representations of the Lie algebra composed of X a and T c . The X a 
is therefore entirely determined by the mixing angles which define the coefficients of the representations 
in the sum. Let the physical states of the low-lying A = mesons, the scalar cr, pseudo-scalar ir, vector 
p and axial- vector a\ mesons, be in the following admixture of the chiral representations |110[ 1112] . 



|tt) = | (iVf, Nf) © (iV f *, N{)) sin 6 1 + | (iV 2 -l,l) © (l, iV f 2 -l)) cos 9 
\ ai (X = 0)) = | (Nf, N£) © (jV?,jV f )>cos0- | (iV f 2 -l,l) © (l,iV f 2 -l))sinfl 

\a) = \{N { ,N?)®(N?,N { )) , 
|p(A = 0)) = |(iV 2 -l,l)ffi(l,iV 2 -l)). 



(74) 



This quartet structure based on the chiral algebraic sum rules is known as a notion of mended symme- 
try [111] . This leads to the decay widths as functions of the mixing angle 9 and masses as well as f n . 
The experimental rate for p — > tt + 7r and the pion decay constant f„ give approximately 9 ~ 45° and 
one then finds 

m p ps m a , m 2 ai ps 2m 2 , (75) 

indicating that the above algebraic consideration indeed yields the reasonable order of magnitude for 
the masses. 

Apparently, those helicity-zero mesons as well as the pion all become massless at a second-order 
phase transition where m a becomes vanishing |lllj . When an external parameter, such as T and p, 
drives a system from broken to unbroken phase, the mixing angle is in general supposed to have some 
intrinsic dependence on the external parameter. Such medium evolution of the mixing angle needs to 
be pinned down so that some signals of the chiral phase transition in hadronic quantities would be 
revealed. 
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3.3.2 The problem of mass 



In introducing baryon fields, there are two alternative ways of assigning chirality to the nucleons. Let 
us consider a pair of chiral partners, ip\ and ip2, that realizes chiral SU(2)l x SU(2)r. One is, 

(Standard Assignment) ip 1L -)> Lip 1L , ip 1R ->■ Rip 1R , 

(7b) 

1p2L -> £^>2L , ^2R -> -R^2R , 

which is anchored on the standard chiral symmetry structure where the entire constituent quark or 
nucleon mass (in the chiral limit) is generated via spontaneous chiral symmetry breaking. The other 
is [TT3] . 

(Mirror Assignment) ipxL -> Lip 1L , ip 1R -)■ iityiR , 

V>2L -> /'' '2I. , '-'21! -> /-' '2H , 

which allows for a chiral invariant mass term, 

£ m = WO (V'2LV'lR - ^2R^1L ~ ^1L^2R + ^IrV^l) , (78) 

with which m remains non-zero at chiral restoration. One should notice that the left- and right- 
handedness of a fermion is the representation of the Lorentz group, whereas chiral symmetry is entirely 
associated with the internal symmetry. Therefore those must be in general dealt with independently, 
which enables us to adopt the above non-standard chiral rotation. In other words, the mirror assignment 
is possible since the nucleon is not an elementary particle but a composite state made from quarks. 

A part of the nucleon mass, mo, must arise from a mechanism that is not associated with spontaneous 
chiral symmetry breaking. Various observables both in the vacuum such as baryon decays |114j and 
pion-nucleon scattering |115[ 1116] and in a medium such as nuclear matter properties |117l 11181 1119] I120j 
and the thermal gluon condensate [121J have been explored based on chiral Lagrangian approaches with 
mirror baryons. The analysis tells us that an mo of a few hundred MeV cannot be ruled out. In the 
non-linear realization of chiral symmetry only the properties coming from the vectorial transformation 
are remained |122j . and therefore the standard and mirror assignments cannot be distinguished. Even 
if indistinguishable in matter-free space, they are expected to start showing their differences due to the 
remaining mo as one approaches the chiral restoration point. 

What is the origin of such a mass mo? Recall that the emergence of Aqcd is signaling the breaking 
of conformal invariance, i.e. the trace anomaly, of QCD in the chiral limit where the theory has no 
dimensionful parameters. Spontaneous chiral symmetry breaking, which gives rise to a nucleon mass, 
and the trace anomaly are intimately linked to each other |123] and dynamical scales in hadronic 
systems are then considered to originate from them. Hence, the origin of mo can be attributed to the 



non- vanishing gluon condensate in chiral symmetric phase |121] (see also discussions in Sec. 3.3.3). 



This is reminiscent of Weinberg's chiral-singlet mass in Eq. (73). The chiral representation mix 



ing (74) would be expected to be resolved when chiral symmetry is restored in such a way that the 7r 
and a mesons form a degenerate multiplet as chiral partners, as presumed from the argument of the 
0(4)-universality class in iVf = 2 QCD. This is in fact realized when the mixing goes like cos (9 — > 
and with which the it and a belong to the same representation \(Nf, N£) © (-/V f *, Nf)) and being massless 
protected by the NG theorem, whereas this is not the case for the p and a\ mesons being degenerate 
in the same representation at chiral restoration and in general having a non-vanishing common mass. 
It is also possible that those lowest-lying mesons form a quartet being all the mesons massless when 
the mixing remains unresolved as demonstrated using the mended symmetry notion in Ref. |lllj where 
fri'P = fa a is predicted. Or even with a resolving mixing, if there are some mechanisms that the p and a\ 
mesons become massless, the massless quartet can be well accommodated to the requirement of chiral 
restoration, such as in the context of the Brown- Rho (BR) scaling |124j . an analysis of RG flows within 
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a class of chiral effective field theory [1251 1126] and of the non-linear sigma model implementing the 
trace anomaly |121j . 

Yet, massive vector bosons with no mixing at chiral restoration are not prohibited by symmetries. 
The mixing angle can be in general a function of temperature and density and therefore would develop 
in a hot/dense medium. Even if the mixing is preferred being m p = m a in matter-free space, it might 
evolve so that cos (9 — y 0, not optionized in Ref. [Ill] , namely = m a = and m p = m ai ^ can be 
allowed. 



3.3.3 Thermodynamics of dilatons 

The trace anomaly is implemented in a chiral Lagrangian by introducing a dilaton (or scalar glueball) 
field x representing the gluon condensate (G^G^) |127[ 1128] . Requiring that the dilatation current 
associated with scale transformation is saturated by a condensed dilation field leads to the potential, 



* 4 Vxo 



i 



In ^ -1 



(79) 



The two parameters, i.e. the bag constant B and a scale parameter xo, are fixed from the vacuum 
energy density £ = \B and the glueball mass M x = 2VB/xo- Effective Lagrangians are written so that 
they possess chiral and scale invariance, mimicking dynamical breaking of those symmetries of QCD in 
low energies. 

In relativistic nuclear physics, a scalar boson plays an essential role as a mediator of the medium- 
range NN interaction as we discussed using the Walecka model in Sec. |3.1 The important lesson 



from the model studies is that this scalar degree of freedom cannot be identified with the scalar meson 
that appears as the chiral partner of the NG boson in a linear sigma model since otherwise no stable 
nuclear ground state is achieved and one cannot reproduce the basic saturation properties |129[ |130j. 
This problem would be resolved in a non-linear representation as in the ChPT, or different chiral 
approaches have been introduced with modified interaction potentials among the chiral fields and the 
dilaton [130. 1131] 1132] . A non-linear chiral Lagrangian with a soft dilaton [133] is derived when a heavy 
component of the dilaton, responsible for the pure glueball dynamics, is integrated out [134J. 

Besides those approaches implementing the dilaton, a parity doublet model based on the mirror 
assignment can also be accommodated to the nuclear ground state [118] . The presence of a chiral 
invariant mass m of the nucleons is essential there and a relatively large m ~ 800 MeV is favored in 
order to have the known nuclear matter properties. (This number becomes reduced substantially when 
a tetra-quark state is introduced and mo in nuclear matter results in ~ 500 MeV [120] .) 

It appears that mo plays a similar role to the dilaton in nuclear matter and thus a non-vanishing 
condensate of the dilaton field is a conceivable origin of itlq [1211 1116] . The scalar meson as the chiral 
partner of the pion generates only the mass difference of the nucleon parity doublers. 



3.3.4 Four-quark condensate in a medium 

Four-quark states were suggested as the lightest scalar mesons |135[ 1136] and this explains naturally the 
degenerate a (980) and / (980) states. How does the four-quark condensate evolve with temperature 
or density? The NG theorem by itself does not dictate a particular operator to fulfill a non-trivial 
commutation relation with the broken charge Qa- Any operators O yielding non-commutativity to Qa 



guarantee spontaneous chiral symmetry breaking as discussed in Sec. |2.2.1| Besides O ~ ipt a ip, some 
operators with higher dimensions such as a four-quark are possible. 

In matter-free space the four-quark condensates are often evaluated assuming factorization, 

$$^) - (V^> 2 • (80) 
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This indicates that the four-quark condensate seems to vanish at the chiral phase transition as (ipip). 



However, the factorized form (80) is justified in the large- N c limit only and not quite accurate in the 
real QCD case with N c = 3. Its validity in a hot/dense medium is even less justifiable. 

Therefore, the issue is that when (ipip) vanishes, some four-quark operators (apart from the chiral 
singlet part) can condense so that it breaks chiral symmetry dynamically. Such a possibility has been 
considered at zero temperature and zero density in Refs. |137| 1138] , where a non-standard pattern of 
dynamical symmetry breaking was suggested. This pattern keeps the center of chiral group unbroken, 
i.e. 

SU(iV f ) L x SU(iV f ) R -> SU(JV f ) v x Z Nt . (81) 

The discrete symmetry protects a theory from condensate of quark bilinears. Spontaneous symmetry 
breaking is driven by quartic condensates which are invariant under both SU(A r f)v and Z^ { transfor- 
mation. However, this pattern is strictly ruled out in QCD both at T = and T ^ but /i q = since 
the axial-vector correlator becomes greater than the pseudo-scalar correlator then, which violates the 
QCD inequality |139j . This argument does not exclude the above unorthodox pattern in the presence 
of a dense medium because the QCD inequalities do not hold then. Several studies have been made, 
such as a similar dynamical breaking in an 0(2) scalar model [140J, constructing a Ginzburg-Landau 



free energy consistent with Eq. (81), and the exploration of its thermodynamic consequences |141j . 



3.4 Hidden local symmetries 

The low-energy pion dynamics is well described by the non-linear sigma model on the manifold G/H 
which is an effective Lagrangian in terms only of NG bosons transforming non-linearly under chiral 
rotation. The approach using the hidden local symmetry (HLS) |142j 1143] is based on the fact that 
the non-linear sigma model is gauge equivalent to a model possessing G g i bai x H\ oca \ symmetry (see 
Ref. |144j for a review and references therein). As demonstrated explicitly in Refs. [1451 I146j . the 
dynamical generation of gauge bosons is a common phenomenon which appears in a variety of systems 
not restricted to particular models. Furthermore, one-loop quantization of the HLS theory based on 
the derivative expansion as made in the ordinary ChPT is established |144j . This enables us to evaluate 
loop effects systematically and thus the low-energy quantities with good precision come out combined 
with the RG analysis. Below, we briefly introduce the HLS approach and discuss several issues in strong 
coupling gauge theories. 



3.4.1 G g iobai x I/iocai model 

The Lagrangian is based on a Global x # local symmetry, where G g i bai = [SU(iV f ) L x SU(A^f) R ] g i oba i 
is chiral symmetry and i?bcai = [SU(iVf)v]iocai is the HLS. The entire symmetry G g i bai x H\ OCSu \ is 
spontaneously broken to its diagonal subgroup SU(2) V . The basic quantities are the HLS gauge boson, 
Vp, and two matrix valued variables £l, £r, which are combined in an N{ x Nf special-unitary matrix 
U = £l£r- There exists an ambiguity in the decomposition and this redundancy is identified with a 
local gauge transformation Hi oca \. These variables are conventionally parameterized as 

£ L R ( X ) = e i"(*)/U f« ) (82) 

where ir = n a T a denotes the pseudo-scalar NG bosons associated with the spontaneous breaking of 
chiral symmetry, and o = a a T a the NG bosons associated with the spontaneous breaking of Hi oca \. The 
a is the longitudinal part of the vector meson and absorbed into the HLS gauge boson through the 
Higgs mechanism, so that the gauge boson acquires its mass. f n and f a are the decay constants of the 
associated particles. 



32 



The main building blocks are defined from the Maurer-Cartan 1-forms as 

l - (D"£ R ■ ^ =F D^ L ■ ej) , (83) 



a. 



2 

where the covariant derivatives of £l,r are given by 

with £ M and IZ^ being the external gauge fields (weak bosons and photons) introduced by gauging 
Ggiobai- The most general Lagrangian with lowest derivatives is given by |142j 

c = + /M<vf) ~ ^HW U ) > ( 85 ) 

where g is the HLS gauge coupling and the field strength is V^ u = d^V u — d u V^ — i\V^ V„]. When the 
vector-meson kinetic term is ignored in low energy, the second term of the Lagrangian vanishes with the 
equation of motion of and the well-known non-linear sigma model on the coset space G/H is derived. 
The leading-order Lagrangian contains two arbitrary parameters, a = f%/f% and g. For a particular 
choice a = 2, the known phenomenological facts are successfully reproduced; universality of the vector 
meson coupling, the Kawarabayashi-Suzuki-Riazuddin-Fayyazudin (KSRF) relation, and vector meson 
dominance of the pion electromagnetic form factor |142j . 

There are other models of vector mesons introduced as the matter field, the anti-symmetric tensor 
field or in the massive Yang-Mills method. It is shown that they are equivalent to the gauge-fixed form 
of the HLS Lagrangian. For details, see e.g. Refs. |143| 1147^ \1AS\ I144j and references therein. 

3.4.2 Modeling strongly-coupled gauge theories 

In principle one can enlarge the gauge group with more redundant variables. With such a generalized 
HLS (GHLS) |143| 1149] . the lowest-lying axial- vector mesons are embedded in a G g i bai x Gi oca ,\ model. 
Other higher-lying vector modes can also be incorporated systematically by introducing enlarged gauge 
symmetries. This idea is employed in so-called dimensional deconstruction which resembles the nearest 
neighbor interactions in condensed matter physics, constructed on a lattice |150[ 1151] 1152] . The theory 
possesses gauge group which is a product of many numbers of the same gauge group. The interaction 
to matter fields is characterized by "theory space locality" |153j where the mixing of left and right 
chirality is generated only through gauge boson exchanges. Taking the continuum limit generates an 
extra dimension. Inspired by those studies, a dimensionally-deconstructed QCD was proposed [154 



The resultant five-dimensional Yang-Mills action is shown to reproduce some of known low-energy 
hadron properties. 

The identical action implemented with the Chern-Simons action is obtained from the string theory 



based on the gauge/gravity duality (see Sec. 3.5.3 for the Sakai-Sugimoto model) [2"3" t I155j . Performing 
the Kaluza-Klein mode expansion yields a four- dimensional action which contains the NG bosons and 
an infinite tower of vector and axial-vector mesons. When all the vector modes except the lowest-lying 
vector meson are integrated out, the 0(p A ) terms as well as the 0(p 2 ) terms of the HLS Lagrangian are 
generated |156j . The higher vector mesons are now accommodated in the 0(p 4 ) part. 

Another type of duality of interest is the electric/magnetic duality in M = 1 supersymmetric Yang- 
Mills theory (i.e. Seiberg duality) |157j . For |iV c < iVf < 3iV c (conformal window), there exists a 
magnetic theory possessing the SU(iVf— N c ) gauge symmetry which is dual to the original SU(iV c ) theory 
regarded as the electric theory. The two theories flow to the same IR fixed point and thus exhibit the 
identical low-energy dynamics. When Nf decreases to Nf = N c , the magnetic gauge symmetry induced 
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at the composite level is spontaneously broken, which corresponds to the system with confinement and 
spontaneously broken chiral symmetry. It should be emphasized that the induced gauge symmetry has 
nothing to do with the original SU(iV c ) gauge symmetry, but it is rather analogous to the HLS |158j . 

Is the Seiberg duality realized in non-supersymmetric QCD under consideration? The answer is 
yes |159[ 1160] : The HLS theory is shown to exhibit chiral symmetry restoration by its own dynamics for 
large Nf ~ 5 (iV c /3), as in the CP N ~ X non-linear sigma model |143j . Therefore the SU(iVf) HLS theory 
is considered to be the Higgsed magnetic theory dual to the electric theory, i.e. SU(iV c ) QCD. Some 
consequences of this "QCD/HLS duality" in hadron physics are discussed e.g. in Refs. [1601 1161] . 

Here we stress that the HLS theory appears as a potentially practical approach for various gauge 
theories in the strong coupling region. Obviously low-energy hadron physics, not only in matter-free 
space but also in a hot and/or dense medium, is a case of great interest. The HLS approach has been 
applied to QCD matter in a wide range of density and temperature. For reviews, see Refs. |162[ I163|, 
11641 11651 11661 1167] and references therein. 



3.4.3 Chiral perturbation theory with HLS 

Beyond the tree-level HLS results, loop effects are systematically evaluated at a given order of the 
derivative expansion as in the ordinary ChPT [168] 1169] . The essential concept is the HLS gauge 
bosons treated as light. The chiral counting of the HLS gauge coupling is assigned of order p, which 
is equivalent to the order assignment by m v ~ 0(p) \170\ 1171] . The pionic coupling is derivative-type 
typically in non-linear realization. With the above counting, the vector mesons and their interaction 
can naturally be put in a perturbative treatment on the same footing with the NG bosons, yielding 
consistent counting rules. In real QCD the vector meson mass is not very small compared to the chiral 
symmetry breaking scale A x ~ 1 GeV. However, there exists a limit of large- N c in which the above 
expansion is certainly justified since 

m?/ 1 , , 

— . (86) 

A 2 X N c 1 ] 

Thus the entire procedure for the loop calculations should be regarded as an extrapolation from large- N c 
to N c — 3 QCD. Indeed, combined with RG equations (RGEs), this works quite well and one sees the 
low-energy constants and the relevant quantities in the vector meson sector in good agreement with 
measurements [144J. 

Quadratic divergences that arise from the pion loops drive a system to chiral symmetry restoration 
(Wigner-Weyl realization). In the HLS theory this happens with the following conditions dictated by 
matching of the current correlation functions to those in the operator product expansion (OPE) |172j : 

a 1 , m y -> . (87) 

The vector meson thus becomes massless toward chiral symmetry restoration where f n = 0, and falls 
into the same chiral multiplet together with the pion. Such somewhat non-standard realization called 
the vector manifestation (VM) is similar to the vector realization. The essential difference is that in 
the vector realization the NG boson pole residue f n remains non-vanishing whereas chiral symmetry 
is "unbroken", i.e. f n ^ 0, (ipip) = 0, therefore the system is not in the Wigner-Weyl phase (see 
Refs. [1 701 1171] for more details). Here one should be cautious of the fact that although the two 
realizations are formulated in the HLS framework, the theory contains a limited number of low-lying 
mesons, pseudo-scalar, and vector mesons, and in particular lack of a scalar particle near the restoration 
point is unsatisfactory from the point of view of building effective theories that describe the chiral phase 



transition. In fact, the two conditions in Eq. (87) correspond to the fixed points of the RGEs separately. 
Thus higher loops or other resonances would modify the RGEs and consequently parametric conditions 
of symmetry restoration. Nevertheless, the power of this effective field theory is that one can examine 
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several scenarios allowed as phenomenological options, including the my 7^ case, within the single 
framework |173j . How the scalar mode emerges from non-linear sigma models is an issue |1741 1121] and 
a naive introduction of an elementary scalar field is known to cause the naturalness problem at loop 
level as in the Higgs sector of the Standard Model. 

An interesting extension of the phase structure study was carried out based on the GHLS at one 
loop \1'25\ I126| . The chiral representation mixing between the pseudo-scalar and axial- vector meson 



states (74) is present in the theory and its RG evolution indicates how the mixing is resolved toward 
chiral symmetry restoration driven by external parameters such as N{, temperature, and density. Three 
associated fixed points are found and their characteristic features are summarized below with the vector 
and axial- vector meson masses my,A' 

Ginzburg- Landau type cos9 — > thus (my/m^) 2 — y 1 VD fulfilled 

Vector Manifestation type sin^— > thus (my/m^) 2 — >■ VD violated 

Hybrid type sin 9 — > ^/l/3 thus (m v /m A ) 2 ^ 1/3 VD violated 

where VD represents the vector dominance. Which realization is eventually chosen is a question of the 
actual QCD dynamics. Violation of the VD in the latter two scenarios is rather strong, in particular 
near the chiral restoration point; violated by 50% in the Vector Manifestation type and 33% in the 
Hybrid type. Its impact on QCD phenomenology is obvious. One example is dilepton yields from 
in-medium p mesons which are suppressed when the VD is violated, whereas the VD is often assumed 



to be valid at any temperature and density in many calculations. This will be discussed more in Sec. 5^2 
Another important finding in the GHLS approach is the spectral function sum rules that are locked 
to theory space locality. The extended condition for locality in the GHLS is shown to be stable against 
the one-loop corrections |125] and this immediately leads to the celebrated Weinberg sum rules in the 
pole-saturated forms |175j . 

fl + fl = fv, fK = fvK , (88) 

with the masses and decay constants of the corresponding particles. In other words, the sum rules 
invariant under the RG evolution are guaranteed by theory space locality. 



3.5 Approaches in the large- N c limit 

We here give detailed explanations on a non-perturbative approach based on the limit of the infinite 
number of colors (see Ref. [T76J for a recent review). The large- N c limit was first considered on the 
diagrammatic level |177[ 1178] , and such diagrammatic considerations are still useful for some kind 
of duality in supersymmetric and non-supersymmetric Yang-Mills theories (i.e. planar equivalence; 
see Ref. |179] for example). For non-perturbative applications we will go specially into the Skyrme 
model |180] and the Sakai-Sugimoto model [23] below. 

3.5.1 Counting rule 

If the number of colors N c is chosen to be infinitely large, the theoretical treatment becomes simplified, 
and it is even possible to extract some model-independent results. The naive limit of iV c — > 00 makes, 
however, all physical observables blow up to infinity. First, let us consider the structure of the leading 
diagrams, and the physically sensible way to take the the large- N c limit. 

It is a common procedure to introduce the double-line notation in which each line refers to the index 
in the fundamental representation. Thus, a quark ipi and an anti-quark ipj are represented by a single- 
line. Because the adjoint representation is composed from the non-singlet part with the fundamental 
and anti-fundamental indices, the gauge field = A a T°- is represented by a double-line, as illustrated 
in the left of Fig. |j 
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Figure 9: (Left) Double-line notation: a quark and anti-quark are represented by a single-line, while 
the gauge fields involve two lines as a combination of the fundamental and anti-fundamental indices. 
(Right) Quark loops are suppressed by 1/N C as compared to gluon loops in the large- N c limit. 




Figure 10: N c counting for the meson two-point function, (a) Tree-level diagram, (b) First loop 
correction g 2 N 2 ~ N c . (c) Next planar loop correction g A Nj} ~ N c . (d) Non-planar loop correction 
g*N c ~N-K 



In the large- N c limit, using the double- line notation, it is easy to make sure that quark loops are 
always suppressed by 1/N C as compared to gluon loops and thus diagrams with only gluon loops are 
dominant (see the right of Fig. [9J. This is a crucial feature in thinking of the thermodynamics of QCD 
in the large- iV c approximation. In other words, whenever we talk about the quark and baryon properties 
in the large- N c limit, we treat such baryonic degrees of freedom as a probe into the gluonic system that 
is not perturbed by the insertion of baryons (i.e. probe approximation). 

The next important observation in the large- N c limit is that the so-called non-planar diagrams are 
suppressed by at least 1/N 2 as compared to the planar diagrams. To see this, let us consider the two- 



point function of the meson composite field M ~ ipip as drawn in Fig. 10 . It is obvious that the tree-level 



diagram connecting two meson fields consists of simple one-loop with two lines, making a contribution 



of 0(N C ) (see Fig. 10 (a)). The next contribution at the two-loop level is shown in Fig. 10 (b), that 



has an extra factor by g 2 N c than the tree-level contribution. Therefore, in this way, every time the 
loop order is incremented, a factor g 2 N c at least is multiplied. Therefore, to construct the meson two- 
point function in a non-perturbative and non-diverging way, the strong coupling constant must scale 
as g 2 ~ 1/N C (or so-called 't Hooft coupling A = g 2 N c is kept fixed), and then the contribution from 
Fig. 10 (b) can be of the same order as (a). The next loop order consists of the planar diagram (c) and 
the non-planar diagram (d) with an overhead crossing. Now that g 2 ~ l/N c , (c) is at the same order 
as (b), while (d) is suppressed by l/iV c 2 . 

The meson two-point function (M(k)M(—k)) is of order of N c as is clear from Fig. 10, that means 
that the meson decay constant scales like f n ~ \/Nc- From this, it is possible to identify the scaling of 
the inter-meson coupling. That is, let us consider the n-point meson correlation function, that involves 
(fn) n ~ N^ 2 and the n-point interaction vertex g n . The whole planar diagrams should scale as N c , 
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and therefore, g n ~ iV c . Hence, the inter-meson vertices are all vanishing for n > 3 in the limit of 
N c — > oo. The baryon, on the other hand, is composed from N c quarks and the baryon mass is thus 
of 0(N C ) (there are N C (N C — 1) pairs of inter-quark interactions, but the interaction strength scales as 
1/N C and then the total scaling is 0(N C )) |178j . Then, in QCD with N c — > oo, baryons are infinitely 
heavy and the theory is reduced to a system with infinite towers of mesons and glueballs. 



3.5.2 Skyrme model 

One could construct the baryon as a bound state of N c quarks. Then, the baryon mass is naturally 
0(N C ) ~ 0(l/g 2 ), while the interaction of low-lying particles is 0(1/N C ) ~ 0(g 2 ). This inverse- 
correlation is reminiscent of the properties of the soliton. Indeed, it is possible to construct the baryon 
as a soliton in terms of the low-lying fields. 



Single skyrmion Here, we shall limit our discussions in the simplest example of the SU(2) case. The 
low-energy effective Lagrangian in the non-linear representation again is 



C= S -f\ X {d li UPrf) + J-j 



\x(djjtf,d v utfy , 



9) 



if we neglect the infinite tower of mesons and pick up only the pion field. A parameter e which is 
0(y/Nc) is dimensionless and takes a real value. The second term (sometimes called the Skyrme term) 
is absent in the leading-order ChPT, and the presence of this term is crucially important; the soliton 
solution would collapse to zero size without this fourth-order derivative term (Derrick's theorem; this 
might be relaxed with inclusion of u meson, see Ref. |181j ). 

Introducing a current, S M = Ud^U\ the expectation of the Hamiltonian leads to the energy, 



E 



d xtr 



J 71 



16e : 



(90) 



and the baryon current is expressed as jg = (l/247r 2 )e / ^ a/3 tr£„£ Q ,£ / g, from which one can define the 
baryon number, B = J d 3 xj^(x). With the following hedgehog Ansatz using a spatial unit vector 
r = r/r for a chiral soliton, 

U{x) = exp[zF(r)r ■ r] , (91) 
the energy can be written functional of F, i.e. 

1 sin 2 F /sin 2 F 



E[F] = 4tt 



dr r 



fl( l9 sin 2 F 
•'" ' F' 2 + 2 



+ 



2e 2 r 2 



+ 2F" 



and the baryon number density is calculated into a form of j^(x) = — F' sin 2 F/(2n 2 r 2 ) 
number is then properly quantized as 

B = — [2F(0) - sin2F(0)] 1 , 
if the boundary condition is chosen as F(0) = ir, where we have used i^oo) — > 0. 



(92) 
The baryon 

(93) 



The concrete form of the pion profile F(r) is determined to minimize the energy (92), 

sin2F sin 2 Fsin2F 



P 2 



+ 2 sin 2 F )F" + -F' + F sin 2F 



. 



(94) 



Here we used a dimensionless variable r = r/R with = 2ef n . The numerical solution that satisfies 
the proper boundary conditions is plotted in the left panel of Fig. 11 For a more realistic description, 
one should include strangeness, quantize the solution to give proper quantum numbers such as the 
angular momentum and the iso-spin, which is beyond our scope of this review. 
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Figure 11: (Left) Profile function F(f) as a function of the dimensionless radius variable. The slope at 
the origin is chosen to satisfy the boundary condition; F(f — > oo) — > 0. (Right) Schematic picture of 
nuclear matter of skyrmions placed on S 3 . 



Nuclear matter as a Skyrme crystal This pion profile describes a single baryon, and it is neces- 
sary to extend the formulation to be capable of including multi-baryons to study the nuclear matter 
properties. It is not easy to construct skyrmions for B 3> 1 directly, but a system at finite density can 
be emulated with periodically placed skyrmions [182] since baryons are infinitely heavy and static in 
the large- N c limit. One can analyze this problem by putting a skyrmion on S 3 , which can be viewed 



as mult i- skyrmions on R 3 as schematically visualized in the right panel of Fig. 11 This idea has been 
extended to the skyrmion in the holographic QCD model |183j . 

Interestingly enough, such a description of skyrmion matter might suggest that nuclear matter in 
the large- N c limit could have an inhomogeneous structure in it. Actually, the Skyrme model is based on 
the non-linear representation of the chiral model as we saw above, and it is impossible to reach chiral 
restoration, and yet, one can argue chiral restoration as a spatial average of (till). This implies that the 
chiral condensate may have inhomogeneity there, and this is indeed one of the characteristic features 



in quarkyonic matter as we see in Sec. 4.5 



Half-skyrmions and the Brown-Rho scaling In Ref. |124] a universal scaling law in dense matter 
was suggested; 

Hp) = ^ = ^ = f , (95) 

m N my f n 

where the asterisks are for in-medium quantities and $ is a scaling function of density p. This scaling 
called the Brown-Rho (BR) scaling is achieved starting with the chiral Lagrangian with the Skyrme 
term and modifying it such that the correct scaling property expected from QCD is recovered. 

The periodic skyrmion configuration is not necessarily optimal in energy and more careful analysis 
on crystallography should be necessary, which results in BCC and FCC crystal structures involving half- 
skyrmions |184tll85] . In fact, an intermediate phase where a skyrmion turns into two half-skyrmions was 
found at some critical density pi/ 2 higher than the normal nuclear density p 118b] I187j . is also observed 
in the holographic QCD models |188j . For an intriguing implication to the QCD phase diagram, the state 
with half-skyrmions indicates the presence of the pseudo-gap phase above the chiral phase transition as 
is suggested in Ref. [189J. 

It turned out that the meson and baryon masses behave differently in increasing density and the 



scaling (95) holds fairly well up to the density py 2 , but above pi/ 2 , it is significantly deviated |19Uj . The 



baryon mass drops at a slower rate than the meson mass with increasing density and stays non-vanishing 



at chiral symmetry restoration, which matches the chiral invariant mass mo discussed in Sec. 3.3.1 Its 
consequences on the tensor forces and the EoS at high density are discussed adopting a modified scaling 
motivated by the above observation in Ref. |190j . 
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Figure 12: Geometries in the confined and the deconfined phases. The roles of the (imaginary) time r 
and the extra coordinate £4 are swapped from one to the other phase. 

3.5.3 Sakai-Sugimoto model 

Another non-perturbative approach developed in the large- N c limit is the holographic technique that 
is becoming more and more familiar among QCD physicists nowadays. There are countless number of 
works using the holographic QCD models, and it is not realistic to make an attempt to cover them all. 
We will pay special attention to the phase diagram disclosed by the Sakai-Sugimoto model, following 
the computational steps of Ref . [22] . 

The basic idea of the holographic approach is based on a hypothetical correspondence between a 
gravity theory in higher dimensions and a field theory of our interest (see Ref. |191j for a review). The 
original conjecture was made for the M = 4 supersymmetric Yang-Mills theory, which at large- N c and 
large- 't Hooft coupling may be equivalently described by the classical solution (anti-de Sitter metric) 
of the super-gravity theory. For the investigation of QCD, however, supersymmetry and conformal 
symmetry are unwanted and should be gotten rid of. To this end, one can compactify one extra direc- 
tion and impose the anti-periodic boundary condition to fermionic super-particles, so that unphysical 



particles become heavy and decouple from the low-lying dynamics (see Fig. 12). This extra direction is 
denoted by the coordinate X 4 here. 

In this subsection, to simplify the notation, we use dimensionless variables rescaled by the radius of 
the AdS space, R. That is, 

I 1 X ~ X 4 U 

R' X = R' X4= ^R' U= R' ^ ' 

where the first (t, x) refer to the ordinary Minkowski coordinates, A 4 is the compact direction, and u, 
together with (t,x), spans the 5-dimensional AdS space, and the Minkowski space-time resides at the 
UV edge, u = 00. In what follows, we omit writing the tilde for notational simplicity. 

Deconfinement phase transition In the confined phase at low temperature, the bulk geometry 
corresponding to the D4-brane background is expressed by the following metric, 

du 2 

(Confining geometry) ds 2 = u 3/2 \-dt 2 + dx 2 + fiu)dxV\ + -. . + u 1/2 dQ 2 4 , (97) 

where f(u) = 1 — (mkk/«) 3 and mkk is fixed by the size of X4 compactification. This above geometry 
is singular at u = mkk and to avoid the conical singularity at u — mkk, i n the same way as the angle 
variable in the polar coordinates, the period of £4 (denoted by 5x 4 here) is uniquely determined that is 
translated to a mass scale called the Kaluza-Klein mass, i.e. 

4tt 2ixR~ l 3 , , , . 

5x A = => M KK = = t;V^ r ■ ( 98 ) 

3m k 7 k 0X4 l 

The physical meaning of Mkk is a cutoff scale above which super-particles could get excited and the 
Sakai-Sugimoto model would be no longer a QCD dual. As we will see soon later, the deconfinement 
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transition temperature is smaller than Mkk by a factor and the dual description of the phase diagram 
is indeed closed within the validity region. 

In the deconfined phase at finite T, on the other hand, the so-called AdS-blackhole solution should 
be energetically preferred with the metric, 



(Deconfining geometry) ds^ 



u 



3/2 r_ 



[-f(u)dt 2 + dx 2 + dx 2 ] + 



du 2 



u 3 / 2 f{u] 



+ u 1 ' 2 dn\ 



4 • 



(99) 



where the singularity in 24 is replaced in the temporal direction, 
function f(u) is changed to f(u) = 1 — (ut/u) 3 with ut fixed by the period as in Eq. p8_ 
space-time with t — it the period of r, that is nothing but the temperature inverse 1/T, leads to 



In Eq. (99) the definition of the 
In Euclidean 



1 

T 



3u 



1/2 

T 



(100) 



Therefore, the location of the horizon at u = Ut corresponds to the physical temperature T. 

One can then evaluate the free energy to identify which of the confining and the deconfining metric 
is more favored and what the critical temperature is. We can easily make a theoretical guess for 



Ta from the corresponding geometries illustrated in Fig. 12 Obviously the energies (or actions) should 
be identical when two periods coincide. That is, the first-order phase transition (Hawking-Page-type 
transition) should take place |192j when 



6x4 



1 

T 



rr 3 ! M KK 



(101) 



On the phenomenological level Mkk is fixed to be 949 MeV to reproduce the vector meson spectrum, 
so that Tj ~ 151 MeV is obtained, which is of the same order as the empirical value presented in 



Sec. 2.3 In the leading order of the large- N c approximation this critical temperature for deconfinement 



is insensitive to the presence of matter and, as discussed in Sec. |3.5.1[ this is simply because quarks are 
suppressed by l/N c as compared to gluons. 

We here make a remark on the weak-coupling extrapolation with non-trivial /3 function. This may 
lead to a possibility of disconnection between deconfinement in the strong-coupling and the weak- 
coupling regimes |193| . 



Flavor branes Let us now turn to chiral symmetry realized in the Sakai-Sugimoto model [23] and the 
chiral phase transition. The same chiral symmetry as QCD, i.e. UiX-ZVf) x Ur(A^), can be emulated by 
putting Nf D8-branes localized at £4 = and Nf D8-branes at x^ = L. This L is a parameter inherent 
in the Sakai-Sugimoto model and its physical interpretation in QCD is still unclear. For the meaning 
of the lower edge u c , see Fig. 13. We will explain how to fix u c when we incorporate the density source 



next. In the large- iV c limit the D4-brane is dominant and one can assume that the classical solutions 
are intact regardless of flavor branes (that is, the probe approximation), which corresponds to quark 
suppression by 1/N C . 

The induced metric on the D8-branes is (apart from the spherical angular part, f^), thus, 



(Confining geometry) ds 2 = i? /2 (-dt 2 + dx 2 ) + 



u 3 / 2 f(u)x' 4 (uf + 



1 



(Deconfining geometry) ds 2 



u 



3/2/ 



(-f(u)dt 2 + dx 2 ) + 



u^x'^uf + 



u 3 / 2 f(u) 

1 

u 3 / 2 f(u) 



du 2 , 
du 2 . 



(102) 
(103) 



Here x±(u) is solved with the equation of motion later. To introduce a finite density, we need a gauge field 
Aq(u) that becomes the quark chemical potential at u — 00. We can then express the Dirac-Born-Infeld 
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Figure 13: Flavor branes (D8 and D8 separated by the length parameter L) in the confined and the 
deconfined phases. In the confined phase chiral symmetry is inevitably broken, while in the deconfined 
phase both chiral broken and symmetric states are realized. 

(DBI) action involving the gauge fields on the Nf D8-brane, after integrating over the four-dimensional 
angular part f2 4 , as 

S™ = ! dt d 3 x du u 1/A J- det{g aj3 + F^) . (104) 

V 4 J " 

Here the dilaton potential is already included in the action and the field strength tensor in the above 
is rescaled to eliminate the string scale l s (see Refs. [231 122] for details). The normalization constant 
including the fi 4 integration is given as M = V 4 ■ (N C N { /3)R 6 /((2it) 5 I®) with V4 = / dtd 3 x. 

In Euclidean space-time at finite T, integrating further over the Euclidean coordinates, we can 
express the action in a form of the one-dimensional integration with respect to u. Hereafter we drop the 
irrelevant normalization M from the action and then the action simplifies, respectively, in the confined 
and the deconfined phases as 



(Confining geometry) S™ = J duu 4 ^ f(u)x' 4 (u) 2 + ^ (/(u)" 1 - a' 2 ) , (105) 



(Deconfining geometry) S^f 1 = J duu A \j f(u)x' 4 (u) 2 + ^(1 — a' 2 ) , (106) 

where ao(u) denotes the rescaled A (u) that eventually translates to the (dimensionless) chemical po- 
tential. It is convenient to use a density variable p(u) = —5S^f l /5a' Q (u), because p(u) turns out to be 
^-independent thanks to the equation of motion. Then, we can easily solve a' (u) from the definition of 
p to find, 



r ■ . , /. . ju 3 f(u) 2 x' 4 (u) 2 + 1 

Confining geometry) a (u) = p \ — — ^— , (107) 

j{u){u + p A ) 



i \ u 3 f(u) xAu) + 1 

(Deconfining geometry) a o\ u ) — P \ s s ; (105 

u b + p 2 



from which the quark chemical potential /i q = ao(oo) is obtained in respective phases. 

Introduction of the density source Here we shall explain how to fix the lower boundary u c dis- 



played on Fig. [13} If the system is in the deconfined and chiral symmetric phase at high T, D8 and 
D8 are parallel and such a configuration can have a finite density without source, and thus u c is not 
necessary and the u-integration starts simply from ut- 

When chiral symmetry is spontaneously broken, on the other hand, there should be a source for the 
density coming from the Chern-Simons action. To accommodate a finite density of baryonic matter, 
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instantons or the D4-branes wrapped on the S A inside of the D8-branes should be placed, as argued in 



Ref. |194j . in a similar way as the skyrmion in Sec. 3.5.2 In the confined phase baryons are the unique 
source for the density, while quark matter is also possible in the deconfined phase. Such a density source 
originating from quarks can be formulated as strings stretched from the D8-branes to the horizon. It 
is known, however, that such a density source leads to a larger action (i.e. larger free energy) than the 
instanton source [22], which means that baryonic matter is always favored than quark matter in this 
model. This observation might be given a deeper insight from the concept of quarkyonic matter that is 



the main subject of Sec. 4.5 



The point where the density source is placed along the w-direction defines u c and its value is dynam- 
ically determined by the zero-force condition. In the presence of the source the flavor-branes are pulled 
and x&(u) has a cusp at u — u c , as illustrated in Fig. 13, and the range of the u-integration should 
be from u c to +oo. The zero-force condition describes a situation that the tension of the D8-branes 
/ds should be balanced by the force from the source D4-brane /d4- These forces can be given by the 
derivative of corresponding actions with respect to u c . 



Deconfined and chiral symmetric phase On the phase diagram we will find a first-order phase 
transition associated with the restoration of chiral symmetry, and then the density p jumps discon- 
tinuously. So, to distinguish the densities in respective phases, let us denote the density in the chiral 
symmetric phase as p sym , while we use Pbroken in the chiral broken phase. 

When chiral symmetry is restored at sufficiently high temperature, D8 and D8 are straight and 



parallel (the far left of Fig. 13). In this case the solution of the equation of motion is as simple as 

ut can take care 
- a (oo), that is, 

(109) 



x' 4 (u) = 0. Then, we do not have to introduce a source term but the horizon at u - 
of the role as a density source. The physical chemical potential is obtained from p q 



Pq = P: 
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du — ^=^= 

UT + P~ 



sym 



from Eq. (108). The thermodynamic potential is then calculated from the action (106) as 
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Deconfined and chiral broken phase The treatment of the chiral broken phase is much more com- 
plicated than the chiral symmetric phase because the configuration of D8-branes should be numerically 



solved from the Euler-Lagrange equation from the action (|106|) as 
x' 4 (u) 



\//K)K + Pbroken«o) 



111) 



U*l 2 ^f(u) [f(u)(u* + pLkcn^ 3 ) - /MK + Pbroken^)] 

where Uq is defined by a point where x' 4 (uq) diverges. Then, in this phase, the chemical potential is 



expressed from Eq. (108) as 



pq Pbroken 



du 



f(u)u 3 



U r 



/(«)(« 8 +PLken« 3 )-/MK + Pb 



2 „.3) 
roken"0J 



;ii2) 



We note that the second term in the right-hand side is added by the D4-brane contribution in the 
presence of density, and this determines the threshold p q for finite density. Such an onset of the density 
corresponds to wn — B/A in nuclear matter; we have already seen p = as long as ps < — B/A (see 
the right panel of Fig. [7]). In other words, if p q is smaller than this threshold value, the only possible 
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solution is Pbroken = and the ground state is matter-free empty. We should concretely figure out the 
lower edge of the integration, u c , to proceed in the calculation. The zero-force condition, / D8 = / D4 , 
yields [22], 

l f\u c )(ul + P 2 M 3) _ f{uo){u 8 + p 2 M 3) = 

8Ur 



P 



f(Uc)u 



u, 



V7( 



Uc 



;ii3) 



The thermodynamic potential in the chiral broken phase is, thus, obtained from the action, which 
amounts to 



n 



broken 



du 



[U)U 3 U 5 



Pbroken^ 3 ) 



Pbrokcn^o) 



(114) 



The above-mentioned setup suffices for the determination of the chiral phase transition point by 



the comparison between Eqs. (110) and (114). The concrete procedure for this is as follows: First, for 
a given density pbroken, the zero-force condition (113) gives u c (u ,T) 

21 



(u)du expressed clS db function of u . 



and this renders the condition 
This enables us to solve «o(pbroken, L, T), so that 



Eq. (112J) directly relates /i q and pbroken- Then we get fibroken(£, T, p q ) finally. In the same (and much 
simpler) way, we can calculate fl sym (L, T, p q ), and we can locate the first-order Hawking-Page transition 
of chiral restoration by looking for T and p q that makes ^broken = ^sym- in this way, particularly at 
vanishing p q , the critical temperature turns out to be T x = 0.154/L, which was first found in Ref. [195J. 



Density onset and the phase diagram The phase diagram is not yet complete. We should consider 
another line corresponding to the liquid-gas phase transition as presented explicitly in the conventional 
description in the right panel of Fig. [7j 

In the confined phase chiral symmetry is always broken. The equation of motion for x' A (u) is slightly 



different from Eq. (Ill) and is given by 



xAu 



f(u )(u S + Pbroken M o) 



U*f(unf(u)(u8 + pLken^) " /(«0)K + Pbroken^)] 



IIS) 



from which the zero-force condition results in a similar form to Eq. (113) in which the right-hand side 



is replaced with d(pu c /3)du c . The onset at which a finite density is turned on can be found analytically 
with pbroken = plugged into the above equations. Then, we simply have u c = u , and the corresponding 
chemical potential is p q = u /3. At the onset we can approximate (115) with f(u) ~ 1 to derive u 
analytically, and the critical p q is then obtained as (47r/3)(r[9/16]/r[l/16]) 2 ~ 0.175 in the unit of 
L. In the deconfined phase, on the other hand, the density onset lies at p q = (u /3)a//(mo), which is 
T-dependent. Then, the phase boundary is not a straight vertical line but it has a minor dependence 
on T. 

Some numerical calculations are needed to clarify the whole phase structure, and the resulting 

which first appeared in Ref. [22]. This is a quite suggestive 



phase diagram is depicted in Fig. 14 



phase diagram, but we shall defer a more detailed comparison between Fig. 14 and the expected phase 



structure of large- iV c QCD till Sec. 4.5 



We here raise some caveats; the above treatment is based on the probe approximation that corre- 
sponds to the quench limit of QCD, and the 1/N C corrections should be considered to take account 
of back-reactions from matter. For more details including other holographic models than the Sakai- 
Sugimoto model, see a recent review |196] . This is not yet the end of the story about the phase diagram 
in the Sakai-Sugimoto model. Recently there has been an important recognition; an instability of gauge 
fluctuations at non-zero momenta was discovered at finite density [1971 1198] . and so, the phase diagram 



in Fig. 14 must be supplemented with additional phase boundaries associated with spatial modulations 
(see Fig. 1 of Ref. |197] for instance). 
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Figure 14: Phase diagram at finite temperature and density in the Sakai-Sugimoto model. The decon- 
finement temperature is chosen at Td = 0.08 and the scale is set by L = 1 according to the choice of 
Ref . [22J . The vertical dashed line represents the onset of baryon density. 

It is quite interesting to notice that this sort of spatially inhomogeneous structure was once discussed 
actively in the context of nuclear physics, and nowadays, is becoming one of the central issues in quark 
matter research. Here, for a while, it should be worth revisiting the possibility of inhomogeneous 
condensation in nuclear matter. 

3.6 Inhomogeneity: pion condensation 

In a medium at finite density there are collective excitations made from particle and hole (p-h) in the 
same channel as pions. The interaction between N and ir is, however, repulsive in the s-wave channel 
and the pion self-energy is positive then. It was recognized later in Refs. |199[ 1200] that the p-wave 
interaction is attractive, which causes a resonant state A in the channel of L = 1, J = 3/2, T = 3/2, 
and it would render the pion energy decrease in matter. This attractive interaction is attributed to 
underlying chiral symmetry. 

The in-medium pion propagator can be written with the self-energy as 

D-\u,p) = u 2 -p 2 -ml-U(u,p) , (116) 

where the self-energy U(u,p) should involve the p-h and A-h contributions; II = — p 2 [Un(u,p) + 
(7a(w,p)]. In the symmetric nuclear matter with Z = N that is of our main interest, the collective 
mode is pushed down at finite p and there appears a condensate of n°, while in neutron matter which 
is relevant to the astrophysical application the collective ir + (Migdal's 7r+) and 7r~ are spontaneously 
generated [200] . 

When the ir° condensation occurs in symmetric nuclear matter at the momentum p c , the expectation 
value should behave as 

(n°(x,p c )) ~ cos(p c • x) . (117) 

The physical implication of such condensate to nuclear matter has been discussed (see Ref. |201j and 
other contributions in this volume). 

One can find the critical density for the 7r° condensation, for instance, adopting the linear sigma 
model including p, n, and A's. Once the density exceeds the threshold, however, the system is unstable 
because only the attractive force is taken into account and the short-range (shorter than the pion 
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Compton length) repulsive interaction between baryons should be considered. In the Fermi liquid 
theory such short-range NN interactions (proportional to 8{xi — x 2 )) are incorporated through the 
Landau parameters in the following form; 



NN 



fNN + f' N N T l ■ T 2 +gNN(7l ■ &2 + ( IdNN^l * ^(ji • T 2 ) . (118) 



Here T\ and cr\ (and r 2 and cr 2 ) refer to the isospin and the spin of nucleon 1 (and nucleon 2, respec- 
tively) in the NN system. In fact this last term is required to cancel Dirac's delta function appearing 
in the one-pion exchange potential (OPEP) (if g' NN = 1/3 is chosen, the attractive delta function in the 
OPEP is exactly canceled). Such a cancellation is a universal phenomenon; the Lorentz-Lorenz effect in 
dielectric media is the most famous example. In the same manner, including the mixing with isobars, 
one should further consider the Landau-Migdal parameters as 
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9na (""i " <)(rt ■ rf ) , F AA = ^ ) 9 aa ("i • ^)(rf • r*) , (119) 



all of which contribute to Un(w,p) and U A (ui,p). Because g' NN , q'nai anc l 9' a a represent the effect of 
repulsive interactions, they would disfavor the pion condensation. Therefore, the precise determination 
of these parameters is fatally important to judge whether the pion condensation can occur or not in 
reality. 

Today it is widely believed that the possibility of the pion condensation in nuclear matter was 
somehow ruled out. This was concluded from the experimental data of the giant Gamow- Teller (GT) 
states which are collective modes with respect to the spin-isospin fluctuations, so that their strength 
can constrain g' NN , g' NA , and g' AA (see Ref. [202] for a review). Since the decomposition is difficult, it 
used to be often the case to assume the universality, g' NN = g' NA = g AA = g', and then the experimental 
data suggests g' = 0.6 ~ 0.9, which is much larger than 1/3 and the pion condensation cannot occur 
in a reasonable range of baryon density. Later on, in theory and experiment, different values of g' have 
been investigated |203] , and as comprehensively reviewed in Ref. [204] . Although the pion condensation 
is unlikely at normal nuclear density po, the fact is that there is still a good deal of chance for the 
p-wave pion condensation to develop in nuclear matter at p > 2p - According to Ref. [204J , the critical 
density p c is calculated with g' NN = 0.59 and g' NA = 0.18 + 0.05g^ A fixed from the experimental data, 
which results in p c = (1.2 ~ 2)po for the undetermined parameter g' AA = ~ 1. This makes a sharp 
contrast to p c under the universality Ansatz, which rapidly grows up to unrealistic density. Definitely, 
the possibility of the pion condensation should deserves more investigations. As we will see later in 



Sec. 4.4, an interesting form of the inhomogeneous condensate analogous to the condensate (117) is a 



likely candidate for the ground state of quark matter. 



4 Quark/ Quarkyonic Matter 

In this section we address the properties and the phase transition in deconfined quark matter assuming 
that the ground state of high density matter is composed from quasi-quarks. First we will make a 
quick review over a phenomenological description of deconfinement, and then continue our discussions 
utilizing chiral effective models. Finally we will closely discuss the possibility of inhomogeneous chiral 
condensate and the general mechanism that would derive inhomogeneity at high density. 

4.1 Statistical bootstrap and the Hagedorn limiting temperature 

In quantum field theories we are still far from satisfactory understanding of confinement and decon- 
finement. In phenomenology, however, it is quite useful (and probably a right picture) to think of the 



45 



Hagedorn spectrum. Let us denote the density of states p(m 2 ) so that p(m 2 )dm counts the number of 
hadronic states in the mass range m ~ m + dm. The idea of the statistical bootstrap model is that 
hadronic states are composed from a convolution of smaller hadrons. This picture is embodied as the 
following statistical bootstrap equation for the hadron mass spectrum [205J, 



oo 1 „ N N 

Hp(p 2 ) = H6( Po )5(p 2 -ml) + J2m ^ " 5>> U H P^ d 'P 

N=2 i i=2 



(120) 



where H represents a proper hadronic volume and m^ n is an input parameter corresponding to the 
"elementary" mass in the bootstrap picture. Interestingly, it is possible to fix the asymptotic form 
of p(m 2 ) semi-analytically from the above equation. Introducing the Laplace transforms of Eq. (120) 
as G({3) = f e~ l3p Hp(p 2 )d 4: p and <p(/3) = f e~ l3p H9(p )5(p 2 — m 2 n )d 4 p, where j3 corresponds to the 



temperature in the Boltzmann factor, the bootstrap equation (120) is rewritten into a simple form; 



Gtf) = <p(0) + e G ^ - Gtf) 



(121) 



Although it is not easy to solve G as a function of <p, Eq. (121 ) already solves ip in terms of G; tp — 2G + 
1— e G , which has a maximum and we shall denote this extremal point by (p$, and this </2 can be converted 
to corresponding (3q. Then, G(/3) can be expanded around @o, i.e. G((3) ~ G{(3q) + (const.) ((3 — (3q) 1 / 2 . 
Such a type of square-root singularity can be reproduced by the density of states, 



p(m 2 ) — Cm 



-3 e m/T H 



(122) 



where Th is a temperature parameter (related to and called the Hagedorn temperature. This 
exponential form of the density of hadronic states is an important observation in understanding the 
physics of deconfinement at the phenomenological level. 

The physical interpretation of the Hagedorn temperature is quite instructive. The partition func- 
tion, J p(m)e~ m / T dm, would diverge when T > Th and thus the physical temperature cannot exceed 
Th as long as the hadronic description is valid. In this sense Th used to be considered as the limiting 
temperature. In Ref. [23] it was argued that Th should be the critical temperature where the hadronic 
description breaks down and the system goes through a phase transition into deconfined quarks. Ac- 
tually the very first prototype of the QCD phase diagram was conjectured in Ref. [23] based on this 
interpretation of Th- In a modern language stemming from the Hagedorn picture, the deconfinement 
phenomenon is correctly captured by the thermal statistical model or the hadron resonance gas model 



(see also Sec. 4.5.2) 



4.2 Chiral effective models 

Once quark deconfinement occurs, the most suited description of the state of matter should be of 
interacting quarks rather than mesons and baryons. It is still challenging to tackle the QCD phase 
transitions from the first principle unless the the quark density is so small that the lattice-QCD is at 
work. Then, for the pragmatic studies on the QCD phase diagram, one can adopt an effective model 
that shares common physical features with QCD, and chiral symmetry is the guiding principle for the 
model building. 



4.2.1 Nambu— Jona-Lasinio model 

At sufficiently high density it is conceivable that the physical degrees of freedom should be quarks rather 
than baryons, and then it makes sense to consider a picture of quasi-quarks as a result of integrating 
out heavy gluonic degrees of freedom. A recent attempt along this line is found in Ref. [206J. 
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In this context the Nambu-Jona-Lasinio (NJL) model is one of the most preferable choices of the 
low-lying QCD models (for reviews, see Refs. [207, I2U8] ). In principle, the integration over gluons would 
lead to an infinite number of interactions in terms of quark fields, which can be expanded according to 
the power of ip and ip as 

C = ${i llx &> - m Q )i; + V) + 42 + " - • , (123) 

where C^ t represents the interaction terms at the quartic order in terms of quark fields, 42 at the sixth 
order, etc. In the vicinity of the chiral phase transition only 42 is predominantly important in the 
mean-field analysis (if it is of second order and the mean- field mass is vanishingly small), and we will 
keep this leading-order interaction only in this article (see Ref. [209] for an extension with higher-order 
interactions). If the above £ is designed to mimic the QCD dynamics, C should be invariant under 
the symmetries of the QCD Lagrangian as exposited in Sec. [2| Motivated by the presence of the chiral 
condensate in the scalar channel, we should retain the interaction involving the scalar channel in 42 > 
and naturally, 42 should contain the interaction in the pseudo-scalar channel as the chiral partner, 
which leads to 

42 = | [(^) 2 + (^l^) 2 ] - g v (^) 2 ■ (124) 

Here the scalar part (ip\ a ip) 2 and the pseudo-scalar part (^i^^Xa^) 2 are both necessary to be consistent 
with chiral symmetry, whilst the second term (ipj^if)) 2 is chiral symmetric as it is. In the NJL model 
study in the vacuum, this vector interaction has only a minor effect except for the vector mesons, but 
at finite density we must not neglect it. As a matter of fact, this vector interaction plays an equivalent 



role as the term (m 2 J /2)u fl u 11 in the Walecka model (45). 



In QCD at the quantum level U(1)a symmetry is explicitly broken due to the axial anomaly, which 



should give rise to the instanton- induced interaction, as already derived in Eq. (28). Because the deter- 
minant in flavor space should be taken, the instanton-induced interaction is the four-quark interaction 
for Nf = 2 and the six-quark interaction for Nf — 3, and written then as 

42 =^[det^(l- 75 )V' + det^(l + 75)V] • (125) 

The model parameters, g s , g v , ga, and the UV cutoff A are fixed so as to reproduce the hadron properties 
in the vacuum. In general these parameter should vary with changing T and u q , but there is no reliable 
strategy for the determination of such medium dependence. Thus, in the model studies, one usually 
assumes that one can keep using the same parameters as fixed at T = /i q = 0, which brings in the most 
unmanageable uncertainty in theory. 

Mean-field approximation in the vacuum The thermodynamic potential is composed from three 
individual contributions in the simplest mean-field approximation of the NJL model. That is, 

f2 f^cond ^zcro ^quark j (126) 

where the three contributions come from the condensation energy, the zero-point oscillation energy, and 
the quasi-quark partition function. Let us take a look at respective contributions in order below. 

The condensation energy comes from the vacuum expectation value of the quark composite operator. 
The four-point interaction vertex of w-quarks, for example, can be approximated as 

(uu) 2 = (uu — (uu) + (uu)) 2 ~ (uu) + 2(uu)(uu — (uu)) = —(uu) 2 + 2(uu)uu . (127) 

Here we have ignored second-order infinitesimals, (uu — (uu)) 2 . Then, in the final expression above, 
— (uu) 2 is the condensation energy and the second term amounts to the dynamical mass term for u- 
quarks, which should be taken into account in the quasi-quark partition function as the constituent 
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quark mass. It should be mentioned that there is no mixing between different flavors like uudd from 
Cfol but the flavor mixing interaction appears solely from the anomaly terms; for N c = 3 for instance, 

3, 



the determinant interaction in Eq. (125) causes the flavor mixing. Including the terms f° r ^Vf 



the condensation energy amounts in total to 



n 



cond 



g s ((uu} 2 + (dd) 2 + (ss) 2 ) + Ag d (uu) (dd) (ss) . 



'128) 



In the NJL model the zero-point oscillation energy is essential to bring about the spontaneous 
breaking of chiral symmetry. Because the NJL model is non-renormalizable, the momentum integration 
needs a UV cutoff A. In a naive momentum cutoff scheme the diverging zero-point oscillation energy is 
expressed as 

" A d 3 p Ei{p) 



(27T) 



(129) 



where the overall factor 2 • 2 is from the spin and the particle and anti-particle degeneracy, and N c = 3 
is the number of colors. For a more sophisticated treatment with non-local interactions that tame the 
UV divergence, see Refs. |210[ 12111 12121 1213] . We also note that the non-local formulation is already 
close to the Dyson-Schwinger approach [214J. The quasi-particle energy e$(p) depends on the flavor 
which is indicated by the index i. In the mean-field approximation, together with the second term of 



Eq. (127), the quasi-quarks should have the constituent quark masses; 



Mi = mi- 2g a (if> i if> i ) - gd^ijk(^jipj)(^kipk) 



(130) 



It is possible to carry out the momentum integration analytically. In particular in the chiral limit (where 
m u = rrid = m s = 0), an expanded form suffices for the investigation in the vicinity of chiral restoration; 



a 



N C N{A 4 



N C N { A 4 
4tt 2 



(i + a + o(e 



(131) 



where we have defined a dimensionless variable £ = M q /A with M q = M u = = M s in the chiral 
limit and 0(£ 4 ) denotes the neglected higher-order terms. The logarithmic term has a special effect 
on the analysis of the order of the phase transition. Together with the condensation energy, the total 
energy is expanded up to the quadratic order as 



cond 



N c NfA 4 
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N c N i9s A 2 



(132) 



We note here that the instanton-induced interaction of 0(£ 3 ) is negligible in the present approximation. 
Then, if the potential curvature is negative, £ or the constituent quark mass spontaneously has a finite 
expectation value, and the symmetry breaking condition reads }215j . 



* A 2 



g s A 2 > g* s A 



NcNt 



(133) 



The existence of the critical value of the coupling strength is presumably an artifact of non-confining 
model. 



The physical interpretation of this condition (133) from the RG flow is quite informative (for a 
review, see Ref. |216j ). For this purpose let us introduce the dimensionless coupling g s ^ = g s ^k 2 . The 
coupling constant in the NJL Lagrangian should be regarded as the bare one at the UV scale, i.e. 
gs = <7s,A at A; = A. The function of g s ^ (denoted by /3 g here) behaves as shown schematically in the 



left panel of Fig. 15 If the dimensionless NJL coupling g s is smaller than the fixed-point g* that makes 
(3 g vanish, j3 g is positive and g s ^ goes smaller with decreasing k. If g s is larger than g*, on the other 
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Figure 15: (Left) Schematic picture of the RG flow of the coupling constant for the four-fermion 
interaction g s . The /3 function crosses zero at g*. (Right) Bubble resummation for the four-fermion 
coupling g s . 



hand, g S} k gets larger and eventually diverges at k — > 0. This diverging behavior of g S fi signals that chiral 
symmetry should be spontaneously broken. As a matter of fact, the RG flow is generated through the 
inclusion of fluctuations and g Sj o should contain all higher-loop diagrams in the four-fermion channel. 
In the one-loop order of /3 g , after integrating with respect to k, g Sj o resums all bubble-type diagrams as 
sketched in the right panel of Fig. 15 Thus, there appear a and ir mesons in the intermediate states 
(dotted line in the right panel of Fig. 15) in g St0 , and g s ^ — > oo at k — > exactly corresponds to the 



(zero-momentum) propagation of massless mesons at the critical point. In this way, g*A 2 in Eq. (133) 
can be given a clear interpretation. 

With gauge field fluctuations, this {3 g is pushed down and if the gauge coupling is large enough, 
flg would be entirely negative, and chiral symmetry is spontaneously broken regardless of the starting 
point of the RG flow |216j. At finite temperate, on the other hand, (3 g is pushed up, and g* moves 
toward a larger value. Then, when g* becomes greater than a given g s at a certain temperature, the 
second-order phase transition of chiral restoration occurs, as we see below. 



Phase diagram in the mean-field approximation At finite temperature and density the partition 
function of the quasi-particle excitations should be added to the thermodynamics potential, i.e. 



n 



quark 



-2N C TJ2 j (^{ ln [ 1 + e_(£<(p) "' iq)/T ] +ln[l + e-^^ + ^/ T ]| , (134) 



which induces a phase transition. In the chiral limit, again, we can find the critical temperature in 
the analytical way. That is, near the critical point, we can assume T 3> M q without loss of generality. 
Then, the high-T expansion leads to, 



N C N{A 4 

4tT 2 
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~ N c N {9s A* 3" I A, 



/fqV 
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(135) 



The potential curvature has a modification from the finite-T effects, and from the condition for the 
curvature to vanish at the critical point, the critical temperature is deduced as 



'3(A 2 



7T Z 



N c N igs 



(136) 



This analytical result leads to a typical shape of the phase diagram in the /x q -T plane as depicted 



in Fig. 16, where the Hatsuda-Kunihiro choice of the model parameters is used; A = 631 MeV and 
g s = 0.214 fm 2 for N { = 2. 

So far, our discussions are limited to the special case of the second-order phase transition, and a 
possibility of the first-order phase transition cannot be ruled out from the above analysis. In fact, 
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Figure 16: Typical shape of the phase diagram in the NJL model. The curve shows the analytical 
result (pi with A = 631 MeV and ^ s = 0.24 fm 2 used for the N { = 2 case. 



the phase diagram like the one presented in Fig. [16] should be more carefully examined and could be 
modified by a first-order phase boundary in the high-density region. This possibility and the underlying 



mechanism will be the subject of Sec. 4.3 



Coupling with the Polyakov loop One of the major obstacles in the application of the NJL model 
to the finite-T study is that this model may allow for unphysical quark excitations even below T c . Such 
unphysical excitations should be diminished by the effect of confinement in reality. In other words, 
the thermal weight e~^^ p - )=FMq - ) / T for quark excitations in Eq. (134) should be suppressed whenever the 



Polyakov loop $ takes a small value. This can be easily implemented by the alteration of Eq. (134) 
into [2171 1218] 

n quark ^-2Tj2 J (0i{^ln[l + Le-^W-^ T ] + tr ln[l + LV^H^] } . 



(27r) 3 I 

+ In [l + 3f e-feOH/^)/ 3, + 3£ e - 2 (^(p)+^)/ T + e -H^(p)+^i)/ T ^ j . (137) 



In fact, in the field-theoretical derivation, one can arrive at the above coupling of Eq. (137) assuming 



that the Dirac operator has a constant A4 background |218j . The integration along the temporal 
direction results in the Wilson line and thus the Polyakov loop as defined in Eq. (J6|. 

This analytical structure itself is quite useful to understand some features of hot and dense QCD, 
especially the sign problem [2H]. One can also find a reason why the strange quark (or heavy quark 
generally) number susceptibility can serve as a deconfinement order parameter. If the quark mass of 



the heavy flavor is sufficiently large, we can expand Eq. (137) and keep only the lowest-order term 



with respect to small Boltzmann weight. Then, differentiating the partition function with respect to 
the heavy-flavor chemical potential /x s (which is eventually taken to be zero), we can derive the heavy- 
quark number n s oc (I — t) e _Ms//T and the susceptibility % s oc (£ + £*) e~ Ma ' T . Since Xs is proportional 
to the Polyakov loop, Xs is qualified as an order parameter, which explains the right panel of Fig. [5] 

To make concrete calculations, one needs to specify the effective potential for <3>. Such a potential, 
^giue[$], is calculable in the perturbation theory that leads to what is called the Weiss potential [219, 
220J. Non-perturbative evaluation is necessary, however, to describe a phase transition. The following 
Ansatz is well motivated from the strong-coupling expansion on the lattice; 



Q g i ue [$] 



b(T) In 1 - 6$$ + 4($ 3 + $ 3 ) - 3($$) 2 



138) 
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where the logarithmic term has an origin in the Haar measure or the Jacobian associated with the 



variable change from L to A 4 . (The potential (138) is a function of the traced Polyakov loop that can 



be expressed using A4, and thus the Jacobian from the group integration of L to A 4 appears.) This 
Haar measure potential has cubic terms reflecting the color SU(3) nature. The T-dependent coefficients 
are chosen so that this fi g i ue [$] can reproduce the thermodynamics of the pure Yang-Mills theory. One 
might wonder that $ or A 4 cannot saturate the full thermodynamics and physical gluons should be 
rather two A^s (transverse components). The point is that Ap's also couple to the Polyakov loop in 
the color adjoint representation and their thermal excitations can be controlled fully as a function of 
$ [2211 [222J, which is also manifest in the non-perturbative construction of f2 g i ue [$] by means of the 
inverted Weiss potential [22311224] . 

Suppose that f2 g i ue [<3>] is known, the Polyakov loop expectation value and the pressure are read from 



9fi glue [4>,$] 
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(139) 



$=$ 0l <j>=ij> 



(1 



We note that $ = (£*) is independent of $ at finite density. For /x q > generally $ > $ is energetically 
favored. Once the expectation values $0 and $0 are known, one can evaluate all other thermodynamic 



quantities from f2 g i ue - For the parametrization of a(T) and b(T) in Eq. (138) the following choice was 
inspired from the strong coupling expansion [225], 



a(T) = 12N*(3TA 3 e 



-aA/T 



b(T) = -(3TA 3 



(140) 



where a = 1.052 and j3 = 0.03 are fixed to set the crossover temperature around T = 150 ~ 200 MeV. 
Another frequently-used parametrization [226] is motivated to fit the thermodynamics of the pure Yang- 
Mills theory as 
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with T being a scale parameter that characterizes the critical (crossover) temperature. Blue curves in 
Fig. [3] show a quantitative comparison between the lattice data and the potential fit from this fl g \ ue - The 
agreement is pretty good in the entire temperature range. One might have thought that the agreement 
is simply a result of the fit, but it is a highly non-trivial question whether the entire T range can be 
covered with only three fitting parameters. 



It is a straightforward extension of the NJL model with this modified fi qua rk of Eq. (137) and the 
Polyakov loop potential fi glue of Eq. (138). Such a model setup is called the Polyakov loop extended 
Nambu-Jona-Lasinio (PNJL) model (there are hundreds of application works; see e.g. Refs. [227] . Here 
let us make a technical remark; we still need to introduce a further approximation in the PNJL model 
for the treatment of Eq. (137). In many calculations the traced Polyakov loops are replaced as £ — > $ 
and £* — > 3>, but this procedure assumes (£ n ) — > and (£* n ) — > $™, which can be acceptable unless 
the Polyakov loop fluctuations become significant as compared to the expectation values. It is possible 
to go beyond this approximation to make use of a matrix model of the Polyakov loop [2281 12291 1230] . 



Figure 17 shows typical results from the mean-field approximation of the (2+l)-flavor PNJL model. 



The order parameters, i.e. the Polyakov loop and the normalized chiral condensate are plotted as 
functions of fi q and T. The NJL model parameters are fixed according to the Hatsuda-Kunihiro choice; 
A = 631 MeV, m u = m d = 5.5 MeV, m s = 135.7 MeV, g s A 2 = 3.67, g d A 5 = -9.29. In the market there 
are several different (and consistent) choices available; for example A = 602.3 MeV, m s = 140.7 MeV, 
g s A 2 = 3.67, gdA 5 = —12.36 [231] leading to qualitatively consistent results. We here make a remark 
that the PNJL model has been applied to color-superconductivity as well |226[ 1232] , but one elaborates 
an improved treatment of the Polyakov loop [230J , otherwise the color density behaves unphysically. 
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Figure 17: Order parameters calculated in the (2+l)-flavor PNJL model in the mean-field approxima- 
tion. Figures adapted from Ref. [225J. 



We note that the trend of two simultaneous crossovers is semi-quantitatively reproduced in the 
whole fi q -T region. There is seen a discontinuous jump in the chiral condensate in the high-density 
region, which signals a first-order phase transition. Actually it is quite non-trivial what is going on in 
the high-density region. One can see that the Polyakov loop stays small at low temperatures, though 
a jump is still visible which presumably arises from the mixing with the chiral condensate. One may 
be tempted to make a conclusion from the smallness of $ that confinement must persist even in the 
high-density region. This smallness is intuitively understood from the interpretation of $ as a single- 
quark free energy f q , i.e., $ ~ e~f q / 7 '. Quark confinement means f q — > oo and so $ — > 0. However, even 
with a finite f q , we can have $ — > in the limit of T — > 0. Thus, $ is no longer a good measure for 
deconfinement at low temperature. The physical picture of confinement and deconfinement of quarks 
still has subtleties, which has inspired theoretical deliberations on quarkyonic matter. 

Finally, we make a remark on the sign problem within the framework of the PNJL model. The quasi- 
quark contribution of Eq. (137) corresponds to the logarithm of the Dirac determinant in the presence 
of constant Aq from which L and U emerge. In the color SU(3) case the traced Polyakov loops, i and 
£*, generally take a complex number. When the density is zero with yU q = 0, the Dirac determinant is 
positive definite. At /i q ^ 0, however, fi quark can be complex, or its real-part can be negative. In the 
mean-field approximation in which trL and trL^ are replaced with $ and i>, the energy becomes real 
and the sign problem seemingly disappears. However, we can show that the thermodynamic potential 
$] is unstable in the direction of growing 5$ = $ — $. The solution of the gap equation turns out 
to sit at the saddle point of the potential [2 33 ] . Such weird behavior of the thermodynamic potential is 
a manifestation of the difficulty of the sign problem within the mean-field approximation. 



4.2.2 Quark-meson model 

In the (P)NJL model it is technically complicated to include the effects of meson fluctuations. In 
such microscopic approaches the meson propagator is dynamically generated through quark loops, and 
thus, the momentum dependence is highly non-trivial. Moreover, the validity of the model is severely 
restricted by the Mott threshold for unphysical meson decay into a quark and an anti-quark. 

Instead of the NJL model, depending on the problem, it would be sometimes more suitable to adopt 
another chiral model including meson degrees of freedom as point particles, namely, the linear sigma 
model. The linear sigma model usually has fermionic degrees of freedom as color-singlet baryons, but 
to investigate the phase transitions at high T and/or /i q , quarks may be able to enter the dynamics, 
and then, a model with quarks and mesons could make more sense. Such a model setup with quarks 
and mesons is called the quark-meson (QM) model. 
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Problems in the mean-field approximation In the two-flavor case, the QM model consists of 
quarks and light mesons, a and tc, with the following Lagrangian density, 



7T 
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T 



cr 



7T — f ) + CCT 



(142) 



With the potential in this Lagrangian density, chiral symmetry is spontaneously broken and the last 
term with c = f^vr^ (which is constrained by the PCAC (23)) makes the potential slightly tilted so 
that cr acquires an expectation value and three 7r's become the nearly massless NG bosons. In the 
chiral limit (c = 0) one can perform the mean-field approximation simply in an analytical way, which 
is, however, troublesome as compared to the NJL model. 

At T = yU q = the chiral symmetry breaking is almost trivial. The potential is so designed that 
a can have an expectation value at v, which turns out to be /„.. The constituent quark mass is then 
M q = go. Unlike the NJL model the zero-point oscillation energy is discarded because the potential 
terms are supposed to take care of it. Then, at finite T, chiral restoration can be caused by meson 
excitations even without quarks, and this fact that one can treat meson loop effects is the advantage 
of the QM model over the NJL model. However, it is known that the simple mean-field approximation 
(or the Hartree approximation) would lead to an unphysical first-order phase transition due to the 
approximation artifact [234J. 

The presence of quarks would make the situation even worth. It is known by now that the con- 



ventional treatment of discarding Eq. (131) is inappropriate for the finite-T investigation especially in 
the chiral limit [235] . To understand the problem, let us take a close look at Eq. (131). We have seen 



that this expanded form in terms of £ = M q /A correctly describes a second-order phase transition in 
the NJL model. Strictly speaking, one should check if the quartic term ~ £ 4 has a positive coeffi- 



cient, and in the case of Eq. (131), indeed, the coefficient of the £ 4 term has a logarithmic singularity 
~ — In | \J\ + £ 2 — 1| — > +oo, and it is always positive for sufficiently small £. We note that this loga- 
rithmic singularity appears from the IR sector of massless quarks. Interestingly enough, at T ^ 0, this 



logarithmic singularity is exactly canceled by terms from Eq. ( 134) because quarks are fermions without 



Matsubara zero-mode at finite T. Therefore, the IR singularity should disappear whenever a finite T is 



turned on. This implies that Eq. (134) alone without Eq. (131) would leave an uncanceled singularity, 
+ ln | a/I + £ 2 — 1| — > — oo, for the coefficient of the £ 4 term. So, a first-order phase transition occurs 
when the £ 2 -term approaches zero. This is obviously an artificial phase transition. 

We can overcome this problem in the following way. One easy resolution is to introduce physical 
quark masses, and then the chiral phase transition is a smooth crossover. Such a treatment might 
be useful for qualitative studies at /x q = 0, but the tendency to favor a first-order phase transition is 
overestimated, which is fatal for the QCD critical point search. For a more serious remedy, the latter 
problem is not so difficult to remove; because the logarithmic singularity is analytically known, one 
can cancel it by hand, or one can simply add the zero-point oscillation energy and readjust the model 
parameters. More cumbersome is the problem of artificial first-order transition; one should improve the 
approximation to go beyond the Hartree resummation. 



Renormalization group improvement It is a tedious calculation to handle the Fock and the 
higher-loop diagrams such as the sunset type [236]. Fortunately, nowadays, a sophisticated formulation 
for non-perturbative resummation has been developed well, which completely resolves the problem of 
artificial first-order phase transition. 

Let us elaborate how to perform the calculation. The starting point is the functional RG equation, 
which is also known as the Wetterich equation |237j . that is concisely written as 

= - M R kGk ) , ^to = ^ ( _^ feW +W, (143) 
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for generic scalar field theories. Here i, j indicate the meson species such as a and rr. The trace denoted 
by Tr involves the momentum integration as well as the sum over i, j. The scale k refers to the RG 
scale and the effective action is coarse-grained up to this IR scale k. The important step to reach 



Eq. ( 143 ) is to introduce a scale-dependent mass that satisfies the following boundary condition, 



Rk(q) -> 





k 2 



{q 2 > k 2 ) 
{q 2 < k 2 ) 



r 

s 



(k ->• 0) 
(k = A) 



(144) 



Then, the theory has /c-dependence through Rk{q), and one can define the generating functional Zk[J], 
and Wk [J] for connected diagrams in the standard method in the quantum field theory. The /c-dependent 
effective action, I 1 *. [</>], is a Legendre transform of Wfc[J]. We note that the regulator term should be 
properly subtract by — \ J d A q Rk(q)(f>(—q)4>(q) in the definition of Tk[4>] (see Ref. |238] for a pedagogical 
derivation). Since Rk{q) vanishes in the limit of k — > 0, the coarse-grained action amounts to the 
standard effective action V with full quantum effects, while I\.^a — > S where S is the classical action 



and the flow start at the UV scale A. Actually Eq. (143) is the exact formula to define quantum field 



theories containing equivalent information as the functional quantization method. Therefore, of course, 



it is an impossible task to integrate Eq. ( 143 ) unless the theory is solvable 



A frequently-used assumption to go into the practical calculations is that the non-locality in the 
effective action is neglected so that r^[0] can be decomposed into the tree- level kinetic term and the 



local potential. Corresponding to the Lagrangian density of the QM model in Eq. (142) the Ansatz in 
this local potential approximation (LPA) is taken as 



r k [a,7T) = f rf 4 x|^[z 7 ^-^(^ + H5T-7r)]^+^(^a) 2 + i(^7r) 2 + ^[(T,7r]| 



(145) 



and then one can convert the RG equation for into that for f^. The RG equation turns out to be 
surprisingly simple under the choice of the so-called optimized regulator [239J, 



Rk(q) = (k 2 - q 2 )d{k 2 - q 2 



(146) 



which is a useful choice for the finite-T calculation, while the covariant form would be more appropriate 
for the T = case. One can immediately confirm that this simple choice satisfies the requirement (144). 
Thanks to the Heaviside theta function, one can carry the g-integration out analytically to find, 



~dk 



127T 2 



E w \2TJ 



— - cothf — ^ 



24 

En 



{l - n F (-E q -/iq) - n F (E q +fi q )} 



(147) 



where the quark part has a coefficient given by twice of the fermionic degrees of freedom, 2N c Nf = 12. 
Here n-p(E) is the Dirac- Fermi distribution function at T, and the energies in the denominator are 
Ei = \/k 2 + Mf with 

M 2 a = nl , Ml = — . (148) 



Ml 



2 2 

9 o 



a 



The a meson mass is naturally the curvature of the effective potential, and the pion mass is also the 
curvature with respect to the ir field, which can be expressed by the a-derivative in an assumption that 
flk is a function of a 2 + 7r 2 (except for the explicit symmetry breaking term ca). 



There are two representative strategies to integrate the RG equation (147) with respect to k. One 
is to expand fifc(cr) around the minimum at each k. The whole shape of O^cr) is not necessary, but the 
meaningful observable is the potential value right at the physical point. Another method is to discretize 
fifc(cr) on the cr-grid and solve hundreds of differential equations numerically. In this latter formulation 
one does not have to find a minimum of the potential at each k. 



We shall make some comments on the RG equation (143). In the traditional application of the RG 



analysis, the critical phenomena associated with the second-order phase transition have been of central 
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interest and the critical exponents have been calculated around the IR fixed point. The applicability 



of the formula (143) is, however, not restricted only around the second-order phase transition. As we 
mentioned, Eq. (143) can be an alternative of the quantization scheme, and the full effective potential 
can be retrieved whether or not the system is far from the critical point. For the purpose to examine 
the first-order phase transition, the grid method is powerful enough to carve the double-well potential 
shape. In practical calculation, one should stop the fc-integration at some small k [60]; otherwise the 
potential becomes flat due to the convexity. 

It is a straightforward extension of the QM model to introduce the coupling with the Polyakov loop 



in the same way as in Eq. ( 137), and this extended model is called the PQM model. In the PQM model 
the resulting phase diagram has turned out to be quite similar to what we have already discussed using 
the PNJL model. In particular chiral restoration in the high density region has a general tendency to 
favor the first-order phase transition. If the existence of the first-order boundary is a robust feature of 
the phase diagram, it ends up with a terminal that is to be identified as the second-order phase transition 
point. Then, around the critical region where fluctuations are enhanced, the (P)QM model should by 
far surpass the PNJL model. For the phase diagram and the critical properties around the QCD critical 
point in the PQM model, see Refs. |240[ 1241] 12421 1243] and references therein. The question is, then, 
whether the first-order phase transition really exists or not. 

4.3 Liquid-gas phase transition of quark matter 

For the first-order phase transition at high density there is no universality argument. This is a dynamical 
problem depending on the details of the microscopic theory. Nevertheless, most of chiral models (without 
repulsive vector interactions) predict the first-order phase transition in the high-density region, and there 
must be some intuitive explanation for this general tendency |244] 1245] . 

Let us start our heuristic argument with a simple setup; cold and dense quark matter in a quasi- 
particle description. This means that we assume a Fermi liquid of quark matter, which should be valid 
for bulk thermodynamic quantities as long as T is small enough and the Landau damping is a minor 
effect. Then, the thermodynamic potential coming from a hot and dense medium is expressed as a 
function of the quasi-particle mass M q as 



r (M q ) = - dfi' Pq (fi') - 4N C N{T J ln(l + e'^' T ) , (149) 



where p q (/x q ) is the quark number density; p q (/i q ) = 2N C N{ J -^Sg [n-p{uj p — // q ) — np(uj p + /i q )] and the 



dispersion relation is u p = ^p 2 + M 2 in the quasi-particle picture. This form ( 149 ) is common in any 
chiral models such as the (P)NJL and the (P)QM models. It is then the vacuum part that is severely 
contaminated by model uncertainties. 

In the spirit of the Ginzburg-Landau expansion, we can postulate the vacuum part as a polynomial 
form, 

fi (M q ) = a(M 2 - M q ) 2 - bM q - cM 3 . (150) 



In the parametrization (150) a is the curvature of the potential at the origin, b represents the effect of 
the explicit symmetry breaking due to current quark mass, and c takes care of the instanton-induced 
interaction for Nf = 3. Here, for simplicity, we take c = 0; we are interested in the reason why a 
first-order phase transition is dynamically favored at high density, and if c ^ the tendency toward the 
first-order transition would be just strengthened by the tree-level contribution. The potential minimum 
is located around M q ~ Mq, so that M is a parameter to be identified as the constituent quark mass 
in the vacuum. 

Although this setup is extremely simple or even primitive, we emphasize that the above Oo + ^matter 
can grasp the essential properties of all chiral quark models, except for the uncanceled logarithmic singu- 



larity in Eq. ( 149 ), which we will ignore here because it is dangerous only near the chiral limit. Moreover, 
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Figure 18: (Left) Potential shapes from Eqs. (149) and (150) at T = and fi q = 370 MeV. fi mat ter is 
fairly model independent, while uncertainty is unavoidable in Qq. The parameters in Qq are chosen as 
Mq = 340 MeV, a = 0.05 ~ 0.08, b = c = 0. (Right) The saturation curve of the energy per baryon as 
a function of the baryon density. From the bottom to the top, the solid curves represent the results in 
the chiral limit without the vector interaction, with massive quarks without the vector interaction, and 
with massive quarks with the vector interaction. The dotted parts are unstable and the dashed curves 
represent the branches corresponding to the chiral spiral. Figures adapted from Ref. [245J. 



in the regime at high T, meson fluctuations may give rise to T-dependent coefficients. Therefore, the 
analysis in this subsection is valid only for /i q 3> T. 



The left panel of Fig. 18 shows the typical behavior of the effective potential. The matter part 
^matter always has a minimum at M q = because the baryon density is the largest when quasi-particles 
are massless. Then, together with Qq that has a minimum at M q = Mo, it is conceivable to expect two 
minima in the total Q, and the absolute minimum jumps from one to the other at the first-order phase 
transition. This is the mechanism why the density effect tends to favor the first-order phase transition. 
At the same time, this argument can explain why the location and even the existence of the first-order 
boundary on the phase diagram are such model-dependent. The largest uncertainty comes from a; in 
Fig. 18 uncertainty associated with a = 0.05 ~ 0.08 (empirical values in the NJL and the QM models) 



is shown. 

In the same way as the left panel of Fig. [7] for nuclear matter, it is useful to make a plot for the 
energy per baryon, s/pb at T = 0, which is presented in the right panel of Fig. 18 As we already 
discussed in Sec. |3.1 the minimum of this curve tells us the saturation density and the binding energy. 
Because we are now working for quark matter, the saturation density can differ from po, and moreover, 
there may not be a minimum at all. 

The existence of the first-order boundary of chiral restoration corresponds to the appearance of the 



minimum in the saturation curve. Thus, if quark droplets are possible [HS], as is the case in Fig. 
for massless 6 = 0, the first-order phase transition and the QCD critical point can be concluded. From 
these analyses it is evident that the first-order phase transition in quark matter has exactly the same 
origin as the liquid-gas phase transition of nuclear matter. Even when droplets are only meta-stable, as 
is the case in Fig. 18 for b = 0.08Mq (massive, without vector), the first-order phase transition occurs. 



We have not mentioned on the impact of the vector interaction as introduced in the Lagrangian 



density 124 The role of this type of the interaction is to add, 



fl vec {M q 



9vp\ , 



(151) 



to the thermodynamic potential at the mean-field level. It should be noted that the mean-field cal- 
culation needs not only the tree-level 
in Sec. 



3.1.1 



together wit 



-g v Pq term but also the shifted chemical potential //< 



as seen 



If expanded by the shift g v p q , this latter effect amounts to 1g v p% and Eq. (151) results 



i the tree-level term. 
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It is obvious from the saturation curve that the repulsive (g v > 0) vector interaction always disfavors 
the first-order phase transition. The energy e is pushed up by oc p| from the vector interaction, and 



eventually the local minimum is completely washed out with increasing g v , as exemplified in Fig. 18 for 
b = 0.08Mq and g v = 0.12/Mq (massive, with vector). These parameters are arbitrarily chosen here for 
demonstration. Once the absence of the local minimum becomes the situation with sufficiently large 
g v , there no longer appears a first-order phase boundary, and the QCD phase diagram has only smooth 
crossover entirely. 

Above mentioned is the simplest and the most intuitive explanation about the impact of the vector 
interaction on the order of the chiral phase transition (at T = 0) that was first addressed in a model 
study [246] and discussed repeatedly |225l 12471 12481 1249]. See also Ref. [?] for an attempt to estimate 
the strength of the vector interaction using the PNJL model. 

4.4 Inhomogeneous chiral condensates 

The existence of the QCD critical point is easily affected by unknown factors, while the inhomogeneous 
chiral condensate is more robust and, to the best of our knowledge, no counter-example against the 
inhomogeneous state has been found. Since the discovery of the relation between the QCD critical 
point and the Lifshitz point |250] , more and more attention is being attracted to inhomogeneous quark 
matter, though there were earlier theoretical studies [2511 1252] 1253] and also the crystalline color- 
superconductors were intensely investigated in the higher density region (see Ref. [19] and references 
therein) . 

The simplest way to introduce inhomogeneity in the chiral condensate is to adopt the Ansatz of the 
chiral spiral, and actually, the chiral spiral is very useful to deepen the understanding of the driving 
force for spatial modulation. Under the 1-dimensional chiral-spiral Ansatz |254j . the scalar and the 
pseudo-scalar condensates form a spiral (along the z-direction in our choice) as 

(ipip) = xcos(2qz) , (^7573^) = Xsin(2gz) (152) 
with a wave-number q. This is nothing but the quark-matter analogue of the p-wave pion conden- 



sation (117) in nuclear matter. Also, this Ansatz is sometimes referred to as the dual chiral-density 



wave 255 



The technical advantage of the chiral spiral is that this inhomogeneity is readily described by the 
redefinition of the basis; ip = e lj5T3qz ip' with a homogeneous condensate \ = (ip'tp'). Such a rotation in 
the chiral limit can remove the /i q -term in the Dirac operator, and then, the quasi-particle dispersion 
relation in the ^'-basis is expressed as 



with dynamical mass M q on the ip' -basis. This type of simple inhomogeneity pattern has been considered 
repeatedly in various contexts such as large- N c QCD [2561 1251] . the Overhauser instability |252] [255j. 



the quarkyonic chiral spiral |257j . and so on. The dispersion relation (153) should be plugged into 



^matter in Eq. (149) to evaluate the thermodynamic potential. 

The most important observation is that a large part of the mass effect can be absorbed if q ~ M q . 
Then, the density p q is not suppressed even with large M q . In fact, the (l+l)-dimensional system is 
this extreme example; p q has an origin in the quantum anomaly and is completely insensitive to M q 
but dependent only on the chemical potential. 



This simple observation from Eq. ( 153 ) can explain how the chiral spiral with q ~ M q can lower the 
total thermodynamic potential. If one makes a plot for n matter (M q , q), one can find that this mechanism 
works efficiently and ^ ma tter(^q, q) goes arbitrarily smaller for q ~ M q — > oo. This is a very simple and 
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Figure 19: Typical phase diagrams with the chiral spiral for two cases. The solid curves represent a 
first-order phase transition. The shaded region in the left shows the chiral spiral for the massless (b = 0) 
case without the vector interaction (g v = 0), and that in the right for the massive (b = 0.08Mq) case 
with the vector interaction (g v = 0.12/Mq) with M being a typical scale in the potential (150). Figure 
adapted from Ref. |245j . 



robust mechanism for the energy gain from finite q. This large energy gain should overcome the energy 
loss from the kinetic term in f2o(M q , q), which can be expanded as 

Q (M q , q) = Q (M q , q = 0) + («M q 2 + (3b)q 2 . (154) 

The first term with a > represents the kinetic term. One could estimate a using a chiral model, but 
one should be careful not to pick unphysical terms up from gauge-variant regularization. The term oc (3 
is also necessary to take account of the effect of the current quark mass that would disfavor the chiral 
spiral. Quantitative details depend on a and (3, but this setup with Q (M q ,q) + fi matter (M q , i?) suffices 
to capture the model-independent essence regardless of the choice of a and (3. 

Figure [T~9] shows typical phase diagrams resulting from this simple description with a choice, a = 0.25 
and (3 = 0.25/Mo. The chiral spiral is favored in the shaded region surrounded by the phase transitions 
of first order or second order. 

We would stress that the vector interaction chosen here (g v = 0.12/Mq) is large enough to make 
the QCD critical point disappear from the phase diagram. Yet, the chiral spiral persists for any g v as 



confirmed in Fig. 19 even for the massive case with strong vector interaction. Therefore, the possibility 



of inhomogeneous states in quark matter is a more robust prediction from the model than the QCD 



critical point. The corresponding saturation curves are overlaid on the right panel of Fig. 18 with dashed 
lines. As noticed in this figure, it is possible that a local minimum is revived, so that a first-order phase 
boundary can come back with inclusion of the chiral spiral structure. 

For more realistic ground states, it is desired to optimize the spatial structure beyond the Ansatz of 
the chiral spiral. This has been done by minimizing the Ginzburg-Landau potential and the solitonic 
solution is derived, which is known to have a lower energy than the chiral spiral in concrete model 
calculations |2581l259j . which is also systematically investigated using the Ginzburg-Landau theory [260J. 
The saturation curves associated with such solitonic solutions are discussed in Ref. [261J. Even though 
the chiral spiral is an oversimplified form as compared to the soliton, it is still useful to exemplify that 
the inhomogeneous state takes over the homogeneous one. More theoretical studies are needed to clarify 
the stability of inhomogeneous states with other interactions (analogous to the effect of g' in the p-wave 
pion condensation) or with meson fluctuations in the RG method. 
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Figure 20: (Left) Schematic phase diagram of large-iV c QCD. The pressure is 0(N®) in the hadronic 
phase, 0(N%) in the deconfined phase, and O(N^) in quarkyonic matter. (Right) Baryon-baryon 
interaction of 0(N C ) in the large-iV c limit. The exchanged quarks should have the same color. If they 
have the different color, one gluon exchange is necessary. 



4.5 Quarkyonic matter 

In the large- N c limit the deconfinement phase transition is well-defined and the phase structure asso- 
ciated with deconfinement is quite simple. As we have seen in Sec. [375 quark loops are suppressed by 



1/N C as compared to gluons, and yet, matter with finite baryon density may appear if the chemical 
potential is large enough. What does nuclear matter look like in the large-iV c limit? We already have 



sufficient ingredients to tackle this question - we have discussed a Skyrme crystal in Sec. 3.5.2 and 



discovered a large- N c phase diagram in Fig. 14 We have also clarified the order parameter behavior in 



Fig. 17 obtained in the PNJL model in which the large- N c limit is implicitly assumed 



4.5.1 Phase diagram and the pressure 

Apart from the fate of chiral symmetry, the phase diagram of large- iV c QCD is as simple as in the left 



panel of Fig. [20j which is equivalent to Fig. 14 without a curve representing T x . The horizontal line at 
T = Td is the deconfinement phase boundary below which only glueballs exist, the pressure of which 
should be of 0(N®). Because there are (iV c 2 — 1) gluons in the deconfined phase, the pressure sharply 
jumps from 0(N®) in the glueball (or hadronic) phase to 0(N%) in the deconfined phase. An onset for 
a finite quark density is located at /i q = m q and, if finite-density matter behaves as a free quark gas, 
its pressure is ~ N c fi^. Hence, the deconfinement phase boundary cannot be affected by quarks unless 

iVc/iq becomes comparable to N%, namely, // q ~ O(Nc^). Eventually, for fi q ~ 0(N^ 2 ), quarks are no 
longer suppressed by 1/N C which is compensated by /i q from the quark loop, and the gluon interactions 
are screened by dynamical quarks. 

The above argument suggests that the pressure in the right-bottom part is of O(N^), which is 
indeed the case if the state of matter is a (nearly) free quark gas and /x q ~ 0(N®). Such an argument 
is, however, too naive; as is clear from Fig. [14] in the Sakai-Sugimoto model, the ground state in the low- 
T and high-/i q region is identified as nuclear matter rather than quark matter. This is also consistent 



with model studies as in Fig. 17 which implies quark confinement with small <£> in this region. In fact, 
gluons are confined below T d , and it would be quite reasonable to assume that quarks are also confined 
there. 

In Ref. [21J it was pointed out that the pressure of large-iV c nuclear matter is of O(N^) and it 
resembles a pressure of quark matter. In nuclear matter at large iV c nucleons are infinitely heavy 
and static, so that their kinetic energy is suppressed by 1/Mn ~ ^(^c _1 ) an d the dominant pressure 
contribution comes from the NN interaction. Such a situation is correctly incorporated as a Skyrme 
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crystal, which should be an approximate description of quarkyonic matter. Baryons at large N c interact 
strongly and the NN interaction energy should be of 0(N C ), which is understood immediately from 



a diagram in the right panel of Fig. 20 the interaction is of order of the combinatorial factor to pick 
exchanged quarks up among N c quarks inside of baryon. One can make sure that not only the two- 
body but also the multi-baryon interactions generally scale as 0(N C ). In summary, in the large- N c limit, 
nuclear matter is a system of strongly interacting baryons, and it looks like knowing quark degrees of 
freedom in it. In Ref. [21] a name was given to such a baryonic and quark-like state of matter, i.e. 
quarkyonic matter. 

If quarkyonic matter is a new state of matter, there must be some order parameter to characterize it. 
This is a natural question but is a source of confusion about the physical interpretation of quarkyonic 
matter. An intuitive picture of quarkyonic matter is the following; quarks deep inside the Fermi sphere 
are weakly interacting, because it is hard to excite these quarks above the Fermi sea due to Pauli 
blocking. On the other hand, the quarks near the Fermi surface with a shell-width ~ Aqcd are not 
affected so much from the Pauli blocking and can interact strongly. Thus, the bulk thermodynamics 
such as the pressure, entropy and so on are dominated by the quarks inside of the Fermi sphere, while 
the physical excitations on top of the Fermi surface are dominated by color-singlet baryons. 

4.5.2 Characterization of quarkyonic matter 

A frequently-asked question is; what is the definition of quarkyonic matter? In other words, what is 
the order parameter for quarkyonic matter? 

Formal definition In accord with the original argument in Ref. [2T], the definition of quarkyonic 
matter could be a state of matter that satisfies the McLerran-Pisarski conjecture. Here we introduced 
a term - the McLerran-Pisarski conjecture - which claims that a system of dense baryons be dual to a 
system of quarks. Thus, in quarkyonic matter, both descriptions of baryonic matter and quark matter 
can work to capture correct thermodynamic properties. 



In Sec. 3.5.3 we mentioned on two density sources; the D4 instantons corresponding to baryons 
and the string sources corresponding to quarks in the Sakai-Sugimoto model. Although the baryonic 
state has a lower energy in the analysis in Ref. [22], such a treatment of the density sources may have 
some subtlety (for example the validity of the DBI action with a cusp singularity, see Ref. [263J ) . If 
the baryonic and quark sources turn out to be indistinguishable after all, it would be the clearest 
holographic example of quarkyonic matter. 

Phenomenological signature This above formal definition is not radical but quite understandable 
in view of the success of the thermal statistical model or the hadron resonance gas (HRG) model to 
capture the QCD thermodynamics [264] . (For details about the model setup, see Refs. |265[ [266J and 
also Ref. |267] for a textbook.) In the HRG model it is assumed that hot and dense matter consists of 
non-interacting mesons, baryons, and all resonating states. This is a modern renovation of Hagedorn's 



old idea (see Sec. 4.1 ). Suppose that all observed states obey the Hagedorn spectrum (122), the partition 
function blows up at T = Th, and the growing behavior of the pressure near T c fits well with the one 
known from the lattice-QCD simulation as in Fig. |4j This observation is quite consistent with the 
behavior of QCD at large N c . As we discussed in Sec. |3.5[ the system is reduced to infinite towers of 
non-interacting mesons in the large- iV c limit of QCD. 

It is also known from the chiral effective models that the thermodynamic properties near T c are well 
understood in terms of quasi-quarks and gluons [224]. Therefore, thanks to a smooth crossover from 
the hadronic phase to the quark-gluon plasma, one can view the state of matter in this transitional 
region as both hadronic and quark-like. In a broad sense, the crossover region of QCD at zero density 
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Figure 21: (Left) Phase diagram with a region that looks like a triple point where the hadronic phase, the 
deconfmed phase (quark-gluon plasma), and quarkyonic matter meet. Figure adapted from Ref. [269J. 
(Right) Phase diagram with interweaving chiral spirals. Figure adapted from Ref. [270J. 



satisfies the McLerran-Pisarski conjecture, and the state of matter is quark-mesonic instead of quark- 
baryonic = quarkyonic . 

In a narrow sense, quarkyonic matter should be distinguished from quark-mesonic matter. Baryons 
interact strongly, and this situation is so different from the description by the HRG model. In the 
thermal statistical model the interaction effects are indirectly taken into account via the Van der Waals 
force with a parameter of the excluded volume [268J. However, in physical processes, the interactions 
should encompass the quark degrees of freedom as depicted in the diagram of Fig. 20 and the excluded 
volume factor is not adequate to go into the regime of quarkyonic matter. Thus, in the beam-energy 
scan of the relativistic heavy-ion collision, if any deviation between the experimental data and the HRG 
model fit becomes more and more appreciable, it would hint quarkyonic matter. 



Chiral property Nuclear physicists would still wonder how quarkyonic matter can be different from 
ordinary nuclear matter in the real world with iV c = 3. There is no clear-cut distinction. In nuclear 
matter kp is as small as the pion mass and the system is dilute (in a sense that one can perform the 
expansion with respect to kp), while quarkyonic matter is dense by definition. 

One possible characterization of quarkyonic matter is how it breaks chiral symmetry. In quarkyonic 
matter with infinitely large N c inhomogeneity is expected |251] 1257] . This spatial modulation is at- 
tributed to the pseudo 1-dimensional nature on top of the Fermi surface at high density. Interestingly, 
as briefly noted in Sec. 3.5.3 the Sakai-Sugimoto model exhibits instabilities with respect to spatial 
modulations [197|, 1198] . However, this sort of inhomogeneity is nothing peculiar for nuclear physics; as 
addressed in Sec. 3.6, the p-wave pion condensation is essentially the same as the (quarkyonic) chiral 
spiral. Then, we have to conclude that quarkyonic matter is not really distinct from nuclear matter, 
or the p-wave pion condensation may be smoothly connected to the quarkyonic chiral spiral and vice 
versa. 



4.5.3 Implications to the phase diagram 

With quarkyonic matter in addition to the hadronic phase and the quark-gluon plasma, the phase 
diagram should take a structure of the left panel of Fig. |2TJ and then there may appear a region that 
looks like an approximate triple-point as speculated in Ref. |269j (see also the lattice-QCD results in 
Ref. [2TT]). 

Later on, in Ref. [270J, this phase diagram was brushed up into the form of the right panel of 



Fig. 21 Here the essential alteration is to add the phase structures associated with the quarkyonic 
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chiral spirals. The chiral spiral is favored due to the pseudo 1-dimensionality, and the Fermi surface 
should be covered with such pseudo 1-dimensional patches. Because the transverse size of each patch is 
expected to be ~ Aqcd, the number of patches should increase with increasing /i q (and increasing area 
of the Fermi surface). This "interweaving" picture naturally leads us to the speculation of successive 
phase transitions associated with the number of patches. In other words, the most optimal crystalline 
structure should depend on the density and the structural changes would divide the inhomogeneous 
region (as indicated by the shaded area in Fig. 19 ) into finer substates named the interweaving chiral 
spirals [270] . 

The most serious concern on reality of inhomogeneous quarkyonic matter is the validity of the 
extrapolation from N c = oo to N c = 3, which is discussed for example in Ref. |272j . We finally point out 
that there exist quite a few theoretical arguments and pictures that are essentially indistinguishable from 
the McLerran-Pisarski conjecture and quarkyonic matter at N c = 3. A related duality (at zero density) 
has already been discussed in Ref. |273] . At high density the quark-hadron continuity in the context 
of color-superconductivity |274[ 1275] is not really different from the McLerran-Pisarski conjecture (for 
the confinement properties in color-superconductor, see Ref. [13] )• At the phenomenological level the 
percolation model may provide us with an interpretation of quarkyonic matter at N c = 3 |276l 1277] . 



5 Experimental Prospects 

There are many theoretical speculations and conjectures on QCD at high density, and some of them 
should be filtered out by experimental data. Here we concisely discuss experimental implications to 
explore the properties of dense QCD matter. 

5.1 Liquid-gas phase transition and the critical point 
5.1.1 Experimental signals in nuclear matter 

For general discussions on the liquid-gas phase transition, see Ref. [92] and references therein. Here 
we enumerate experimental signals for the critical point in nuclear matter. The extensive study of 
nuclear multi-fragmentation has been stimulated by the suggestion that this process might be a critical 
phenomenon |278] 12791 280, 281J. In particular the produced fragment distributions follow a power law 
of the form P(A) ~ A~ T with the atomic number A, which suggests that a class of "universality" exists 
behind it [282]. Available multi-fragmentation data in nucleus-nucleus and hadron-nucleus reactions 
provide a distinct change of its character. Given these signals, the observation of the liquid-gas phase 
transition in a finite nuclear system is strongly supported |283] 12841 12851 12861 1287] . 

The caloric curve, i.e. the heat versus temperature, is one of the important observables to distinguish 
the critical phenomena. Temperature of the system has to be reconstructed from observable quantities. 
In chemical and thermal equilibrium the temperature is extracted from the yields of isotope ratios [288J . 
The observed isotope temperature as a function of the excitation energy per nucleon E/A shows a 
plateau over the range of relatively low E/A [289] I29U] . The observed caloric curve agrees qualitatively 
with predictions of the multi-fragmentation model [280] which is an event generator as a generalized 
liquid-drop model for hot nuclei, and the properties are consistent with a liquid-gas phase transition. 

Another experimental signal of the phase transition is the appearance of a negative specific heat 
capacity: The entropy has typically a convex structure when a system experiences a first-order phase 
transition and this results in the specific heat C being negative. This effect has been extracted from 
analysis of the energy fluctuations and the presence of a negative C branch was indeed found |291j . 

The location of the critical point T c for the nuclear liquid-gas phase transition comes out from the 
observables. At this point the isotherm in the phase diagram has an inflection point. The major source 
of experimental information for T c is the fragment yield. However, the procedures to extract T c are 
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highly scheme-dependent and the obtained values lie in a wide range T c ~ 5-20 MeV, indicating a severe 
model-dependence to be resolved (see e.g. Ref. |292] ). The critical exponents have also been extracted 
from the moments of the fragmentation charge distributions [293L 12941 1295] . The analysis concludes 
that the phase transition is of second order and finds the critical exponents consistent with the Z 2 
universality class. 



5.1.2 Toward the QCD critical point 

First-order phase transitions for cold and dense quark matter have been predicted in several approaches 
using chiral models [2251 1223 E3], Dyson-Schwinger equations [22E1 12221 EDO], and lattice-QCD in the 
strong coupling limit |3U1[ I3U21 1303] 1304] !3U5j . Given the observation of a crossover at zero chemical 
potential from lattice QCD computations, this might suggest an additional critical point other than 
that of the nuclear liquid-gas in the QCD phase diagram. 

The order of the QCD phase transition at low temperature and high density is not established yet 
and thus the existence of a QCD critical point by itself remains an issue under debate. If it exists, 
the critical point must exhibit the critical exponents of the 3-dimensional Ising model that belongs to 
the Z2 universality class |306] . A potential signal for the QCD phase transition is modifications in the 
magnitude of fluctuations or the corresponding susceptibilities. In particular, fluctuations related to 
conserved quantities (baryon number B, quark number q = 3B, electric charge Q, strangeness S, energy 
E, etc) play an important role since they are directly accessible in experiments [307J. 

A phase transition can be probed with response of a medium to temperature and chemical potential. 
The n-th order cumulant of a conserved quantity X, such as B, q, Q, S, etc, is computed from the 
partition function Z via 

x 1 d n In Z 
Xn ~ VT* d(fi x /T) n ' (155) 

By taking the following cumulant ratios a volume factor can be eliminated; 

Rim = ^ " ( 156 ) 

Am 

The ratios R\ 2 and R® 2 are shown to be sensitive to quark deconfinement at /i q = and they exhibit 
qualitatively different behavior from those of the HRG model |308j . At finite /i q higher order cumulants 
of baryon number are particularly important since they diverge stronger on the chiral phase boundary 
in the chiral limit as well as at the critical point [309J. Higher moments and their signs have also been 
proposed as sensitive probes to the phase transitions [310J. 

Recently it has been shown that the fourth-order cumulant and the ratio R^2 (kurtosis) become 
universally negative when the critical point is approached from the crossover side [31 lj. Here the 
variance ax, the skewness Sx, and the kurtosis k x are defined as 

a 2 x = VT 3 xZ, S x a x = ^, n x a x = 4 • (157) 

A2 A2 

A study using scaling functions shows that the fourth-order cumulant Xa changes its sign around the 
crossover region from positive at lower temperature to negative at higher temperature, as a generic 
feature [312] . The temperature at which xf becomes negative is not universal since it depends on the 
regular part of the partition function. This suggests that the negative kurtosis is not a unique indicator 
of the critical point. 

The STAR collaboration reported the ratio Rf 2 |313j . which is consistent with the HRG model 
result |314] 1315] . More recent data of the moment products, Kxa x and Sx&x, of the net proton 
number and electric charge distributions |316] 1317] 1318] differ from the HRG model result. This might 
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be due to such simplified modeling of the QCD medium, although the thermal statistical model or the 
HRG model works well with various particle number ratios. So far, no experimental indication of the 
critical point is found. 

If any signal of a first-order phase transition were observed, this would become an indication of the 
presence of the critical point as in the case of nuclear matter. In this respect, it is also challenging 
in heavy-ion collisions to observe the proposed consequences of the spinodal decomposition associated 
with a would-be first-order phase transition at high density, such as an enhancement of baryon and 
strangeness fluctuations |319l 13201 1321] and a negative quark number susceptibility [322J . 



5.2 Dilepton measurements and chiral symmetry restoration 

Dileptons are considered to be promising probes to study for changes of hadron properties in matter 
since they pass through the fireball created in heavy-ion collisions without hadronic interactions. The 
short-lived vector mesons like the p mesons are expected to decay into dileptons inside the hot and 
dense matter. An enhancement in the dilepton rates below the p/u resonance - indicating the medium 
modification of the vector mesons - has been observed in many experiments. Although a major number 
of experiments detected a strong width broadening and no "mass shift" , no perfect agreement among 
different experiments is reached yet |109[ 13231 1324j . A shift of the p/u peak in the in-medium spectral 
function is often regarded as a "mass shift" of the vector meson. However, this is clearly an oversimplified 
interpretation since complicated reaction processes, numbers of resonances in a medium, and their 
interactions are involved and those complications are integrated into the final spectra. Therefore, 
picking up a maximum of such broad distributions does not readily give us the in-medium mass of the 
vector meson. 

Below we will briefly argue some theoretical issues related to the dilepton measurements. 



5.2.1 Is vector dominance fulfilled in a medium? 

A lepton pair is emitted from the hot and dense matter through a decaying virtual photon. The 
differential production rate in the medium for fixed temperature T or baryon density pb is calculated 
with the imaginary part of the photon self-energy Imll via 

^{q ,q;T,p B ) = -g^ -— [ _ - Imn(g , g; T, p B ) , (158) 

where a = e 2 /47r is the electromagnetic coupling constant, M is the invariant mass of the produced 
dilepton, and q p = (go, q) denotes the momentum of the virtual photon. We will focus on an energy 
region around the p-meson mass scale in the following argument. In this energy range the photon 
self-energy is expected to be dominated by the two-pion process and its imaginary part is related to the 
pion electromagnetic form factor J-"(s; T) through 

ImII( S ; T,p B ) = —i^-^^j \F(s;T,p B )\ 2 (159) 

with the pion mass m n and the form factor given by 

TV T \ / rp \ | 9p(T, Pb) ■ g P im(T, p B ) ftari\ 

m 2 p (T, p B ) - s - 6(s - 4m£) im p {T, p B W P {s; T, p B ) 

In matter-free space, a single photon mainly couples to two pions via the vector meson exchange, 
known as vector meson dominance (VD). This is assumed to hold in a hot and dense medium in many 
model calculations for the in-medium vector spectrum. One should, however, keep in mind that this is 
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not a priori justifiable. A reliable theoretical framework to handle the VD is a hidden local symmetric 



(HLS) approach as discussed in Sec. 3.4 the vector meson properties are essentially controlled by one 



parameter a and a particular choice, a = 2, leads to vanishing photon-pion coupling, 

9^n{T = P b = 0) = . (161) 

A consistent achievement of chiral symmetry restoration at a given set of T and pe requires that the 
parameter a should evolve toward unity with T and p-Q |325[ 13261 1327] . The immediate consequence is a 
strong violation of the VD in matter: The coupling g 77r7r is an increasing function of thermal parameters 
whereas the vector-meson-photon coupling is a decreasing function of them toward chiral symmetry 
restoration as 

1 . . 

fiW -> ^ > 9 P -+ • ( 162 ) 

This indicates more dilepton rates from two-pion annihilation not involving the in-medium vector mesons 
when the system approaches the chiral phase transition. Thus the dilepton yields become suppressed 
there [328] . 

5.2.2 Has the BR scaling been excluded by the NA60 dimuon data? 

The NA60 collaboration reported a strong in-medium broadening of the p meson and no mass shift in 
dimuon measurements [329, 330J. The data were compared with the theoretical predictions based on 
phenomenological Lagrangians (width broadening due to hadronic many-body effects versus dropping 
p mass), and the broadening with no mass shift is apparently favored. The measurement has excluded 
the scenario with a naive dropping mass as in the form, 

m p (n) ~ mJn = 0) ( 1 — a — ) , (163) 

where n represents either temperature or density, a a constant, and n an appropriate scale for normal- 
ization. 

If partial restoration of chiral symmetry modifies the vector meson property, it would appear not 
only in the mass but also all the interactions involving the vector meson. In principle those modifications 
should be incorporated systematically within a given effective theory for relevant hadrons. Therefore, 



a naive replacement of the vector meson mass with Eq. (163) in the quantities such as the spectral 
function is incomplete. Also, the VD is shown to be violated at a certain n once the vector meson, 
if its mass drops, feels partial chiral symmetry restoration [325, 1326[ 1327] . On top of the in-medium 
broadening and the suppression of the vector-meson-photon coupling due to the VD violation, it was 
argued |3311 1332[ 1167] that the signal of the Brown-Rho (BR) scaling is totally hidden in matter-free 
contributions essentially coming from the vector mesons which do not feel partial restoration and thus 
peaked around the vacuum mass, m p = 770 MeV. Thus, a comparison of the data with the "dropping 
mass scenario" needs to be made with great caution. 

Aiming to characterize a broad spectrum without clear resonance structure, an "average mass" can 
be introduced as |333j [334J 

- 2 _jdssp VA 
J dsp V:A 

At order T 2 the vector meson pole mass does not receive any thermal corrections |3351 1336j and the 
numerical analysis for fhy using QCD sum rules was shown to be consistent with this statement. The 
vector meson "mass" fhy stays almost constant up to T ~ m n , whereas the axial- vector meson "mass" 
ffiA decreases with rising temperature. Therefore, in this range of the temperature a dropping p mass 
is not supported |334] . 
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On the other hand, in-medium p mesons in cold nuclear matter exhibit a different feature: The 
calculated my making use of the in-medium spectral function follows 



^L^^^l_0.15^. (165) 
m™ U Po 

This suggests that the tendency of the BR scaling is indeed visible with the in-medium spectral function 
which shows a quite broad distribution. This observation is, contrary to the first impression, a good 
example of the "hidden" BR scaling in the broad spectrum in matter. 

5.2.3 Chiral mixing in vector and axial- vector spectral functions 

There is no doubt that dilepton measurements have observed in-medium modifications in the vector- 
meson channel. Yet, it remains unclear how such broadened vector spectra would eventually be linked 
to chiral symmetry restoration. Although the axial-vector spectral function is hardly measurable, axial- 
vector mesons come in to the vector spectrum via a mixing to pions in hot/nuclear matter [335, 337J. 
The mixing between the vector and axial-vector spectral functions, py and pa, is derived in a model- 
independent way at low T and pe, arranged in a low-temperature/density theorem, 

p v {s; n) = (1 - e)p v {s; 0) + ep A {s; 0) , 
p A (s; n) = (1 - e)p A (s; 0) + ep v (s; 0) , 

where chiral mixing is characterized by the parameter e = T 2 /6f% or Ap-BO~ 7T ^/3f 2 m 2 with the pion- 



nucleon sigma term a^. The spectral functions with the mixing Eq. (166) at finite temperature has 
been explored in the QCD sum rules |338[ I334[ I339[ 1340] where one finds an indication of the expected 
tendency of the spectra py and pA becoming degenerate. 

Properties of the vector and axial-vector mesons constrained by current algebra are summarized in 
the sum rules for the spectral functions, known as Weinberg sum rules |175j . 

ds r°° 

-[p v (s) - PA (s)] =fl, \ ds[ Pv (s) - p A (s)] = . (167) 
s Jo 

Those sum rules extended to a finite temperature have been explored |341] where the following relation 
was deduced assuming the pole-saturated forms of py and pa', 

f 2 (T) ( m 2 n ml \ 

JA ) - 2Z p (T) —f- - —j— (168) 



ft PV 7 WOO <(T) 

with the pole residue Z p . Clearly, when m p (T) = m ai (T), chiral symmetry is restored with the vanishing 
order parameter f n (T) — > 0. The sum rules by themselves do not distinguish how the masses behave 
with increasing temperature. Possible scenarios are (i) both decrease, (ii) both increase, (iii) m p (T) 
increases whereas m ai (T) decreases. Using an Nf = 2 linear sigma model, m p (T) is shown to go up 
by chiral restoration temperature if the VD holds in hot matter [342] . (See also Refs. [343L 1344] for 
the linear sigma model with global chiral symmetry, which can be extended to non-zero temperature.) 
Non-linear sigma model calculations also yield the same tendency at low temperature where the VD is 
well satisfied |327] . If the VD is violated, there is no unique prediction unless one imposes a certain set 
of conditions for the model parameters being consistent with restoration of chiral symmetry. 

The spectral functions py and pa become degenerate when the mixing parameter reaches e = 1/2. 



However, such a naive extrapolation of Eq. (166) to higher n requires the vector and axial- vector meson 



masses staying their vacuum values and thus they are not necessarily degenerate in contradiction to 



the desired property of chiral partners. An approach which does not rely on Eq. (166) is to use chiral 
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effective theories of the vector and axial-vector mesons explicitly. The mixing and its evolution in 
matter are generated from processes involving the tt pax interaction. Due to the mixing the axial- 
vector meson appears as a bump in py around y/s = m ai — m n and this makes py broadened. When 
the system evolves toward chiral symmetry restoration with temperature, it was shown that the chiral 
mixing vanishes at a critical point T c in the chiral limit |173j . Therefore the "maximal mixing" scenario, 
e = 1/2, does not happen but instead, the resonance peak and the bump tend to be degenerate forming 
a single peaked structure which obviously indicating the degenerate m p and m ai . It is also interesting 
that even with explicit symmetry breaking the p and a\ masses are predicted to be well degenerate 
around a pseudo-critical temperature. 

At finite chemical potential chiral mixing is generated at tree level since charge conjugation invariance 
is lost, which allows the following term that mixes vector and axial- vector fields [345J, 

£ mix = 2Ce 0vXfT tx[d v V x ■ A a + d v A x ■ V a ] (169) 

with the coupling C and the total anti-symmetric tensor e 0123 = 1. One can deduce the above term 
from the u pa\ part in the Wess-Zumino-Witten term with replacing with its expectation value, 

(W ) = QujN N ~W ■ (170) 

ml 

Therefore, the coupling C is 

C = g wpai {uo) , (171) 

and at normal nuclear density po one finds C ~ 0.1 GeV. The mixing term is also derived from the 
reduction of five-dimensional Chern-Simons action to four dimensions in a holographic QCD model [346J 
where C ~ 1 GeV at pb = Po was found. Such a strong mixing is, however, unrealistic since otherwise 
it yields vector meson condensation slightly above po- Thus, a large C is supposed to be an artifact in 
large N c and/or of the specific model construction applied to a dense system. 



The mixing term (169) modifies the dispersion relation resulting in 

1 



2 -2 

Po ~ V =2 



m v + m\ ± J (m\ — m v ) 2 + 16C 2 p 2 



;i72) 



which describes the propagation of a mixture of the transverse p and a x mesons. The longitudinal 
polarizations follow the standard dispersion relation. Consequently, the spin-averaged current correlator 
is superposition of the longitudinal mode peaked at my and the transverse modes with two bumps 



peaked at shifted my and m^ obeying Eq. (172). The presence of mixing makes the entire spectral 
function broadened and this influences over dilepton production rates. Its impact on the rates crucially 
depends on the size of C. Let us suppose that the lowest u dominance giving rise to C ~ 0.1 GeV 
at pb = Po is favored. Mixing does not have significance at this density. With increasing density, i.e. 
increasing C, the structure of the spectral function starts to show a change from the vacuum distribution, 
and such a modification becomes distinct when density reaches ~ 3po [345J. 



The mixing (169) is chirally symmetric and thus does not vanish at chiral restoration in contrast 
to the vanishing chiral mixing at finite temperature without baryon density. In the presence of both 
mixing effects, a tendency of chiral symmetry restoration signaled by degenerating p and a\ mesons 
would become more obscure. Therefore, it seems a hard issue to observe a clear signal of chiral symmetry 
restoration in dilepton measurements. At least from the theoretical side, the in-medium vector spectral 
functions known to agree with data can be compared with the axial-vector spectra calculated in a given 
model with those mixing. Unless the description for the hadrons involved in terms of local fields breaks 
down totally, such a comparison could offer us an indirect evidence for the symmetry restoration, with 
attention to the notion of in- medium mass. 
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5.3 Astrophysical Implication 



Historically speaking, the research on dense nuclear matter and the microscopic calculation of the EoS 
were motivated in the context of the structures of compact stellar objects (see e.g. Ref. |347j for a 
textbook). 

There was an astonishing news in the year of 2010; it was reported in Ref. [25] that the Shapiro 
delay measurement for the binary pulsar PSR J1614-2230 results in the pulsar mass of (1.97±0.04)M Q . 
This discovery has a tremendous impact on nuclear physics: Indeed, the abstract of Ref. [25] reads, 
11 effectively rules out the presence of hyperons, bosons, or free quarks at densities comparable to the 
nuclear saturation density." Once the EoS is known, one can plug p(p), the pressure as a function of 
density, into the Tolman-Oppenheimer-Volkoff (TOV) equation to solve the relation between the mass 
and the radius of compact stars. The mass cannot be arbitrarily large but bounded by an upper limit 
corresponding to p(p) that should sustain the star against the gravitational collapse. The EoS should 
be harder, i.e. p should be larger for a given p, to support a larger mass. 

The problem is that the EoS tends to be softer with exotic compositions, with which it may not be 
consistent with the existence of the neutron star as heavy as almost twice of the solar mass, though 
it can be; see Refs. |348[ 1349] . There is not really a final statement on this issue, and actually it is 
quite difficult to establish any reliable derivation of the EoS of quark matter (see e.g. Ref. [350] for a 
recent attempt), and thus it can become harder than believed (see Ref. [351] for example). Also, if the 
repulsive vector-interaction as in Eq. (124) is strong enough, it would easily make the EoS sufficiently 
hard. In other words, to explain the neutron star with ~ 2M , a substantially large vector interaction 
should be expected, and this implies in turn that the QCD critical point is disfavored from the QCD 
phase diagram [248J. 

Without assuming the state of matter, it is also a feasible program to compile astrophysical ob- 
servations to constrain the EoS of dense matter directly. For recent progresses in this direction, see 
Refs. [3521 13531 13541 1355] . 

The cooling rate is also informative on the state of matter in the neutron stars [356] . and it is an 
interesting question how to constrain the possibility of the pion condensation [357] . We must keep 
in mind, however, that the environment of the neutron star is not necessarily the same as that we 
consider for the QCD phase diagram or for the heavy-ion collision. In the neutron star /3 equilibrium is 
achieved, while it is usually symmetric nuclear and quark matter that is relevant to the phase diagram 
research. Besides, even if exotic contents such as the pion condensation, the kaon condensation, and 
quark matter were ruled out from the neutron-star observation (which is difficult from the mass-radius 
relation alone [348] ) . they might be still possible in symmetric matter and/or at higher densities. 



6 Summary and Outlook 

In this review we have addressed physical properties of nuclear and quark matter and given pedagogical 
descriptions on the tools used in theoretical research. On the phase diagram with two environmental 
parameters, the temperature T and the baryon chemical potential only a small portion has been 
understood; the QCD phase transitions of deconfinement and chiral restoration at (nearly) zero density, 
and the nuclear liquid-gas phase transition that is inevitable from the saturation property of nuclear 
matter. For the zero-density crossover, a pile of experimental and lattice-QCD data have been accu- 
mulated. Also the hadron resonance gas model works nicely, with which a physical picture has been 
established. In the nuclear physics territory the correct understanding is guided by experimental data, 
but the direct application of QCD to nuclear matter is still a big challenge. In practice one cannot 
avoid modeling the QCD dynamics in a form of the effective description. 

In nuclear matter not far from the normal nuclear density chiral symmetry cannot be totally restored, 
but partial restoration has been confirmed both in theory and in experiment of the deeply bound pionic 
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Figure 22: Two representative scenarios for the QCD phase diagram; one with the QCD critical point (s) 
adapted from Ref. [18] and the other without the first-order phase transition at all. "Quarkyonic Matter" 
is intentionally replaced with "Quarkyonic Regime" (see an explanation in the text). 



states. Mesons and baryons should be modified accordingly, and the hidden local symmetry model 
provides us with the most consistent description on these issues. It is, however, still hard to make 
conclusive statements on the changes of the hadron spectra and the mixing pattern at finite density. 
To grasp the qualitative character of dense matter, it is useful to increase the number of colors to 
infinity; N c — > oo. Then, the Skyrme model and also the most promising holographic QCD model 
(Sakai-Sugimoto model) tell us that the ground state of nuclear matter at N c — >■ oo takes a crystal 
structure and spatial modulations should appear at high density, which is reminiscent of the old idea 
of the p-wave pion condensation. 

Interestingly enough, this sort of inhomogeneity is favored by a robust mechanism in quark matter. 
Using a generic quasi-quark picture, we have demonstrated a simple calculation how the first-order 
phase transition and the inhomogeneous states can emerge in such a way induced by the density effects. 
This simple argument is helpful in understanding what part is robust and what part is not in theoretical 
predictions. We have demonstrated that the first-order phase transition and thus the QCD critical point 
can disappear from the phase diagram if the repulsive vector interaction is substantially strong, which is 
also suggested from the existence of the two-solar-mass neutron star. With the same parameter set, on 
the other hand, the inhomogeneous (chiral spiral) phase can persist, and we can say that the appearance 
of inhomogeneity is quite robust unlike the QCD critical point. 

Then, it is conceivable that such inhomogeneous states of quark matter may be connected to the 
p-wave pion condensation. More generally, it might be even possible that strongly interacting nuclear 
matter has a dual description in terms of quark degrees of freedom. In fact, such "duality" is suggested 
in the large- N c counting, which defines a new regime on the phase diagram called quarkyonic matter. It 
was once said that quarkyonic matter could be a state with confinement but chiral symmetry restora- 
tion, but this characterization was misleading. The crucial aspects of quarkyonic matter are; (1) strong 
interaction of baryons, (2) inhomogeneous chiral condensates, and (3) quark degrees of freedom. Model 
buidling with these properties correctly equipped along the similar line to the hadron resonance gas 
model is an urgent theory problem. 



What should the QCD phase diagram look like? In the market of the QCD phase diagram research, 



there are many variants of the schematic figures. We pick two representative scenarios up here in Fig. 22 



The left figure (taken from Ref. [18]) is similar to the conventional type of the QCD phase diagram. It 
has a first-order boundary for the chiral phase transition and accommodates the critical point (s). The 
right one is a rather unconventional but these two phase diagrams are both viable enough. The QCD 
phase transitions, namely, deconfinement and chiral restoration, can be only crossovers and no critical 
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point (apart from the nuclear liquid-gas transition) exists on the whole QCD phase diagram. Even in 
this case the inhomogeneous states around the regime of quarkyonic matter should be robust, as pictured 
with the dashed curve, in a sense that no model has falsified yet. In particular, in this latter scenario, the 
quarkyonic chiral spiral, the p-wave pion condensation in nuclear matter, and possibly the crystalline 
color-superconductivity may well be linked without sharp phase boundary, which is collectively referred 



to as "crystalline state" in the right panel of Fig. 22 In these figures we intentionally use a term 
"Quarkyonic Regime" instead of quarkyonic matter because it is important to note that quarkyonic 
matter is not a distinct state of matter. Quarkyonic matter should be understood as a transitional state 
between dense nuclear matter and quark matter, having both aspects of them. 

With available theoretical knowledge and experimental data, it is still far from feasible to be able to 
select out the genuine one among many proposed diagram. We must solve QCD to establish a unique 
picture of the QCD phase diagram, and at the same time, we need more experimental data of the 
beam-energy scan from the relativistic heavy-ion collision. The sign problem inherent in the Monte- 
Carlo simulation is extraordinarily difficult to handle in a satisfactory manner, and the QCD phase 
diagram research awaits some other technical breakthrough not relying on the importance sampling. 

As an extrapolation from nuclear matter toward higher density, one need to resort to effective 
approaches. Applying effective chiral Lagrangians in terms of a given set of mesons and baryons to 
dense matter is based on the assumption that those hadrons keep their particle identity and can still be 
expressed as local fields. It is not obvious whether this approximation is legitimate, but rather it can be 
taken as a working hypothesis. In fact, the spectroscopic factor for various nuclei indicating a deviation 
from the fully occupied mean-field orbits clearly shows that single-particle-ness of excitations reaches 
~ 70% [358]. This would encourage modeling dense matter in a quasi-particle picture. In particular, 
when partial restoration of chiral symmetry sets in, light degree (s) of freedom that will be degenerate 
with the NG bosons must appear in a system. Thus, if one chooses "appropriate" degrees of freedom, 
such effective Lagrangian approaches with the local-field approximation would give a good description. 

This speculative consideration is also supported by our knowledge on condensed matter physics 
where strongly correlated systems are well described by effective field theories and the notion of quasi- 
particles has been successful. A considerable example is a phase transition between the Neel (antiferro- 
magnet) and valence bond solid (VBS, paramagnet) phases. Topology of those phases is well captured 
by the non-linear sigma models and a skyrmion configuration splitting into two half-skyrmions has been 
known [359, 360] 1361] . The effective theory contains an emergent gauge field that mediates interactions 
between effective degrees of freedom [362] . This is quite suggestive and a similarity would be expected 
in nuclear many-body systems in high densities. Indeed, the concept of induced gauge fields, which arise 
when "fast" modes separated from "soft" modes are integrated out, have been recognized in various 
fields of physics |363j . In the context of QCD, non-Abelian gauge potentials (Berry phases) appear 
naturally in topological chiral bags under any adiabatic rotation 



On the QCD phase diagram we must widen our perspective outside of nuclear matter toward higher 
temperature and/or density, and the validity of effective approaches is questionable then. To this end, 
there is no bypass but solving QCD is the one-track way to reveal the QCD phase diagram completely. 
This used to be an impossible problem, but nowadays, it is becoming tractable owing to the technical 
developments in the functional methods. The QCD application of the Dyson- Schwinger equation has a 
long history and the mechanism of color confinement has been understood based on the gluon and the 
ghost propagators in the Landau gauge (see Ref. [365] for a review). Also, the Wetterich equation as we 



discussed in Sec. 4.2.2 is the exact formulation and has an equivalent content as the Dyson-Schwinger 
approach [366] . The QCD-based understanding of confinement is successfully transferred to the finite-T 
studies on the QCD phase transitions |223[ 13671 13681 369] . In these approaches toward QCD matter at 
high density, the most severe obstacle lies in the description of baryons and quarks on the equal footing. 
There is no model known in the market that can cope with baryons melting into quarks within a single 
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framework. In principle, such a unified treatment of baryons and quarks may be formulated by means 
of dynamical hadronization with quarks and diquarks. Further progresses in the Dyson-Schwinger and 
the functional RG methods with inclusion of diquarks are indispensable for future studies. 

The completion of the whole QCD phase diagram requires mutual collaborations extending over 
nuclear physics, high-energy particle physics and super-string theories, assisted by condensed-matter 
physics and astrophysics. This is such an interdisciplinary research area and new ideas are always 
desired. 
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